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ABSTRACT 


The objective of the integrated Digital Flight Control System is to provide 
rotational and translational control of the Space Shuttle orbiter in all phases of 
flight: from launch ascent through orbit to entry and touchdown, and during pow- 
ered horizontal flights. The program provides a versatile control system structure 
while maintaining uniform communications with other programs, sensors, and con- 
trol effectors by using an executive routine /functional subroutine format. The 
program reads all external variables at a single point, copies them into its 
dedicated storage, and then calls the required subroutines in the proper sequence. 
As a result, the flight control program is largely independent of other programs 
in the GN&C computer complex and is equally insensitive to the characteristics 
of the processor configuration. 

Section 1 of Part 1 describes the integrated structure of the control system 
and the DFCS executive routine which embodies that structure. Section 2 
deals with the input and output including jet selection. The succeeding sections 
show specific estimation and control algorithms for the various mission phases: 
cruise (including horizontal powered flight), entry, on-orbit, and boost. 

Part 2 describes attitude maneuver routines that interface with the DFCS. 



FOREWORD 


This document presents the Space Shuttle digital flight control design base 
developed by the Draper Laboratory Control and Flight Dynamics Division and 
the JSC Guidance and Control Systems Branch. The flight control and attitude 
control of the Shuttle from launch through touchdown and for powered horizontal 
flight is addressed. The subject matter reflects different levels of design 
maturity and completeness of documentation for the several different mission 
phases. It is anticipated that this report will be utilized to establish an initial 
baseline for the digital flight control software design. It is essentially the same 
as the NASA document Digital Flight Control Software Design Requirements, 

Rev. A . JSC-07759, July 6, 1973. 

This revision reflects more than just the additional design activities per- 
formed in the past half-year by the groups mentioned above. Honeywell Aerospace 
Division completed an extensive review of the original version, and some changes 
resulted from the constructive Honeywell review comments to NASA. Additionally, 
treatment of control modes in this revision reflects information exchanges with 
Rockwell and the Flight Crew personnel. 

There are numerous important design areas which have not been covered. 

It is anticipated that they will be addressed as baseline additions at an appropriate 
time. These areas include the treatment of redundant sensors and redundant 
actuators. Interface with pilot's displays has not been considered in this revision. 
Manual modes have been included only for the Cruise phase. All design to date 
has assumed rigid bodies; the filter structure, however, has been so specified 
as to facilitate the handling of bending and slosh when details are available. The 
Roll-out control phase has not been addressed. 

In this document, a number of areas have been addressed in considerable 
detail. The following paragraphs summarize the scope of the document and the 
design activity which it reports. 

The unified structure which is presented for the Digital Flight Control System 
(DFCS) is a useful, flexible framework for the evolving Space Shuttle control 
system. The implementation of input-output processing, on the other hand, can be 
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expected to change as the Shuttle avionics becomes better defined. 

Two designs for jet selection are included. The first, which is based upon 
linear programming, is flexible and powerful. The second, which is based upon 
a table look-up scheme, is more restricted. The latter will be considered the 
baseline design until the Shuttle requirements are found to go beyond its capabilities. 

The design for the Cruise phase addresses the requirement for a variety of 
manual, semi-automatic and automatic control modes during the developmental 
horizontal flight tests and the ferry flights as well as during the final stages of 
operational vertical flights. The control gains have been calculated for the 
Rockwell 0089B vehicle, but only for a single flight condition; the scheduling of 
gains for a full range of conditions is currently in progress. 

A number of features have been added to the design for the Entry phase 
since this document was originally issued. The applicable flight regime has been 
extended to subsonic, thus eliminating the separate Transition phase. A state 
estimator has been introduced. Rudder augmentation has been added to the control 
during the early part of entry. The gains have been calculated for the Rockwell 
0089B vehicle. A capability for higher maneuver rates in response to large bank 
commands has been added. The elevator-trim algorithm during early entry has 
been changed from open-loop to closed-loop. 

The design for control using the reaction control jets during the On-orbit 
mission phase has been developed and implemented, as has been the design for 
the attitude maneuver routine during the phase. The design does not reflect 
the availability of vernier jets or the unavailability of a clock for on-off commands 
other than once per minor cycle; these factors must be considered in redesigning 
firing policies. The design for control during major maneuvers in orbit (in- 
cluding Insertion) is in development; only a functional outline is presented in 
this document for that segment of the DFCS. 

The control for the Boost phase, which was developed outside the DFCS 
framework, is now incorporated into it. 

All the designs presented in this document have been set up for the Rockwell 
International 0089B Space Shuttle orbiter except that the ACPS thruster configuration 
is assumed to be that included in the Rockwell International ATP orbiter. 

There is some lack of commonality and uniformity in the material presented 
in this document as a result of various design efforts not yet being fully integrated. 
For example, the phase plane developed for on-orbit ACPS control may well be 
usable in the Entry phase as well. Such areas need to be addressed. 
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A version of the DFCS which incorporates the designs for the Cruise, Entry 
and Boost phases has been implemented in FORTRAN for testing on the SSFS 
(Space Shuttle Functional Simulator). A second version has been programmed 
in HAL using structured programming techniques. It incorporates the designs 
for the On-orbit phase, and runs on the engineering simulator (FSIM) which is 
being developed at the C.S. Draper Laboratory for the testing of flight programs 
that are written in HAL. 
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SECTION 1 


INTEGRATED STRUCTURE; EXECUTIVE ROUTINE 


by 


Peter S. Weissman 


1. 1 Integrated Structure 

The DFCS program consists of an executive routine, several major (functional) 
subroutines, and a number of utility subroutines. The executive routine performs 
"housekeeping" which is required regardless of the mission phase and calls the major 
subroutines in the appropriate order. There is a major subroutine for each of the 
following functions: Backup Control System (BCS) implementation, "pad load" manip- 
ulation, input manipulation, filter and parameter initialization, filter update (or meas- 
urement incorporation), control logic, command transmission, filter pushdown, (or 
state extrapolation), and parapieter estimation. Most of these routines are divided, 
at least in part, into separate segments for the various mission phases. 

This structure of executive routine and functional subroutines gives modularity 
to the program; it allows complete flexibility in designing the state estimation and 
control algorithms for each of the mission phases, designating new sensors or con- 
trol effectors, and indeed, redefining the mission phases themselves. At the same 
time, however, the unified structure which groups together the coding for similar 
functions encourages economy both in computer storage, and in engineering design 
and programming effort because it is so easy for different mission phases to use 
common coding whenever it is appropriate to do so. Another advantage to this struc- 
ture is that the calls to the functional subroutines indicate the completion of one set 
of calculations and the beginning of another; hence, they are suitable break points 
for external interrupts, and they facilitate the testing of new coding. For example, 
if modification to the orbit thrust vector control law is contemplated while other seg- 
ments of this phase's computations are unchanged, the new coding can be added as a 
totally separate subroutine, leaving the old coding intact. During testing, it is then 
possible to switch between old and new control laws following completion of the first 
segment of state filtering. This feature also could be of utmost importance for 
emergency or backup coding that is required to be functionally separated from the 
primary coding. 

The unified DFCS structure has additional, managerial advantages that are 
valuable in an enterprise as complex as the Space Shuttle program. The use in dif- 
ferent mission phases of the same assumptions about input/output, timing and other 
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interfaces is encouraged if not enforced. Furthermore, commonality in conventions, 
nomenclature and formats is encouraged, which simplifies and clarifies both testing 
and documentation. 

Not all operations of a control system need to be performed with the same 
frequency. The current version makes allowance for 3 sampling intervals, the long- 
er intervals being multiples of the shortest interval. In the current cruising flight 
phase, for example, the "fast" rate (10 sps) is used for angular control, the "medium" 
rate (2 sps) is used for velocity control, and the "slow" rate (1/2 sps) is used for 
altitude control and parameter updating. The medium and slow computations are 
offset, in order to prevent them from occurring on the same control cycle. Allow- 
ance is made for different sets of sampling intervals in different mission phases. 

1.2 DFCS Executive Routine 

An overview of the logic flow through the DFCS Executive Routine is shown in 
Fig. 1-1; Fig. 1-2 illustrates the implementation of this logic. The major subrou- 
tines of the DFCS are denoted in these figures by the hexagonal blocks. Although the 
executive serves all the mission phases, many of the algorithms within the subrou- 
tines are specifically designed for a single phase. Table 1-1 defines the mission 
phases for which the DFCS makes provision. The inputs to the Executive Routine 
are listed in Table 1-2; the outputs are listed in Table 1-3. The operations performed 
by the Executive Routine are discussed in the following paragraphs. 

The phase indication is changed from "entry" to "cruise" when the Mach number 
drops below 0. 9. All other determination of phase is currently assumed to be done 
outside the DFCS. 

When the Backup Control System (BCS) is requested, almost all the control 
system functions are performed using algorithms in a special subroutine that is sep- 
arate from the rest of the DFCS. The executive transfers control to this subroutine 
if the appropriate flags have been set. The BCS mode has been designed only for the 
cruise phase at this time. 

If there is an external restart request, or the phase or mode has changed since 
the previous pass, the executive sets flags to 

1) call for a DFCS initialization pass 

2) call for initialization of parameters for any modes that are new 

3) cause the computations that are done at the slower sample rates to be 
performed this pass. 

If the reason for initialization was an external restart request or a phase change, 
the Executive Routine establishes the counters and timers that reflect the desired 
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Fig. 1-1 Logic Flow of DFCS Executive 
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Performed 4 times with (x) 
replaced by U, E, A and R 
for the speed, elevator, 
aileron and rudder channels 
respectively. 



Fig. 1-2 Detailed Logic Flow of DFCS Executive 











Input Interface 
Subroutine 
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Control Output 
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Parameter Estimation 
Subroutine 




ISTART 
= 1 


Parameter Estimation 


Subroutine 


IM = IM + 1 modulo ITM 
IS = IS + 1 modulo ITS 


Return to External 
Scheduler 


Note: modulo operation results in the 
indices IM and IS being set to 
zero (which is the signal for the 
corresponding computations to 
be executed) when they reach the 
values ITM and ITS respectively. 


Fig. 1-2 Detailed Logic Flow of DFCS Executive (Cont) 


1-5 











Table 1-1 DFCS Mission Phases 


Rollout 

Cruise 


Entry 

On- orbit, 
TVC 

On- orb it, 
AC PS 

Insertion 

Boost 


Transition from aerodynamic control to steering and braking 
with landing gear; after touchdown; not currently implemented 

Approach and landing following entry phase (currently below 
Mach 0. 9) and horizontal powered flight (take-off through land 
ing). 


From after the deorbit maneuver to cruise phase; from total 
reliance on jets (ACPS), high angle of attack and hypersonic 
to primary control with the aerodynamic surfaces (ACS), low 
angle of attack and subsonic. 

Thrust vector control during maneuvers with the Orbi :al 
Maneuvering System. 

Control with the jets during coasting flight. 

Insertion into orbit following the Boost phase; not currently 
implemented, but similar to on-orbit TVC control. 

Flight from launch almost to insertion in orbit using the 
Solid Rocket Motors and/or the Main Propulsion System. 
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Table 1-2 Inputs to DFCS Executive Routine 

External request for DFCS initialization (1); none (0); 
cleared in Executive Routine 

Mach number 

MODA, MODE, Current control mode - aileron, elevator, rudder, 

MODR, MODU speed loops 

MODAG, MODEG, Current guidance mode - aileron, elevator, rudder, 
MODRG, MODUG speed loops 

MODAGR, MODEGR, Requested guidance mode - aileron, elevator, rudder, 
MODRGR, MODUGR speed loops 

MODAR, MODER, Requested control mode - aileron, elevator, rudder, 
MODRR, MODUR speed loops 

PHASER Flight phase flag (external) 

Table 1- 3 Outputs of DFCS Executive Routine 

IM, IS Medium- and slow- sampling- rate indices: 0 = perform 

calculations this pass 

ISTART Initialization pass (1), not initialization (0) 

ITM, ITS Number of fast passes per medium (slow) pass 

* MODB BCS mode (1), not BCS mode (0) 

MODACH, MODECH, Indicate mode change - aileron, elevator, rudder, 
MODRCH, MODUCH speed loops : 1 = no change; 2 = guidance mode change; 

3 = control mode change; 4 = phase initialization 

* PHASE Flight phase flag (internal copy of PHASER) 

Tp» Tjyp Tg Fast-, medium- and slow- sampling- rate sample 

periods for current phase 

Also an input to the next Executive Routine pass 


ISTRTR 

m 
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sampling rates for the current phase, and calls the major subroutine which transfers 
the set of pad- loaded constants required for the current phase from the pad- load lo- 
cations into working storage. The "pad-load" denotes those parameter values that 
are determined prior to a flight and are loaded into the flight computer on the launch 
pad . 

Regardless of whether or not it is an initialization pass, the Executive Routine 
calls the Input Interface Routine, which reads all the dynamic input required by the 
current phase and mode: sensor measurements, navigation estimates, guidance 
commands, and manual commands. All the dynamic data is read at the same point 
to ensure that all the values are registered with respect to the other values. The 
parameters are differenced, rotated into other coordinate systems and otherwise 
manipulated so as to be in the proper form for use in the other routines. During 
the appropriate phases, an Air Data Routine is called which manipulates some aero- 
dynamic parameters. These two routines are described in Section 2. 

At this point, the logic flow is detoured in an initialization pass to go through 
the Filter and Parameter Initialization Routine and the Parameter Estimation Routine. 
The former assigns values to those variables that need to be initialized when the 
phase (or mode) is entered. The latter determines the current values of parameters 
that will be regularly updated in the course of this phase (or mode). 

The flow then returns to the normal path. The subroutines entitled Filter - 
Part 1, Control Algorithms, Control Output and Filter - Part 2 are called in sequence. 
These subroutines comprise the core of the control system. After they are executed, 
the Parameter Estimation Routine is called; it is not called at this point in an initial- 
ization pass because it was executed earlier in the pass. The state filters are divided 
into two parts in order to minimize the amount of computation between data input and 
control command. The first part is limited to updating the filter estimate according 
to the most current measurement. The second part, which occurs after the control 
is exercised, is used to pushdown data (in the case of classical digital filters) or to 
propagate the state (in the case of discrete- time modern control estimation). The 
Control Algorithm Routine scales and combines the error signals from the filter rou- 
tine to form control commands for the aerodynamic surfaces, reaction control thrust- 
ers, etc. The Control Output Routine, which is described in Subsection 2.2, pre- 
pares the control commands for the DFCS interface with the control actuators. The 
Parameter Estimation Routine updates the gains required for estimation and control. 
The functions performed in this routine may include scheduled or adaptive gain adjust- 
ment. Much of the Parameter Estimation Routine is executed only at the slow sam- 
pling rate. 

After these routines have been executed, the Executive Routine increments the 
indices that indicate the "medium-" and "slow- sampling- rate" DFCS passes. If it is 
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the end of an initialization pass, the "slow- sampling- rate" index is given an extra 
increment. This extra increment will smooth the load on the computer by offsetting 
future "slow" passes from "medium" passes, provided that one sampling rate is an 
integer multiple of the other. The Executive Routine then returns control to the ex- 
ternal executive. 

The algorithms incorporated in the routines pertaining to input and output are 
presented in the following section of this document. The contents of the other rou- 
tines are presented phase by phase in the succeeding sections. 
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SECTION 2 


DFCS INPUT AND OUTPUT 

b y 

J. Edwin Jones, Peter S. Weissman, Craig C. Work, 

Edward T. Kubiak (JSC/EJ3) 

This section describes the interface between the DFCS and other components of 
the Space Shuttle avionics. Subsection 2. 1 discusses the DFCS inputs. It lists all the 
inputs; i. e. , state measurements, guidance commands, and pilot commands. It in- 
cludes a description of a routine which processes some of these quantities prior to 
their use in the filtering and control algorithms. It also describes a routine that fills 
the function of an air data computer. 

Subsection 2. 2 lists all the outputs from the DFCS to the control actuators and 
describes the limiting performed on these commands. Finally, Subsections 2. 3 and 
2. 4 present two candidate jet selection schemes; the first is powerful and flexible, but 
the second may suffice for the Space Shuttle requirements. 

Much of what is presented in this section depends heavily on many features of 
the Shuttle avionics that have not yet been specified, much less modelled in a simu- 
lator. Therefore, in many cases convenient assumptions have been made for the 
sake of proceeding with the DFCS design. These assumptions include the following: 

1) The inputs indicated in Subsection 2. 1 are, in fact, available to the DFCS. 

2) The actuators indicated in Subsection 2.2 are, in fact, available to the DFCS. 

3) The DFCS does not deal with redundancy management of sensors or actuators. 

4) The DFCS does not create outputs for pilot displays or telemetry. 

5) The DFCS can schedule software tasks to turn ACPS jets off and a 
sufficiently fine clock granularity is available for such scheduling. 

It is expected that each of these assumptions will, in time, be replaced by more 
realistic specifications, at which time these designs will have to be revised. 

2. 1 Inputs Required by the DFCS 

Table 2-1 lists all inputs currently required by the DFCS in the various mission 
phases. They are presumed to be available from sensors (or sensor- interface modules), 
the guidance and navigation module, and the pilot (or command- interface modules). 

In this interface, and throughout the DFCS, English units (viz feet, pounds, and seconds) 
and radians are used. 
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Table 2-1 Inputs to the DFCS 


Phase 

Mnemonic 

Symbol in Text 

Definition 

2,4 

ACC BOD 

5 b 

Linear accel measurement (body coord) 

2,4 

AG 

a G 

Angle-of -attack guidance command 

CD 

<T» 

AT V (1 ) 

0 

Platform outer gimbal angle 

00 

AI V (2) 

8 

Platform inner gimbal angle 

8, 6 

AIV (3) 

>P 

Platform middle gimbal angle 

8 • 

AI VC ( 1 ) 


Platform outer gimbal angle command 

8 

AI VC (2) 


Platform inner gimbal angle command 

8 

AIVC (3) 

*c 

Platform middle gimbal angle command 

2 

AY M 

ay 

Measured lateral acceleration 

2 

AZM 

Sz 

Measured normal acceleration (+zbody) 

4 

BANKG 

0 VG 

Bank guidance command 

2 

BG 

^G 

Sideslip guidance command 

2 

CG 

C G 

CAS guidance command (QG or CSTARG) 

2,4 

CRTB 

c b r 

Ref erence-to-body direction cosine matrix 

2 

DAG 

6» g 

• 

Aileron guidance command 

2 

DDG 

6 g 

TAEH guidance interface variable 

2 

DEG 

6e G 

Elevator guidance command 

8 

DELORB 

6 orb 

Prestored orbiter engine deflection command 

8 

DELSRM 

6 srm 

Prestored SBM engine deflection command 

2 

DBG 

6r G 

Rudder guidance command 

2 

DSG 

6s g 

Speed brake guidance command 

2 

DTG 

6t G 

Throttle guidance command 

2,4 

DYNP 

q 

Dynamic pressure estimate 

2 

HG 

h g 

Altitude guidance command 

2 

HM 

*sj 

H 

Measured altitude 

2 

IAGAIN 


AGAIN switch position (1 thru 4) 

2 

IATRIP! 


Roll stick trim activity: 0, * 1 (rgt) , - 1 (lef t) 
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Table 2- 1 Input to the DFCS (Cont) 


Phase Mnemonic 
2 IEGAIN 

2 IETRIH 

2 IGEAFX 

IN GUI D 
2 IPHM 

2 IRCHAI 

2 IRCHEI 

2 IRGAIN 

2 IRTRIM 

rSTRTP 
2 MODAGR 

2 HODAR 

2 HODBR 

2 HODEGR 

2 HODER 

2 HODRGR 

2 HODRR 

2 HODOGP 

2 MODUF. 

2 PG 

PHASER 


Symbol in Text 


INGUID 


ISTRTR 


P G 

PHASER 


2 PHSEL ^SEL 

2,4 PH pr 


Definition 

EGAIN switch position (1 thru 4) 

Pitch stick trim activity: 0, +1 (fwd) ,-1 (aft) 

Extend gear (1) , retract gear (0) 

Flag specifies guidance command set 
PHM valid (1), invalid (0) 

Reguest change IAINT state (1), no reguast(O) 
Reguest change IEINT state (1), no request (0) 
RGAIN switch position (1 thru 4) 

Pedal trim activity: 0,*1 (lft.) ,-1 (rgt) 

Request initialization (1 ) , no request (0) 

Requested values for HODAG 

Requested values for HODA 

Requested values for MODB 

Requested values for HODEG 

Requested values for HODE 

Requested values for MODRG 

Requested values for MODR 

Requested values for HODtIG 

Requested values for MODU 

Roll rate guidance command 

Flight phase flag (external) : 1 =rollout phase 
2=cruising phase 
4=entry phase 
5=orbit TVC phase 
6=orbit RCS phase 
7=insertion TVC phase 
8=booster TVC phase 
Pilot-selected values for PHREP 
Heasured roll rate 
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Table 2- 1 Input to the DFCS (Cont) 


Phase 

Mnemonic 

Symbol in Text 

2 

PSG 

>I>g 

2 

PS M 

* 

2 

PSSEL 

^SEL 

2 

QG 

q G 

2, a 

Q* 

q 

2,4 

RM 

r 

2,4 

ROLL 


2 

ROLLG 

0 bG 

4 

RR 

5rel 

=T 

* 

CM 

TEMPRK 

tk a 

2 

THG 

0 G 

2 

THM 

0 

2 

THSEL 

0 SEL 

6,8 

TNOW 


2 

UG 

U G 

2 

USEL 

U SEL 

2,4 

oxv,uiv,az V 

UV ,UV , 
“ x’ - y* 

4 

VBODY 


4 

VMACH1 


4 

VR 

V . 
— rel 

8 

MB (1,1) 

RRATEj^ 

8 

MB (2,1) 

PRATE jL 

8 

WB (3,1) 

YRATE t 

2 

X AM 

xa 

2 

XEM 

5e 

2 

XFM 

xf 

2 

XRM 

St 


Definition 

Heading guidance command 

Measured heading 

Pilot-selected values for PSREF 

Pitch rate guidance command 

Measured pitch rate 

Measured yaw rate 

Computed roll angle 

Guidance commanded roll angle 

Position vector in reference coordinates 

Input air temperature in deg. Kelvin 

Theta guidance command 

Measured theta 

Pilot-selected values for THEEF 
Present time 

Velocity guidance command 

Pilot-selected values for OREF 

Unit vectors along the velocity coordinate 

system, in reference coordinates 

Air-relative velocity in body coordinates 

Lateral switch parameter 
Air-relative velocity in reference 
coordinates 

Individual roll rate gyro signal 
Individual pitch rate gyro signal 
Individual yaw rate gyro signal 
Measured roll stick deflection 
Measured pitch stick deflection 
Measured flap controller deflection 
Measured pedal deflection 
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Table 2- 1 Input to the DFCS (Cont) 


Phase 

Mnemonic 

Symbol in Text 

Definition 

2 

XSM 

/v 

xs 

Measured speed brake controller deflection 

2 

XTH 

xt 

Measured throttle controller deflection 

8 

YACIN (I) 

«• 

Y. 

i 

Individual y accelerometer signal 

8 

ZACIN (I) 

• •' 
Z i 

Individual z accelerometer signal 


Table 2-2 Inputs to the Input Interface Routine 


Phase 

Mnemonic 

Symbol in Text 

Definition 

2,4 

ACC80D 

a* 3 

Linear accel measurement (body coord) 

2,4 

AG 

“G 

Angle-of-attack guidance command 

4 

BANKDO (2) 

a4 

Bank rate filter extrapolated value 

4 

BANKG 

0 VG 

Bank guidance command 

4 

CAT 

cosa^, 

Cos (AMEAN) 

2.4 

CRTB 

cj 

Reference-to-body direction cosine matrix 

2,4 

DYNP 

q 

Dynamic pressure estimate 

2 

HM 

H 

Measured altitude 

2,4 

PM 

<v 

P 

Measured roll rate 

2 

QG 

q G 

Pitch rate guidance command 

2,4 

QM 

3 

Measured pitch rate 

4 

RDOTC 

\c_ 

Bank rate command 

2,4 

RM 

r 

Measured yaw rate 

2,4 

ROLL 


Computed roll angle 

2 

ROLLG 

0 bG 

Guidance commanded roll angle 

4 

RR 

R , 
— rel 

Position vector in reference coordinates 

4 

SAT 

sinarp 

Sin (AMEAN) 

4 

VR 

—rel 

Air-relative velocity in reference coordinates 
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2. l. l Input Interface Routine 

The purpose of the Input Interface Routine is to preprocess measurements and 
compute state- estimator inputs. Inputs to the routine are listed in Table 2-2. Out- 
puts are listed in Table 2-3. Interface- routine operations are governed by the dis- 
cretes INGUID, ISTART, PHASE, and MODA in the manner summarized in Table 2-4. 
The remainder of this subsection is a description of the logic flow of the Input Inter- 
face Routine, which is illustrated in Fig. 2-1. 

The routine is called by the DFCS executive. If the phase is not either entry 
or cruise, it returns immediately. If it is a DFCS initialization pass, some variables 
internal to this routine are set to nominal values. 

The section of the routine for preprocessing of measurements is then entered. 

The angle of attack, the sideslip angle, and the air speed are computed using the 
air-relative velocity and the direction-cosine matrix. The Air Data Routine is called 
if the phase is in the aerodynamic regime. Some aerodynamic "measurements " are 
generated in the Air Data Routine, which is described in the next subsection. The 
computation of bank angle (under the indicated conditions) concludes the measurement 
preprocessing segment of the interface routine. 

The next section of the routine computes some nominal state values. Nominal 
values are computed for those state components of the entry phase state estimator 
(Section 3) for which nominal values are not specified by guidance. Nominal body rates 
are computed to provide turn coordination in the entry phase. A nominal (commanded) 
bank rate is also computed for use in the entry phase control logic in the specified 
flight regime d LAT = 0, q a q 3 ) when the bank error is sufficiently large. 

The filter input computation section is now entered. When in cruise phase, 
either the current roll angle or the bank angle measurement is put into the cell PHI 
for processing by the cruise phase filters. When in entry phase, the state estima- 
tor inputs (namely, errors in the angles of attack, sideslip and bank, and errors 
in the body rates) are computed. Additionally, estimator extrapolated values are 
adjusted for changes in the state-variable nominal values. 

Finally, a check is made to ensure that INGUID is negative. This avoids re- 
processing non-updated G&N inputs in subsequent DFCS passes. Control is then re- 
turned to the DFCS executive routine. 

2.1.2 Air Data Routine 

The Air Data Routine is called by the Input Interface Routine during the atmos- 
pheric flight phases. It computes several parameters required for aerodynamic 
control, serving the purpose of an air data computer. 


2-6 



Table 2- 3 Outputs from Input Interface Routine 


Phase 

Mnemonic 

Symbol in Text 

Definition 

4 

M (1) 

Ac? 

Angle-of-attack filter input 

2 

AM 

a 

Measured angle-of-attack 

4 

AO (2) 

Aa' 

Angle-of-attack filter extrapolated 

4 

BANK I (1) 

A0y 

Bank filter input 

4 

BI 0) 

At 

Sideslip angle filter input 

2 

BM 

e 

Measured sideslip 

2 

HM 

/V 

H 

Measured altitude 

2 

PHG 

0 G 

Bank angle guidance command 

2 

PHM 

0 

Measured bank angle 

4 

PI (1) 

Ap 

Roll rate filter input 

4 

Ql (1) 

A? 

Pitch rate filter input 

4 

HKD 

A? 

Yaw rate filter input 

2 

on 

tr 

Measured velocity 
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Table 2-4 Discretes Governing Logic Flow in Input Interface Routine 


ISTART 

DFCS initialization pass flag. 

= 1: initialize variables this pass. 
= 0: do not initialize this pass. 

PHASE 

Mission phase. 

= 2: cruise phase. 

= 4: entry phase. 

(Routine not entered in other phases. ) 

INGUID 

Flag from guidance. 

= i 2 ^guidance commands are pitch rate and 
roll (about body x-axis); terminal part of cruise 
phase. 

= ±1: guidance commands are angle-of-attack 
and bank (about velocity vector); entry and early 
cruise. 

< 0: as above, but no new guidance commands to 
process this pass (set negative internally). 

MODA 

Mode flag for cruise phase aileron channels. 

= 4;: automatic mode; compute roll (or bank) 
filter inputs. 

# 4: not automatic mode. 

*LAT 

Mode flag for entry phase lateral control. 

= 0: early entry, compute desired bank rate. 
= 1: late entry 


Note: The parameters , r^, r ( j,2‘ r <j> 3 are used in the I n P ut Interface 
Routine as limits. They are introduced in the discussion of the entry phase 
control logic in Section 5. 













Fig. 2-1 Input Interface Routine Logic Flow 
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Fig. 2-1 Input Interface Routine Logic Flow (Cont) 
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The inputs to this routine are the following measurements which are obtained 
from outside the DFCS or are computed by the Input Interface Routine: 


TK A - 

temperature in degrees Kelvin 

a - 

angle- of- attack (radians) 

H 

barometric altitude (ft) 

U 

air- relative velocity (ft /sec) 


The outputs are: 


O' rp 

stnc^ 

COScyn 


U 


M 

n 


mean (trim) angle- of- attack (radians) 

sine (cosine) of 

mean airspeed (ft/sec) 

Mach number 


The means (which are the outputs of first-order filters with time constants 
equal to the slow sampling period) and the sine and cosine are updated every pass. 

The means are set equal to the values of their respective measured inputs in an 
initialization pass. 

The other outputs are only computed on slow- sample- rate passes. The com- 
putations are shown in Fig. 2-2. 

It is expected that, in the future, when more realistic models of the sensors 
are available, this routine will compute the dynamic pressure and perform other 
computations. 

2. 2 Outputs to Control Effectors 

The DFCS computes commands for the control effectors listed in Table 
2-5. The question of how the commands are transmitted to these effectors is con- 
sidered to be outside the scope of this document; the following simple implementation 
is currently assumed: 

1) Commands are transmitted by the DFCS to the actuators (or a software 
module that interacts with the actuators) at the end of each execution of 
the control routine . 

2) The DFCS limits the magnitude of each output, but not its rate of change. 

3) Timing of jet firings is obtained by means of a computer clock that can 

be set to call a DFCS routine which will turn jets off after the required time. 
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Fig. 2-2 Computations in Air Data Routine 


2 - 








Table 2-5 Control Effectors Used by the DFCS 




Phases in 
Which Used * 

Aerodynamic Control System (ACS) 

elevons, rudder, speed brake, 
flaps 

1, 2, 4, 8 

Attitude Control Propulsion System (ACPS) 

40 jet thrusters 

4, 5, 6, 7 

Air Breathing Engine System (ABES) 

throttle 

1 , 2 (horizontal 
flight only) 

Solid Rocket Motors (SRM) 

c 

2 engines; gimballed at 45 from 
body axes 

8 

Main Propulsion System (MPS) 

3 engines; gimballed 

8 

Orbital Maneuvering System (OMS) 

2 engines; gimballed 

5, 7 


Phase designations defined in Table 1-1. 


Table 2-6 DFCS Limits for Commands to Control Effectors 


Effector Limit (deg ) 


aileron (differential elevon) 

6a max 

= 

± 15. 0 

elevator (tandem elevon) J 

’ R max 
1 6e 

, 6e min 


+ 15. 0 
- 40.0 

rudder 

. max 
6r 

= 

±15. 0 

speed brake j 

; . max 
| os 

i r min 
6s. 

= 

+ 70. 0 
- 17. 6 

flap ] 

1 

- „max 

l 6f . 

1 6f mm 

= 

+ 10. 0 
- 12. 0 

throttle 1 

j6t max 

= 

18, 000 lb 

1 

| 6t min 

= 

500 lb 

SRM and MPS gimbals 

5g max 

= 

± 8. 0 
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The assumed sign conventions for the aerodynamic control surfaces and the associated 
pilot controller motions are illustrated in Fig. 2-3. 

All the commands are limited within the DFCS prior to being issued. The 
limits magnitudes are shown in Table 2-6. The aileron and elevator commands are 
treated separately even though they apply to the same pair of actuators, the elevons, 
and are, therefore, interrelated; this question will be dealt with in the future. 

The limiting and blending (for the elevons and MPS gimbals) of the actuator 
commands is performed in the Command Output Routine, which is called by the DFCS 
executive after the Control Routine. The Command Output Routine serves as the 
single point from which the actuators (or the software modules that interact with the 
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RESPONSE 

PRIMARY FORCES 
AND MOMENTS 

PILOT CONTROLLER 
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Fig. 2-3 Sign Conventions for Aerodynamic Control Surfaces 
and Manual Controllers 
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2. 3 ACPS Jet Selection Subroutine 

This section describes a scheme for converting angular-rate-change and 
translation requests into timed jet firing policies for the ACPS (Attitude Control 
Propulsion System). This routine has been tested with the DFCS program. 

2.3.1 Thruster Geometry 

The thruster configuration currently incorporated in the jet selection routine 
is the North American Rockwell ATP design. 

As shown in Fig. 2-4, the thrusters for the ATP configuration are grouped 
in 16 locations. There are 40 thrusters in all, and these are distributed among 
the 16 locations as shown in Tables 2-7 and 2-8. There are 16 thrusters in the 
nose of the fuselage and 24 in the "tail". 

2. 3. 2 Jet Selection Logic 

The DFCS Control Routine determines a commanded velocity change and calls 
the jet selection subroutine, which handles angular velocity components and trans- 
lational velocity components simultaneously with a single algorithm. A flag can be 
set instructing the jet selection subroutine to ignore translation, simplifying and 
speeding up the jet selection procedure. Furthermore, the jet selection logic can 
be instructed to deal with jet clusters or with the individual thrusters. 

Figure 2-5 is a functional flow diagram of the jet selection logic. A dis- 
cussion of this logic follows. 

The ACPS jet control authority subroutine is only fully executed when the 
jet selection routine is called for the first time in a mission phase. The control 
authority subroutine calculates the ACPS control authorities on the basis of total 
vehicle mass, the vehicle inertia matrix, and the center of mass location. The 
control authorities are calculated for the jet clusters and for the individual thrust- 
ers. The control authority is the angular and translational acceleration in vehicle 
coordinates generated by firing a jet or cluster of jets. In addition, a cost per 
unit firing time is assigned to each jet. Currently, all these costs are equal. 

When a jet is disabled (or restored) and the jet- selection logic is operating on 
clusters, the jet control authority subroutine is partially executed, appropriately 
altering the corresponding acceleration vector and cost element. 

The jet selection subroutine next checks its inputs to see whether or not 
there is a new velocity change request. If not, the jet selection subroutine call 
must be for the purpose of turning off jets which have completed their firing 
intervals. If the velocity change request is a new one, the linear programming 
subroutine is called, which expresses the velocity change vector as a linear 
combination of jet acceleration vectors, the linear coefficients being the firing times. 
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Fig. 2-4 Thruster Locations - NR ATP Baseline Orbiter 


Table 2-7 Distribution of ACPS Thrusters 


NR ATP Baseline Orbiter 


Thruster Location Number of Thrusters Thruster ID 


1 

2 

10 , 

12 



2 

2 

14, 

16 



3 

2 

13, 

15 



4 

2 

9, 

11 



5 

2 

6, 

8 



6 

2 

5, 

7 



7 

4 


19, 

21, 

23 

8 

4 

18, 

20, 

22, 

24 

9 

3 

32, 

34, 

36 


10 

3 

26, 

28, 

30 


11 

3 

25, 

27, 

29 


12 

3 

31, 

33, 

35 


13 

2 

2, 

4 



14 

2 

1 , 

3 



15 

2 

38, 

40 



16 

2 

37, 

39 




40 total 
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Table 2-8 Parameters of Individual Jets 


ID 

Jet Location 

— Body Station (inches) 

Thrust 

Comp 

onents 


Fabrication Frame 


Vehicle 


Frame 

1 

318. 0 

- 96.0 

374. 0 

-1.0 

0.0 

0 . 0 

2 

318. 0 

96. 0 

374. 0 

-1.0 

0 . 0 

0 . 0 

3 

318. 0 

- 96. 0 

349. 0 

-1.0 

0.0 

0 . 0 

4 

318. 0 

96. 0 

349. 0 

-1. 0 

0 . 0 

0 . 0 

5 

318. 0 

-109. 0 

368. 0 

0.0 

1.0 

0 . 0 

6 

318. 0 

109. 0 

368. 0 

0 . 0 

-1.0 

0 . 0 

7 

318. 0 

-109. 0 

355. 0 

0 . 0 

1.0 

0 . 0 

8 

318. 0 

109. 0 

355. 0 

0 . 0 

-1.0 

0 . 0 

9 

318. 0 

-109. 0 

346. 0 

0 . 0 

0 . 0 

-1. 0 

10 

318. 0 

109. 0 

346. 0 

0 . 0 

0 . 0 

-1. 0 

11 

330. 0 

-109. 0 

346. 0 

0 . 0 

0 . 0 

-1. 0 

12 

330. 0 

109. 0 

346. 0 

0 . 0 

0 . 0 

-1. 0 

13 

318. 0 

- 80. 0 

367. 0 

0.0 

0 . 0 

1. o 

14 

318. 0 

80. 0 

367. 0 

0 . 0 

0 . 0 

1. 0 

15 

330. 0 

- 80.0 

367. 0 

0 . 0 

0.0 

1. o 

16 

330. 0 

80. 0 

367. 0 

0.0 

0.0 

1. o 

17 

1536. 5 

-158. 0 

425. 0 

0.0 

1. o 

0 . 0 

18 

1536. 5 

158. 0 

425. 0 

0 . 0 

-1. o 

0 . 0 

19 

1532. 5 

-158. 0 

438. 0 

0 . 0 

1. o 

0 . 0 

20 

1532. 5 

158. 0 

438. 0 

0 . 0 

-1.0 

0 . 0 

21 

1528. 5 

-158. 0 

450. 0 

0.0 

1. o 

0 . 0 

22 

1528. 5 

158. 0 

450. 0 

0.0 

-1.0 

0 . 0 

23 

1525. 5 

-158. 0 

464. 0 

0.0 

1.0 

0 . 0 

24 

1525. 5 

158.0 

464. 0 

0.0 

-1. o 

0 . 0 

25 

1541. 5 

-130. 0 

480. 0 

-0. 234 

0.0 

0. 980 

26 

1541. 5 

130. 0 

480. 0 

-0. 234 

0 . 0 

0. 980 

27 

1555. 5 

-130. 0 

484. 0 

-0. 234 

0 . 0 

0. 980 

28 

1555. 5 

130. 0 

484. 0 

-0. 234 

0 . 0 

0. 980 

29 

1567. 5 

-130. 0 

487. 0 

-0.234 

0 . 0 

0. 980 

30 

1567. 5 

130. 0 

487. 0 

-0. 234 

0 . 0 

0. 980 

31 

1546. 5 

-130. 0 

416. 0 

0. 407 

0 . 0 

-0. 915 

32 

1546. 5 

130. 0 

416. 0 

0.407 

0 . 0 

-0. 915 

33 

1559. 5 

-130. 0 

420. 0 

0. 407 

0 . 0 

-0. 915 

34 

1559. 5 

130. 0 

420. 0 

0.407 

0 . 0 

-0. 915 

35 

1572. 5 

-130. 0 

428. 0 

0. 407 

0 . 0 

-0. 915 

36 

1572. 5 

130. 0 

428. 0 

0. 407 

0 . 0 

-0. 915 

37 

1579. 5 

-149. 9 

468. 0 

0. 980 

0.0 

0. 234 

38 

1579. 5 

149.0 

468. 0 

0. 980 

0.0 

0. 234 

39 

1579. 5 

-149. 0 

455. 0 

0. 980 

0 . 0 

0. 234 

40 

1579. 5 

149. 0 

455. 0 

0. 980 

0. 0 

0. 234 

c. g. 

1096.5 

0. 1 

367. 9 
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Fig. 2-5 ACPS Jet Selection Functional Flow Diagram 
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Generally, many different jet combinations could satisfy a given velocity change 
request; if there is more than one such combination, the total cost is minimized 
in the selection process. Any jet which is flagged because of failure or mission- 
related constraints is omitted from consideration. 

The selected jets are turned on (if they are not already firing), and all others 
are turned off and removed from the internal jet- off list. The interrupt clock is 
set to return control to the jet selection subroutine when the next jet- off is due, 
and the "new velocity request" flag is zeroed. 

The interface parameters for the jet selection package are: 


INPUT: 


a) 

Velocity change command vector 


b) 

Flag for new velocity vector. 


c) 

AC PS jet exclusion list. 


d) 

Translation inclusion flag. 


e) 

Phase flag. 

OUTPUT: 


a) 

Jet-on/off commands. 


b) 

Interrupt clock settings. 


c) 

Zeroed flag for new request. 


Figure 2-6a presents features of the linear programming subroutine, CCWLJS, 
and Fig. 2-6 b describes its logic flow. 
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Linear Programming Subroutine 

PURPOSE : Represents input vector W as a linear combination of linearly 

independent columns from matrix B, choosing columns from B 
to minimize cost of solution. 

INPUTS : LROW - Number of components in W. (Dimension of requests) 

LCOL - Number of columns in B matrix. (Number of jets or clusters) 

W - Vector (size LROW) to be analyzed as a linear combination of 
columns from B. 

COST - Array (size LCOL) of strictly positive cost coefficients, one 
associated with each B column. 

B - Matrix (LROW by LCOL). Each column represents the accelera- 
tion associated with a one-second firing of the corresponding 
jet in vehicle coordinates. Components may include rotational 
or translational acceleration elements, or both. 

LALARM - Non-zero indicates that previous jet selection may not be taken 
as a starting basis, for one of the following reasons: 

1) Jet exclusion list has been changed 

2) 1$, LROW, or LCOL have been altered since the preceding 
pass 

3) Preceding delta- V request was not successfully resolved into 
jet firings. Initialized to 1 on first pass. 

JFAIL - Array (size LCOL) of flags, J th element non-zero to indicate J th 
jet failed. 

# — — 

BI, CF, LQ - Can be inputs as well as outputs. See definitions below under 
OUTPUTS. 

PARAMETERS - SAVING, MAXPAS, MAXITR, TOLRNC, URGNCY can be chosen 
to match the precision of the computation and the accuracy of the 
data. 


Fig. 2-6 a Features of Linear Programming Subroutine, CCWLJS 
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Linear Programming Subroutine (Cont) 

OUTPUTS : BI - Transformation matrix (LROW by LROW) going from vehicle 

coordinates to coordinates based on selected columns of B. 

V - Solution vector (size LROW)_of coefficients for selected columns 
of B matrix, satisfying: Bl W = V 

CF - Vector (size LROW) of cost coefficients selected from COST 
array, corresponding to chosen columns from B matrix. 

Z - Cost of solution, satisfying: Z = CF • V. 

LQ - Array (size LROW) giving the indices and sequence of selected 
B columns. 

MCOUNT - Number of basis element substitutions performed in current 
execution. 

LALARM - Non-zero if W was not successfully resolved into jet firings before 
MCOUNT exceeded MAXITR . 

ASSUMPTIONS: 1) Cost array entries are all strictly positive, i.e., greater 

than zero. 

2) Negative LQLIST outputs and negative V components are to be 
excluded from any implementation of the solution. 

3) COST array entries are small compared to URGNCY. 

4) BI matrix is never altered by any other program. 

5) LALARM will be input as non-zero whenever changes occur 
in LROW, LCOL, B, or the jet exclusion list. 

ENTRY : Call CCWLJS (B, W) 

EXIT: RETURN to caller. 


Fig. 2-6 a Features of Linear Programming Subroutine, CCWLJS (Cont) 
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For I = 1 to LROW, then go to 840 

* 

i. - Bid, i) 


-I - LQd) 



MCOUNT counts iterations on this pass 


If BI cumulative roundoff is not bad, 
if previous solution was valid, and 
if S, LROW, LCOL are unchanged, 
then start with solution from previous 
pass. Otherwise initialize completely. 

If W is unchanged from previous pass, 
use same solution. 


NPA3S counts iterations since last 
complete initialization. 


Negative LQ values identify vehicle 
frame coordinates, which are non- .jets. 


For I = 1 to LCOL, then go to 100 

I 


I - LRCORD(I) 



Positive LRCORD values correspond 
to non- selected jets. 


Fig. 2- 6b Flow Chart of Linear Programming Subroutine, CCWLJS 
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V contains latest version of solution 
in terms of current basis. 




© 



Cost target vector is_initially a 
weighted version of V (or a signum 
(V) vector). 


URGNCY factor places premium 
on jets which approximate W or 
which fire in directions represented 
so far by non- jets. 


Y(I) is weighted cost of generating 
the Ith original (vehicle) basis 
vector as a linear combination 
of current basis vectors. 

Initialize the search for the next 
jet to be included in basis. 


Fig. 2-6b Flow Chart of Linear Programming Subroutine, CCWLJS (Cont) 
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I 



Search for next jet to be 
introduced to basis, bypassing 
any which are failed or already 
in basis. 


X Is the weighted price penalty 
for simulating a one-sec firing 
of jet, using the basis rather 
than the jet itself. 


Greatest saving per second is 
H; corresponding jet is K. 


If no jet was selected for 
inclusion, then exit. 


K iet is to be included. Y 
is K>h jet acceleration vector 
described in current basis. 


Fig. 2-6b Flow Chart of Linear Programming Subroutine, CCWLJS (Cont) 
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For 



I = 1 to LROW, then go to 330 



Choose J basis member for 
replacement so that the new 
V(J) is £ 0. , and no other V 
component changes sign. 


Fig. 2-6b Flow Chart of Linear Programming Subroutine, CCWLJS (Cont) 
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Bad result. 


Fig. 2-6b Flow Chart of Linear Programming Subroutine, CCWLJS (Cont) 
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450 ) 


EXIT 



Calculate _cost of solution: 
Z = CF • V 

Save W for comparison 
In next pass. 


Negative firing times represent 
requests which cannot be met 
with currently non-failed jets. 


NPASS reflects quality of BI 
matrix. 


Return to caller. 


Fig. 2-6b Flow Chart of Linear Programming Subroutine, CCWLJS (Cont) 
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2.4 ACPS Jet Selection Logic for the ATP Configuration 


A jet selection routine is presented in this section that is more restricted 
than the design of the preceding section. It is essentially a table look-up; 
the firing policy for each possible command, as well as the individual jets 
and the control authority corresponding to each jet cluster are pre-determined 
and stored for access by this on-line routine. The North American Rockwell 
ATP design for the ACPS, which has been presented in Fig. 2-4 and Tables 
2-7 and 2-8. was assumed in the off-line computations. 

The flow chart of Fig. 2- 7a describes in detail the jet selection logic 
for the ATP configuration. Figure 2-7b defines the flowchart symbols. A 
discussion of some of the features of the logic follows. 

General - To provide attitude control with no translational disturbances, jet 
couples must be used. However, pure jet couples exist in only one axis, the 
roll axis. Jet couples to be used for pitch control (e. g. 3 & 1 2 ) and yaw 
control (e. g. 5 & 7), also effect roll. To achieve pure pitch or yaw two 
thruster couples must be fired to null the roll disturbance; for example, 3 and 1 
12 with 2 and 10 for pitch and 5 and 7 with 3 and 1 for yaw. Hence, couple 
control can require as many as four thrusters to fire whereas non- couple 
control may require only one thruster firing (e. g. yaw control). Additionally, 
the multiple thruster firings required for couple control produces large minimum 
impulse rate changes and these cause fast limit cycles and inefficient propellent 
usage. This logic has the option of attitude control with either coupled or 
non- coupled jets. 

Nose versus Tail Thrusters - for Non- coupled Jet Operation - The logic allows 
the crew to select either nose or tail thrusters for control of each axis. 

This serves two purposes: 

1) Comparing nose thrusters to tail thrusters, the nose thrusters have 
shorter roll moment arms, but longer pitch and yaw moment arms. Hence, 
the crew can select the preferable moment arm length for each axis. 

2) This is a simple means of accounting for thruster failures. Should 
all the thrusters in any single location fail, control can still be maintained 
(Note: Individual failures in any jet location are accounted for by placing another 
jet of the same location in "common" with the location, see page 16 of the 
flowcharts). 

Phase Plane Jet Selection Logic Interaction - (Non- coupled Operation) - 
Thrusters in locations 1 through 4, and 9 through 12, produce both pitch and 
roll torques. Hence, it takes two thrusters to provide either pure roll or pure 
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pitch. However, smaller minimum impulse rate changes are possible by firing 
only one thruster. Therefore, even if there is a phase plane command from 
only one axis (roll or pitch), the logic will attempt to do the latter (fire only 
one thruster) if the resulting disturbance in the "non- command" axis reduces 
the rate error (called RR^ or RRy >- n the logic). 

Use of the tail thrusters for yaw control causes a roll disturbance; should 
the roll disturbance exceed a rate limit, the logic will consider the disturbance 
part of the phase plane command and cause corrective roll thrusters to fire. 

Translations - No combination of thrusters will produce Y or Z translations 
with negligible rotational disturbances. Thrusters must be modulated to null 
rotational disturbances. This modulation could be performed either open 
loop in the jet selection logic, or closed loop through the phase plane logic 
(the crossing of a switching line) and jet selection logic. The latter method 
was chosen because it will not produce unnecessary modulation (thruster 
firings) during small Y or Z translations (i. e. , when rotational disturbance 
is not long enough to cause switching line crossing). 

During +X- translations, it is possible to control yaw by modulating the + X- 
thrusters; this option is available in the logic. 

Jet Coupling and Logic Flow - With respect to thruster control actions, there 
are two basic types of thrusters: 

1) Those that produce Y translation and yaw. 

2) Those that produce Z translation and roll and pitch. 

Sequentially, the logic first calculates commands for Y translation and yaw and 
then issues commands for Z translation, roll, and pitch control. 


2-31 




Fig. 2- 7a Flow Chart Jet Selection Logic for ATP Configuration 
(Flow Chart continues for 22 pages) 
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THCX,Y,Z 

^ C’®C'^C 
a i, 0,©, or tf 

J,K 

MI 

RDFY 

BRL 

HRy 

PI 

ST 

M,N 

t CM 

E 

*C 

ROC 

RRLT 

ST 1 

**1 

CL 


Commands for translation, either positive, none, or 
negative, depending on +1, 0, or -1 

Commanded rate changes from switching logic. 

1H 

Angular Acceleration for i jet for either roll, 
pitch, or yaw. 

Value of 1 (not 0) signifies request for long moment 
arm thrusters 

Indices for jets 

Minimum impulse time 

Roll rate disturbance for yaw and/or Y translation firings 
Roll rate limit before RLFY is added directly to 0 
Roll rate error from phase plane logic when 0 p is zero 
Index indicating presence (l) or absence (0) of roll command 
Sample time 
Jet indices 

Maximum roll compensation time available by unbalancing 
pitch firing times 

Small number to account for round off in comparisions 
(say. 10 _i +) 

One jet roll correction time 
Net roll command, i.e., RR + 0 C 

Roll rate limit before RR^is added to 0 C during translation 
Sample time plus E 

Pitch rate error from phase plane logic when © is zero 
Couple option 


Fig. 2- 7b Definition of Symbols for Jet Selection Logic 
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SECTION 3 


STATE ESTIMATOR (ENTRY PHASE) 


by 


J. Edwin Jones 


3 . 1 Introduction 

The function of the state estimator is to obtain estimates of controller state 
variables in an uncertain environment. This uncertainty arises from noise in the 
measurement process and noise driving the states directly, e.g., due to wind gusts. 

The state estimator design must meet certain prescribed requirements in the 
face of externally-imposed constraints. These requirements, many of which are 
not yet specified quantitatively, pertain to the degree of accuracy and smoothness 
in the estimated state variables, the time delay between reading the sensors and 
issuing computed control commands, and the effective phase lag associated with 
filter dynamics . External constraints are imposed by fixed sensor hardware and 
configuration characteristics and computer timing and sizing limitations. 

Meeting all requirements while satisfying all constraints may necessitate 
design compromises. The design should be comprehensive to ensure prescribed 
accuracy by making best use of available information, but design sophistication 
is limited by computer constraints. Noisy inputs may require heavy attenuation to 
ensure prescribed smoothness, but the low gains associated with heavy attenuation 
must not cause excessive phase lag or excessive sensitivity to modeling errors. 
These compromises are central in the development of the design presented below. 

The current state estimator is designed specifically for the entry flight phase 
(and some of the parameters used by the estimator are defined in the entry control 
description. Section 5). The structure used, however, is applicable to any flight 
phase. Extension of the state estimator to other flight phases is under study. 

3 . 2 General Description 

The design presented for the state estimator is linear and therefore implies 
approximation, namely linearization of the state dynamics about a set of operating 
points and linear modeling of control effector and sensor operation. The linear 
assumptions facilitate both design and analysis. Nonlinear effects will be noted as 
development continues and possibly included. 
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The design is cast in the state-space formulation and uses more than one 
measurement to produce each state-variable estimate (multi-input, multi- output) . 

It is comprehensive in the sense that it uses all pertinent information. However, 
approximations are made to decrease computer requirements; namely, only a 
reduced number of state -variables are estimated, and decoupling of longitudinal 
and lateral dynamics is assumed. An important feature of the design is that control 
effects are fed into the estimator's state -extrapolation step, keeping the estimator 
response fast even while noisy inputs are being heavily attenuated. 

Time delay between reading sensors and issuing control commands is mini- 
mized in the sense that it provides for incorporating the current measurements and 
extrapolating the state t in two separate steps with the latter performed after the 
commands are generated. Furthermore, the slow-varying filter parameters are 
updated separately at a slower repetition rate. Details of these computer timing 
considerations are given in Implementation, Subsection 3.3. 

The estimator gains are constant with respect to time at any particular flight 
condition; however, their on-board computation is dependent on the current estimates 
of dynamic pressure and Mach number. When the dynamic pressure is below a 
prescribed level, tf the estimator gains are set equal to predetermined fixed values. 
When dynamic pressure is above this level, the gains are interpolated from pre- 
determined tabulated values with Mach number as the sole independent variable. 

The fixed gain values used at very low dynamic pressure were determined 
classically on the basis of desired filter bandwidth and damping ratio . m The 

j}c >}c 

scheduled gain values used for high (above q^ or q^) dynamic pressure were found 
by analogy to the steady-state Kalman (Wiener) filter for the reduced, decoupled, 
linearized state assuming locally stationary noise statistics at each operating point. 
Measurement noise strengths were chosen to reflect relative sensor quality; driving 
noise strengths were chosen to keep the gains at a selected level, sufficiently low . 
to achieve desired smoothness but sufficiently high to exclude extreme sensitivity 
to modeling errors. 


* Analogous to "filter update" and "filter pushdown" in classical filters. 

tt The prescribed dynamic pressure level isq^ for the longitudinal channel and q<> 
for the lateral channel. Current values for q* and q| are 20 and 30 lb/ft^, 
respectively. 

Ttt Currently, the low-dynamic-pressure estimator bandwidth is 5 rps and the 
damping ratio is 0.7 sec" 1 . 
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3. 3 Implementation 


3.3.1 Input-Output 

Estimator state variables (outputs) are linearized variations about nominal or 

A 

commanded states. These state variables are the error angle-of-attack An and its 

A A ^ 

time derivative Ad, the error sideslip angle A0 and its time derivative A/ 3, and the 
error bank angle Ad y and its time derivative Ai p^. (The symbol [ A ] signifies an 
estimated quantity.) 

To obtain variational outputs, the measurements must be differenced with 
their respective desired values to form variational inputs. These variational inputs 
are the measured errors in angle-of-attack A a, body pitch rate Aq, sideslip angle 
A/3, bank angle Ad , body roll rate Ap, and body yaw rate Ar. (The tilda [~] 
symbolizes a measured quantity.) The measurements are either read directly from 
sensors or computed (preprocessed) from the sensor data. The body rate measure- 
ments are assumed to be directly accessible from body- mounted rate gyros, but 
the angle measurements are derived from IMU (Inertial Measurement Unit) and 
body-mounted accelerometer information. Commanded values corresponding to each 
of the measurements are either received from the guidance subsystem or computed 
internally. A detailed description of estimator input computation is presented in 
Subsection 2. 1. 

3.3.2 Measurement Incorporation 

Current measurements are incorporated into the state estimate immediately 
after they are obtained and preprocessed. This is immediately before the control 
subroutine is called, in order that control commands computed using the most recent 
information may be issued as soon as possible after the DFCS is called. Subsequently, 
the state estimate is extrapolated in preparation for the next DFCS pass, as shown 
in Subsection 3.3.3 below. 

Measurements are incorporated according to the equations 


A/3 


Ad 



+ K 


Ion 


(A a - A a). 


(Aq - A a ) 


r ~ 

(A/S - A(S ) 


(3-1) 








(Ap - A/3 sin - Ad y cos a,p) 


~ S' ' ^ ' 

(Ar + A/3 coS - Ad y sin a^) 




(3-2) 
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where the prime (') indicates a state variable propagated from the previous sample 
instant, is the 2x2 longitudinal gain matrix, and is the 4x4 lateral gain 

matrix. The gain matrices are updated at the slow repetition rate in the parameter 
estimation routine, which is performed after the state extrapolation step in the 
normal DFCS cycle. 

The flow charts in Fig. 3-1 and 3-2 show the implementation in the DFCS 
program of the measurement incorporation equations and the computation of the 
estimator gain matrices. When dynamic pressure is below q. for the longitudinal 
channel and q 2 for the lateral channel, the constant values KC1 and KC2 are loaded 
into the estimator gain matrices as shown; the current values for these parameters 
are 0.45 and 2.5 respectively. For dynamic pressure greater than q* and q* for 
longitudinal and lateral channels respectively, the estimator gain matrices are found 
by linear interpolation on Mach number; the values are given in Tables 3-1 and 3-2. 

3.3.3 State Extrapolation 

After control commands are computed using the current state estimate (from 
Eq (3-1) and (3-2)) and issued, the estimator state is extrapolated to the next sample 
according to the equations: 



r *\ 


r 

a / 

A/3 


/% 

A/3 

“f ' 


A 

A/3 


A/3 

A ' 

> ■ < 

> 

A 

A <t> 

V 


A0 v 

aJ; 


-> 

A <t> , 

V. J 

k+1 

v 

V- J 


k 


0 

(sin o^.) 
0 


(cos a T ) 


0 



(sin » T ) 


(3-3) 


(3-4) 


CL CL 

where and <t>^ ^ are longitudinal (2x2) and lateral (4x4) state transition 

matrices*; Ap^, Aqj, Arl are the average changes in body angular rates about the 
x, y and z axes due to ACPS jet activity over the previous sample interval, and 
« T is the mean angle-of-attack. 


♦ cl 

The transition matrices are designated <J> (i.e., closed-loop) to indicate that 

closure of the aerodynamic control loop is taken into account in their derivation: 


.CL . OL , „aero ,-iaero 
cj> = cj> + r c 

where is the open-loop state transition matrix, r aer0 is the aerodynamic- 
control-to-state matrix, and c aer0 are the linear-feedback aerodynamic control 
gains . 
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Fig. 3-2 Computation of Estimator Gain Matrices 
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Table 3-1 


Longitudinal Estimator Gains Scheduled vs Mach Number 



Mach = 0.9 

Mach = 2 

Mach = 5 

Mach =10 

Mach =20 

(1. 1) 

0.0731 

0.0769 

0.0680 

0.0568 

0.0401 

<1, 2) 

0.0408 

0.0411 

0.0554 

0.0523 

0.0586 

(2, 1) 

0.0538 

0.0353 

0.0329 

0.0207 

0.0116 

(2,2) 

0.4320 

0.4289 

0.2351 

0.2616 

0. 1909 


Table 3-2 Lateral Estimator Gains Scheduled vs Mach Number 



Mach =0.9 

Mach = 2 

Mach = 5 

Mach =10 

Mach = 20 

(1. 1) 

0.2119 

0.0578 

0.0427 

0.0314 

0.0222 

(1,2) 

-0.0297 

-0.0794 

-0.1495 

-0.1098 

-0.0942 

(1,3) 

-0.0115 

-0.0282 

-0.0234 

-0.0210 

-0.0203 

(1.4) 

-0.0337 

-0.0456 

-0.0637 

-0.0605 

-0.0659 

(2, 1) 

0.2410 

v 0.0238 

0.0214 

0.0115 

0.0065 

(2, 2) 

0.0257 

0.0155 

0.0016 

0.0177 

0.0075 

(2,3) 

0.0298 

0.0696 

0.0701 

0.1256 

0.0900 

(2, 4) 

-0.4127 

-0.4213 

-0. 2223 

-0.2261 

-0. 1650 

(3, 1) 

-0.0027 

-0.0008 

-0.0015 

-0.0011 

-0.0009 

(3, 2) 

0.4924 

0.4372 

0.3569 

0.3213 

0.2440 

(3,3) 

0.0148 

0.0187 

0.0198 

0.0224 

0.0256 

(3,4) 

0.0010 

0.0030 

0.0072 

0.0122 

0.0142 

(4, 1) 

-0.1164 

-0.0212 

-0.0181 

-0.0133 

-0.0070 

(4, 2) 

1.355 

1.130 

0.8698 

0.6984 

0.4032 

(4,3) 

0.3703 

0.3846 

0. 1853 

0.2043 

0.1473 

(4, 4) 

0.0387 

0.0725 

0.0783 

0.1229 

0.0903 
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where the parameters ce d and C d specify the desired longitudinal response. (See Sub- 
section 5.3. 1. 1). Equation (3-7) is used for dynamic pressure below q^, when the 
elevator is essentially used only for trim. Equation (3-8) is used for dynamic pres- 
sure above q ^ . 

CL 

The lateral state transition matrix is found by Eq (3-5) from F^ at , which takes 
on either of two forms: 



3-8 



where the parameters a^ through b ^ specify the desired lateral response 
(see Subsection 5. 3. 2. 1. 2). Equation (3-9) is used for dynamic pressure below 
when the aileron is essentially not used. Equation (3-10) is used for dynamic 
pressure above . 

The flow charts in Figs. 3-3 and 3-4 show the implementation in the DFCS 
program of the state extrapolation equations and the computation of the state transi- 
tion matrices. The numerical values of the parameters used in the state extrapola- 
tion are indicated in Table 3r3. The control gains C'. , C' , C'. c' will 

6a 6 a p 

be obtained from the entry control module in the future. Currently, however, they 
are independently scheduled as functions of Mach number, as shown in Fig. 3-5. 



Fig. 3-3 State Extrapolation for the Entry State Estimator 
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Fig. 3-4 Computation of Estimator State Transition Matrices 
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Fig. 3=4 Computation of Estimator State Transition Matrices (Cont) 
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Table 3-3 Parameter Values for State Extrapolation 
Aerodynamic and Inertial Parameters 


s 

2690. 

b 

78. 056 

I 

XX 

805000 

I 

zz 

6067000 


Fixed Control Gains 


a 2 

a 3 

a 4 

a ll 

a 12 

a 21 





0. 2401 
6 . 0 
2. 94 
2. 9008 


See Table 5-1 


Scheduled Control Gains 


See Figs. 5-16 through 5-21 





See Fig. 3-5 
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Fig. 3 5c C * vs Mach Number 
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3. 4 Preliminary Performance Summary 


Simulated maneuvers have been performed on the Space Shuttle Functional Sim- 
ulator using the state estimator in conjunction with the entry DFCS and the NR00089B 
(February version) vehicle. Both large and small step changes in angle of attack and 
bank were commanded under a wide range of flight conditions (Mach = 2, q = 150 lbs /ft 
to Mach = 25, q = 0. 25 lbs/ft^). 

Comparisons were made between simulations in which the DFCS controller oper- 
ated on state estimator outputs (FILTER = 1) and corresponding simulations in which 
the DFCS controller operated directly on measured data (FILTER = 0). When random 
noise sources were not simulated (perfect measurements), comparisons between 
corresponding simulations were very good in every case. When random noise sources 
were simulated, corresponding simulations exhibited similar mean behavior, but aero- 
dynamic surface commands were smoother for those cases which used the estimator. 
Figure 3-6, for example, shows actual angle of attack, roll angle, and aerodynamic 

surface position time histories for corresponding cases at Mach 14 and a dynamic pres- 
2 

sure of 108 lbs/ft . 
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Fig. 3-6b Longitudinal and Lateral Responses with Noise Sources 
Simulated and FILTER = 1 (estimator in loop) 
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SECTION 4 


CRUISE PHASE 


by 


Albert G. Engel, Franklin H. Moss, Alexander Penchuk, Nicole Lamotte 


4. 1 Introduction 

The cruising flight phase of an orbital Space Shuttle flight begins at about 
Mach 0. 9, ends at touchdown, and follows the high-angle-of-attack/high-Mach-number 
entry phase. During the cruise phase the control is accomplished by moving the aero- 
dynamic control surfaces: elevons, rudder, rudder flare (speed brake), and body 
flap. In addition, the cruise phase DFCS has been extended to deal with Shuttle pow- 
ered horizontal flight, which encompasses developmental testing and ferry missions. 
During powered flights, thrust is an additional control variable. 

For the cruise phase the DFCS contains three primary modes of operation: 
manual, semiautomatic and automatic. Manual mode consists of two types of control 
systems: Stability Augmentation System (SAS) and Command Augmentation System 
(CAS). By means of a side arm controller and rudder pedals, the pilot usesjeither of 
these two systems to manually control the orbiter. The pilot can reduce his workload 
while controlling the orbiter by selecting certain automatic functions provided by outer 
loops. This is the semiautomatic mode of operation. In the automatic mode of opera- 
tion guidance algorithms provide the commands to the digital flight control system, 
and the pilot's control role is reduced to monitoring the overall performance. One 
further system available is the Backup Control System (BCS) which is nearly identical 
to the SAS. This system is completely manual, and following its initial flight test, 
is intended to remain unchanged throughout the Space Shuttle program. 

Speed control is accomplished by using either speed brake or throttle. The 
settings can be controlled in an open loop fashion either by the pilot in the manual 
mode or by the guidance algorithms in the automatic mode. Alternatively the settings 
can be obtained from closed loop speed control loops. In this case, depending on the 
mode of operation, the speed commands come either from the pilot (select/hold mode) 
or from the guidance. For the purpose of pitch trimming, the body flap's setting is 
scheduled as a function of Mach number, which can be supplemented by a manual flap 
command. 

Stabilization of the vehicle in the cruise phase is accomplished with compensa- 
tion filters which appear in forward loop, feedback, and crossfeed paths. It was 
found that for control of the rigid body many of these filters are not needed and good 
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control performance can be obtained with simpler compensations than allowed for. 
However, these filters will become necessary and may, in fact, become insuificient 
in size as control of the flexible vehicle is considered. A table of gains suitable for 
control of the rigid body 89B vehicle at a particular flight condition (Mach 0.6, 
dynamic pressure 300 psf, altitude 6000 ft) is provided (Table 4-1). The gains which 
must be varied as a function of some flight parameter (such as dynamic pressure), 
are marked by an asterisk. 

The description of the cruise phase autopilot is presented in the order of 
increasing detail, starting with a discussion of simplified block diagrams, proceed- 
ing to detailed functional block diagrams, and then finally to detailed flow charts. 
Subsection 4.2 describes manual, semiautomatic, and automatic control modes as 
well as speed and body flap controls. The description of the BCS is followed by that 

of mode selection and mode change trim initialization. Filter implementation is 

th 

discussed in Subsection 4.3 where a generalized n -order digital filter implementa- 
tion is detailed along with low -pass and washout filters, and integration functions. 
Subsection 4.4 describes preflight initialization of gains and filter coefficients, as 
well as limits, deadbands and other items. Subsection 4. 5 presents summaries of 
input/ output quantities for various channels such as elevator, aileron, etc. The 
detailed functional block diagrams in Subsection 4. 6 are preceded by a list of conven- 
tions in order to facilitate their presentation. Subsection 4. 7 contains functional 
flow charts for initialization subroutine of the cruise DFCS, and then detailed flow 
charts are presented in the Appendix, Subsection 4A. 

4. 2 Control Modes 

4.2.1 Manual Modes of Operation 

For manual operation the two types of control systems are a command augmenta- 
tion system (CAS) for normal operation, and a stability augmentation system (SAS) for 
lower performance downmode operation. The pitch rate and roll rate feedback loops, 
with simple gain compensations in the feedback paths, are employed in SAS (Fig. 4-1). 
The manual commands to the elevator and aileron are summed with these feedback 
signals and sent out to the control surfaces. The manual rudder command is summed 
with two feedbacks and one crossfeed signal from the aileron before it is sent to the 
rudder. Washed-out yaw rate feedback augments the dutch-roll stability. The turn 
coordination in SAS is achieved by either a combination of aileron to rudder crossfeed 
and lateral acceleration, or sideslip feedback. For selecting the best control response, 
pitch rate, roll rate, and yaw rate gains can be scaled by the pilot from his control 
panel. 

A simplified block diagram of CAS is shown in Fig. 4-2. The longitudinal and 
lateral control systems are pitch rate (or C )loop and roll rate loop closures having sim 
ilar forward path compensation filters followed by a proportional path and an integral 
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* Gain- switch- sensitive gains 


Fig. 4-1 Simplified Functional Block Diagram of The Stability Augmentation System 
(also Backup Control System) Mode 
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Fig. 4-2 Simplified Functional Block Diagram of the Command Augmentation System Mode 










bypass. The integral bypass provides automatic trimming and assures unit-gain rate 
tracking. This function can be switched in or out from the control panel, and the pitch 
rate/C loop has an inner damping loop and an outer feedback loop enabling a selection of 
either C or pitch rate. C is a blend of pitch rate and normal acceleration. The 
pitch rate feedback can be fed back directly, passed first through a washout filter, or 
fed back with an added steady turn compensating component. This component is pro- 
portional to the yaw rate and tangent of the roll angle which is measured or derived 
from the measured roll rate. Depending on the type of feedback (C or pitch rate), 
the manual command is scaled according to either C or q for the control variable. 

The roll rate loop does not have an inner damping loop and thus uses roll rate feed- 
back. The rudder channel is the same as in the SAS. 

4. 2. 2 Semiautomatic Modes of Operation 

Certain selectable outer loop closures provide semiautomatic modes of opera- 
tion, the construction of which is shown in Fig. 4-3. For longitudinal control the in- 
puts are either pitch attitude or altitude references with the controlled variables being 
pitch attitude and altitude. For lateral control the references and the controlled 
quantities are the bank angle and heading. In both cases, engaging the control stick 
opens these outer loops and enables the control stick steering function. Therefore, 
returning the stick to null reengages the selected semiautomatic mode once the body 
rates fall below a threshold value. 


4.2.3 Automatic Modes of Operation 


For one particular cruise phase guidance scheme the automatic mode is sub- 
divided into two regions of operation according to the following table: 


Name of 
guidance phase 


Approximate Guidance 

altitude range commands 


Approach 

Terminal 


40000 to 7000 feet 
7000 feet to touchdown 


a G' 

q G' 


0 

^ V 


G 


0 


b 


G 


During the approach phase the longitudinal control is of the angle-of-attack. In the 
terminal phase, longitudinal control is accomplished by the pitch rate loop which is 
identical to the CAS loop described in Subsection4. 2. 1 (including pitch rate compensa- 
tion during steady turns). In both cases the lateral control is of the roll attitude. 
However, the roll angle commands are different in the two guidance cases: roll 
angle is taken about the velocity vector in the approach phase, and about the body 
(vehicle x-axis) in the terminal phase. For both cases the form of the lateral control 
system and the gains are the same. 

As illustrated in Fig. 4-4, guidance commands other than those shown above 
can also be accommodated. Thus the longitudinal control system can be configured 
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Fig. 4-3 Simplified Functional Block Diagram of Outer Loops for Semiautomatic Modes 
of Operation 






4 ^ 
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Fig. 4-4 Simplified Functional Block Diagram of the Automatic Mode 









as either an altitude loop or a pitch attitude loop, both of which are outer loops closed 
around the inner CAS loops. Similar to the CAS, an additional longitudinal control 
mode has C feedback. As well as the roll angle loop already discussed, the lateral 
control system is provided with a heading loop closed around the roll attitude loop. 

The rudder loop is implemented similar to the SAS or CAS mode, except that 
a provision is made to accept sideslip guidance commands and direct or washed-out 
yaw rate feedback for automatic decrab maneuver. Rudder pedal commands enable 
this maneuver to be made manually. 

It should be mentioned that, although not shown on the simplified block diagram, 
additional guidance commands can be applied in the automatic mode. These are 
elevator, aileron, and rudder commands applied directly to the appropriate summing 
junctions, and roll rate guidance command to the roll rate summing junction. 

4.2.4 Speed Control 

Control of indicated airspeed for the cruise phase is accomplished by two control 
effectors, the speed brake (also commonly known as "rudder flare"), and the air 
breathing engine system (ABES). Use of the ABES is anticipated only for purposes 
of powered horizontal flight operations, but the appropriate control algorithms re- 
sulting in throttle command are instituted in the current comprehensive cruise DFCS. 
Simplified functional block diagrams of the speed control modes are presented in 
Fig. 4-5. 

In BCS and manual modes, which are functionally identical, the pilot's manual 
commands are passed directly to the control actuators. Although not indicated in 
Fig. 4-5, speed brake and throttle commands in BCS and manual modes, as in all 
speed modes, are crossfed to the elevator for the purpose of trimming out resulting 
pitching moment. Note that BCS and manual control may employ speed brake or 
throttle or both. 

The IAS-hold mode is functionally similar to the automatic guidance mode which 
provides IAS commands, and is a velocity loop closure with forward loop compensa- 
tion. Currently, the speed brake loop is implemented with proportional plus-integral 
compensation, and the throttle loop has lead-integral lag compensation. In this mode 
either speed brake or throttle is utilized and the nominal elevator mode selection is CAS. 

The automatic (direct guidance) mode issues commands to the speed brake 
and throttle based directly on guidance issued speed brake and throttle commands. 

This may occur for either speed brake or throttle or both. 

Finally, the automatic Terminal Area Energy Management (TAEM) guidance 
mode receives the TAEM guidance rate command and first passes through switching 
logic to decide whether the command is to be issued to speed brake or throttle. If 
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BCS AND MANUAL MODES 
(Speed Brake or Throttle or both) 


Manual . | ) 

Speed Brake /' . — . 

Command " ! 

I I 

Manual I i 

Throttle - 

Command ] I 


Speed Brake 
Command 


Throttle 

Command 


IAS-HOLD OR AUTOMATIC MODES 
(Speed Brake or Throttle) 


IAS or 
Guidance 
Velocity 
Command 


IAS or 
Guidance 
Velocity 
Command 



AUTOMATIC (DIRECT GUIDANCE) MODES 
(Speed Brake or Throttle or both) 


Guidance | 

Speed Brake S— 

Command I 

Guidance j 

Throttle - 

Command ( 


I Speed Brake 

j Command 

I ^ Throttle 
j Command 


AUTOMATIC -TAEM GUIDANCE 


TAEM 

Guidance 

Command 



Speed Brake 
Command 


Throttle 

Command 


Fig. 4-5 Simplified Functional Block Diagrams - Speed Control 
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an indicated airspeed decrease is called for and the throttle is presently opened 
above its minimum, the command is issued to decrease throttle, or otherwise, the 
speed brake is opened. On the other hand, if an indicated airspeed increase is 
called for and the speed brake is presently opened above minimum, the command is 
issued to close the speed brake, or otherwise, the throttle is opened. These com- 
mands are scaled appropriately, integrated, and issued to the controller. 

4.2.5 Body Flap 

Control of the body flap, an aft pitching moment generator, is achieved via 
commands issued according to a piecewise linear body flap versus Mach number 
profile, biased with a direct pilot manual command (see Fig. 4-6). In addition, 
appropriate command biases are instituted to provide for full control surface deflec- 
tion range. 

4.2.6 Backup Control System (BCS) 

The BCS mode is an autonomous, digitally implemented, multi-axis manual 
control capability which serves as pilot requested backup in the event of primary 
DFCS failure. Although the design presented here is only of a preliminary nature, 
it is intended to demonstrate the capabilities currently envisioned for BCS mode 
plus techniques for their implementation. 

Whereas the primary DFCS is likely to be subject to modifications throughout 
the flight testing and orbital flight operations, the intent for BCS is to obtain an 
acceptable design prior to the first flight and to leave it unmodified thereafter. As 
the only BCS coding which is integrated into the primary DFCS structure is the BCS 
call flag from the executive routine, possibility of subsequent faulty modifications is 
virtually eliminated, thus providing the pilot with a continuously reliable alternative. 

BCS mode is functionally similar to primary manual mode and shares the same 
simplified functional block diagrams (Figs. 4- 1, 4-5, and 4-6). However, there are several 
fundamental differences which are not indicated: The axis-by-axis, multi-routine 
format of the primary DFCS is replaced by a single routine which processes all 
axes in succession; and the filter update, control and filter pushdown (wrap-up) 
functions are imbedded in the proper sequence. Whereas the primary system will 
have some gains scheduled as a function of dynamic pressure, a single set of fixed 
gains which provides suitable performance over as large a cruise flight spectrum 
as possible will be utilized for BCS. Finally, pilot signal preprocessing is not in- 
cluded for elevator, aileron, and rudder, and there is no Mach number profile for 
body flap. 
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Fig. 4-6 Functional Block Diagram of Body Flap 


Body-Flap 

Command 
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4.2.7 Mode Selection 


Mode select /indicate interfaces with the pilot are helpful in understanding 
the capabilities of the cruise phase DFCS, especially where manual control and a 
multiplicity of modes complicate the picture. Figure 4- 7a is compatible with simul- 
taneous multi-axis mode selection in contrast to axis-by-axis mode selection. Thus 
it is compatible with earlier versions of the cruise phase DFCS wherein axis-by-axis 
mode changes were not considered. Multi-axis mode selection is not unlikely for the 
ultimate cruise phase DFCS version, considering such things as control loop design 
failure modes, pilot training, and computer input /output simplifications. The gain- 
select switches and integrator moding select /indicators are on a separate panel. 

This panel might be associated with the ABES console since these capabilities are 
intended primarily for the flight test phases of development . 

Figure 4- 7b, on the other hand, suggests an axis-by-axis mode selection 
capability compatible with the present cruise phase DFCS design, including all 
peripheral moding capability. All indicated functions are implemented, except for 
the pilot selectable pitch turn- coordination- override capability. Gain- select switches 
are active only in the SAS mode. Pilot -selectable integrator- moding- override capa- 
bility is provided only in the CAS mode (pitch and/or roll). 


4. 3 Digital Filters and Their Implementation 

aj. 

4.3.1 Generalized n -Order Filters 

th 

Control filters for the cruising mode are classical, recursive, n -order 
digital filters, whose Z transform is D(Z): 


D(Z) = 


a l + a 2 Z + 


+ a 2n Z 


-n 


1 + a 3 Z _1 + 


+ a 2n+l Z 


The general form of their implementation is presented in block diagram format 
in Fig. 4-8. The symbol x represents the controlled variable and the a's are either 
fixed or scheduled gains (currently fixed) which correspond to the gains of form 
"Gx" which appear in the flow charts and block diagrams presented later in this 
section. 
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Fig. 4-7a Mode Select/Indicate Panels, Type 1, Cruise Phase DFCS 




Fig. 4- 7b Mode Select/Indicate Panels, Type 2 (Axis-by-Axis) Cruise Phase DFCS 
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nd 

The current filter implementation provides for no greater than 2 -order filter 

capability, but extension to higher order filters (if needed) is straightforward. The 

Z * calculation blocks are achieved in the DFCS programming by the relative sequence 

of instructions, which permits the DFCS cycle time lag to occur at the point so desig- 

nd 

nated. In the coding, each 2 -order state filter is divided into two parts in order to 
minimize the computational transport lag between data input and control command. 

The first part, which occurs in the Filter Update Routine is limited to updating the 
filter estimate according to the most current measurement: 

x , = a, x. + x. T1 
out 1 in N1 

Where x Nn represents the n^-node, or summing junction, of the filter. Note that 
exactly n nodes are required to implement an n^-order filter. 

The Filter Update Routine is sequentially followed by the Control Routine which 
issues control commands based upon the most current filter output value x Filter 
data "pushdown", which is actually the tabulation of the current value associated with 
each node, is then accomplished in the Filter Pushdown Routine by the following two 
equations for a 2nd-order filter: 

x. T , = a„x. - a„ x , + x„ T „ 

N1 2 in 3 out N2 

x. TO = a . x. - a c x , 

N2 4 in 5 out 

This sequence of instructions results in the following equivalent difference 
equation for the filter output: 

x . = a, x. +a„x. + a . x. - a 0 x , -a c x , 

out 1 m 2 in_ 1 4 m_ 2 3 out^ 5 out_ 2 

Where subscripts -1 and -2 denote past values, one and two samples previous. 


Control loops which utilize 2 nc *-order filter capability are available for control 
of the following quantities: 


BCS 

Manual 

Modes 

Automatic and 
Semi-Automatic 
Modes 

Feedback Quantities 

Others 

roll rate 
pitch rate 
yaw rate 

roll rate 
pitch rate 
yaw rate 

roll angle 
pitch angle 
yaw angle 

angle-of-attack 

sideslip angle 

indicated air- 
speed 

altitude 

measured pitch rate 

measured yaw rate 

measured lateral 
acceleration; 

normal acceleration 
(with gravity) 

turn coordination 
rudder crossfeed 
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4.3.2 Specialized 1 st - Order Filters 
st 

Recursive, 1 -order washout filters are utilized in the pitch rate and yaw 
rate feedback paths of the primary system and in the yaw rate feedback path of BCS. 
The Z transform of the filter is: 

1 - Z _1 

D(Z) = — — — — ry 

1 + G ,„Z 
xW 

A block diagram and the update and pushdown equations for the washout filter 
are presented below 



Fig. 4-9 Digital Washout Filter 

Filter Update: = x in + x^ 

Filter Pushdown: x NJ = -x. n - G xW x out 


s t 

Recursive, 1 -order low pass filters are utilized for the filtering of the fol- 
lowing quantities: 

1) pilot stick inputs 

2) crossfeeds to elevator (from speed brake, throttle, and body flap) 

3) elevator trim tracker 

4) yaw rate for pitch turn coordination 

The Z transform of the low pass filter is 

1 + G 

D(Z) = 

1 + G Z 1 

X 
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and the block diagram and update and pushdown equations are presented below 



Fig. 4-10 Digital Low Pass Filter 

Filter Update: x out = x in + G x x^ 

Filter Pushdown: x^ = -x Qut + x. n 


4.3.3 Integrators 

Integration is achieved by both trapezoidal and rectangular integrators. Trape- 
zoidal integrators are employed for the integral bypass functions in the elevator and 
aileron loops, while rectangular integrators are utilized for functions in which trape- 
zoidal integration proves to be of no additional advantage. Thus, rectangular inte- 
grators are utilized for the integration of pilot trim discretes (elevator, aileron, and 
rudder in primary system and BCS) and for the integration of TAEM guidance rate 
commands (speed brake and throttle). The Z transforms of these integrators are 

X X “ 1 

Rectangular: D(Z) = — r Trapezoidal: D(Z) = - 1 + - 

l-z ' 1 2 x-z " 1 


and the block diagrams and update and pushdown equations are as follows: 

limiter limiter 



Rectangular Integrator Trapezoidal Integrator 


Fig. 4-11 Digital Integrator Filters 

Rectangular Integrator Trapezoidal Integrator 



Rectangular 

Trapezoidal 

T 

F 

Filter Update: 

x out " x out + T F x in 

x out 2 x in X N1 

T 

F 

Filter Pushdown: 

none 

x ivri = x - + X , 

N1 2 m out 
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Note from the block diagrams that the output of the integrators is limited; in 
practice this occurs in the filter update routine. Also notice that rectangular inte- 
gration requires no pushdown equation, thus substantiating its use whenever possible. 

4 . 4 Preflight Parameter Specifications 

Basic parameters for the cruise phase DFCS include filter coefficients, gains, 
limits, deadbands, logic-control flags, etc. Each of these basic parameters must 
be numerically specified prior to the cruise phase initialization pass. A representa- 
tive preflight parameter specification is presented alphabetically: Table 4- la for 
non-subscripted parameters and Table 4- lb for subscripted parameters. All numeri- 
cal data is based upon a single-point design at Mach 0.6 and 300 psf (about 6000 feet) 
utilizing aerodynamic data for the February 1973 AER016 version of the 89B orbiter. 
The RCAH version of the CAS elevator mode is indicated (ICSTAR = 0); the rudder 
sideslip feedback option is chosen (IBFDBK = 1); pitch turn coordination is bypassed 
(IQTRN = 1) . 

Most of the basic parameters are utilized directly, but a number of "working" 
quantities are employed, such as sample- interval- multiplied trim rates or a - guidance 
values of CAS filter coefficients, limits and loop gains. These are computed on the 
initialization first pass (Fig. 4-24) or upon mode initialization as required (Fig. 4-25 
through 4-28). 

Earlier orbiter cruise phase DFCS designs have utilized scheduled gains and 
coefficients to provide desirable response characteristics over the cruise phase 
flight envelope (including Powered Horizontal Flight). The scheduling parameter 
has been dynamic pressure. The scheduling operations (when enabled) are carried 
out every slow pass by the Parameter Estimation Routine (Fig. 4-12), which operates 
on preflight- specified tables of parameter value vs dynamic pressure. Typically 
seven data points are used: 100, 150, 225, 300, 425, 525, and 650 psf. The thirteen 
parameters corresponding to the scheduled quantities of earlier designs are identi- 
fied by asterisks in Table 4-1 . Tabular specifications for these gains and coefficients 
will be added to the preflight parameter specification tables, as they become avail- 
able . 

4 . 5 Input / Output Summaries 

Figures 4-13a through 4-13e detail those quantities which are input to and 
output by the five basic control channels (speed, elevator and body flap, aileron, 
rudder and BCS). Input includes information transferred across the various inter- 
face boundaries — executive, guidance, and pilot, as well as miscellaneous moding 
information, crossfeeds from other loops, and sensor measurements. Output con- 
tains control effector commands, pertinent mode status for display, and crossfeeds 
to other loops. 
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Table 4-la Preflight Parameter Specifications for the Cruise Phase 

{Non- subscripted Parameters) 


AINTMX 

.2400 

ETf-M 

-.1745 

GF SB 

. 0 . 

CCAMN 

-.0873 

ETMX 

.0654 

GET 

0. 

CCAMX 

.3491 

* GAINT 

.707 

GETB 

0. 

CCCMN 

-.0873 

* GAP 

.288 

GETF 

-.9900 

CCCMX 

.3491 

GAPHB 

.288 

GEX 

.2327 

OACMX 

.4800 

GAX 

.2400 

GEXB 

.2327 

DECMN 

-.6981 

GAXB 

.2400 

GF X 

.7679 

DECMX 

.2618 

GCFOM 

12.422 

GF XB 

.7679 

OFCMN 

-.3840 

GCH 

.3 

* GPHPS 

4. 

DFC.MNR 

-.3840 

■* GCTh 

1 .4 

* GPPH 

1.19 

ORCMX 

.2618 

GCX 

-.1163 

GPX 

.1745 

OSCMN 

-.307? 

GOFC 

-.9048 

GOW 

-.9802 

OSCHNH 

-.307? 

GO SC 

-.3679 

GRAY 

0 . 

DSC MX 

1.2217 

GOTC 

-.9753 

* GRB 

.314 

D5UMX 

.1257 

* GFCA 

-.5 

GRCPH 

0 . 

DTCMN 

500 

* GECC 

-1 . 

GRDAC 

0 . 

OTCMNB 

600 

GEF 

0 . 

GRDACB 

0 . 

OTCMX 

1 HOOD . 

GFFB 

0 . 

GRM 

-.9802 

DTD MX 

H 7 5 0 . 

* GEINTA 

.4 

GRRE 

0 . 

DTDMXH 

8750. 

* GE1NTC 

.707 

* G«RWA 

1.19 

F I NT MM 

-.3491 

* GFQMA 

.2 

GRRWB 

l.)9 

F 1 N T to X 

. 1309 

GFOMR 

1 . 

GRW 

-.980? 

FTDB 

.0087 

* GFOtoC 

• 0 . 

GRWB 

-.980? 

ETDMX 

.0174 

GFS 

0 . 

GRXAMT 

.2094 


Scheduled Gains 
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Table 4-la Preflight Parameter specifications for the Cruise Phase (Cont) 


(Non-subscripted Parameters) 


GRXB 

.2094 

ICSTAR 

0 

X F MX 

.5406 

GRXSAS 

.2094 

I PHRC 

1 

XETMN 

1 

o 

GSDDG 

1. 

1QTRN 

1 

xetmnr 

-0.75 

GSU 

1. 

IRWP 

. 1 

XETMX 

0.75 

CSX 

2.03B5 

PCM X 

.3491 

XETMXR 

0.75 

GSXB 

2. 0385 

PHCMX 

.5236 

XETR 

0.075 

GT DOG 

152.25 

RMOtS 

.0873 

XETRB 

0.3 

GTU 

H. 3X10000 

OCTMX 

.0367 

XFMMN 

0 . 

GT X 

35000 

TFB 

• i 

xfmmnr 

0. 

GT XB 

35000 

X A R 

0. 

XRB 

0 . 

GXA 

-.6065 

XADB 

.0? 

XRDB 

.02 

GXAX 

1. 

XAOMX 

.5 

XROMX 

.25 

. GXE 

-.6065 

XAMDIS 

.05 

XRMX 

1. 

GXEX 

2. 

X A MX 

1. 

XRTMX 

.3125 

GXR 

-.6065 

XATMX 

.5 

XR TMXR 

.3125 

GXRX 

.25 

XATMXR 

.5 

XRTR 

.0156 

GXSX 

.25 

XATR 

.05 

XRTRB 

.0625 

gxsxb 

.25 

. XATRR 

.2 

XSMDIS 

. 1 

GXTX 

.5 

XEB 

0. 

XSMMN 

0. 

GXTXR 

.5 

xerjh 

.02 

XTMDI S 

. 1 

I ABES 

1 

XFDMX 

.3 

XTMMN 

0 . 

I ANGRV 

3 

XErtDI S 

.05 



I BFDBK 

1 

XEMN 






* 

Scheduled Gains 



Table 4- lb Preflight Parameter Specifications for the Cruise Phase 

(Subscripted Parameters) 

OFCMCHU.J) ( —#2094 * .1), (-.2094, .2), (-.2094, .5), 

(-.0873, .52), (-.0873, 6), (.1745, 8), (.1745, 9) 

GAN ( I ) 1., 0., 0., 0., 0. 

GAPMS(I) .144, .2304, .288, .3456 


GAYMI 1 ) 

1 • t 

0. 

, 0. 

. 0. 

, 0. 

GBE ( I ) 

1 • t 

0. 

t 0. 

. 0. 

. 0. 

GC EA ( I ) 

1 • f 

0. 

♦ 0 . 

. 0. 

. 0. 

GCEC ( I ) 

1., 

0* 

, 0. 

. 0. 

. 0. 

GOACI I ) 

1. , 

0. 

* 0. 

. 0. 

. 0. 

GEOMSI I ) 

.5, 

.8 

. 1. 

. 1 . 

2 

GHE(I) 

1., 

-. 

9905 

. .01. 0., 0. 

GPE ( I ) 

1., 

0. 

, 0. 

. 0. 

. 0. 

GPHE ( n 

1. , 

0. 

, 0. 

. 0. 

, 0. 

GPHT ( I ) 

. 1 , 

0. 

♦ ”• 

9y9 , 

0 • t 0 • 

GPM ( I ) 

1., 

0. 

, 0. 

* 0. 

. 0. 

GPMB ( I ) 

1 . , 

0. 

. 0. 

. 0. 

. 0. 

GPSE ( I ) 

1 • * 

0. 

1 0, 

. 0. 

. 0. 

C , QM ( I ) 

1., 

0. 

, 0. 

. 0. 

A 0. 

GOMB( I ) 

1. , 

0. 

, 0. 

. 0. 

, 0. 

GRE ( I ) 

1., 

0. 

♦ 0. 

. 0. 

A 0. 

GREB ( I ) 

1. , 

0. 

. 0. 

A 0. 

A 0. 

GRRWS ( I ) 

.595, 

.9S2 

. 1 . 

19 A 1.4PR 

GSIJE ( I ) 

1. , 

-. 

93, 

-} A A 

0 A t 0 • 

GTHE ( I ) 

1., 

0. 

, 0. 

. 0. 

A 0. 


GTtlE(I) 1., 


-1.9713, -1.905, .9716, .9048 




Fig. 4-12 Parameter Estimation Routine for Cruise Phase 
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velocity guidance command (Uq) 
speed brake guidance command (fis^) 
throttle guidance command (St^) 
TAEM guidance command ( 5stj_, ) 


mode change indicator (MODUCH) 
BCS mode indicator (MODB) 


speed mode requests (MODUR) 
manual speed brake controller (xs) 
manual throttle controller (xt) 
speed selection (Ug) 


air breathing engine 

system indication (IABES) 

none 

measured earth-relative 
velocity (U) 



speed brake command (Ss^) 
throttle command (St^) 

speed loop mode (MODU) 
speed guidance mode (MODUG) 


6s c and 6t^ to elevator loops 


Fig. 4-13a Speed Control Input/Output Interface Variables 
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elevator guidance mode requests (MODEG) 
altitude guidance command (Hq) 
pitch attitude guidance command (0 q) 
angle of attack guidance command (o'q) 

CAS guidance command (c^ = or c *) 
elevator guidance command 

pitch mode change indicator (MODECH) 
BCS mode indicator (MODB) 

pitch mode requests (MODER) 

control stick steering override (ECSS) 

measured pitch stick & body flap controller (xe, $£f) 

pitch trim indicator (IETRIM) 

pitch attitude selection (0g) 

pitch integrator control (IEINOK) 

pitch SAS gain select indicator (IEGAIN) 


<fi - validity indicator (IPHM) 
gear extension indicator (IGEARX) 

speed brake command (Os^-,) 

throttle-command (Otj-,) 

bank angle for pitch turn coord (^) 

measured pitch rate (qf) 
measured pitch attitude (d) 
measured altitude (H) 
measured angle of attack (a) 
measured z-body acceleration (az) 
measured yaw rate (r) 
measured bank angle i$) 
measured Mach number (W) 



elevator command (be^,) 
body flap command (Of^) 

elevator loop mode (MODE) 
elevator guidance mode (MODEG) 

Ge^-, to elevon commands 
bank angle for turn coordination 
(^ T ) to rudder loop 


Fig. 4- 13b Elevator and Flap Input/Output Interface Variables 
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aileron guidance mode requests (MODAG) ] 
heading guidance command (4 *q) 
bank angle guidance command ($q) 
roll rate guidance command (Pq) 
aileron guidance command (Sa^-.) 


GUIDANCE 


aileron mode change indicator (MODACH) \ EXECUTIVE 
BCS mode indicator (MODE) J 


roll mode request (MODER) 

control stick steering override (ACSS) 

measured roll stock (xa) 

roll trim indicator (IATRIM) 

roll attitude selection (^g) 

heading selection (^g) 

roll integrator control (IAINOK) 

roll SAS gain indicator (IAGAIN) 


PILOT 


gear 


extension indicator (IGEARX). 


none 


> 

> 


MODING 


CROSSFEEDS 



CONTROL^ 


OUTPUTS " 

AILERON 

STATUS 

LOOPS 



CROSSFEEDS 



^ l l 



aileron command (Sa^) 


{ aileron loop mode (MODA) 
aileron guidance mode (MODAG) 


to elevon commands 
^ bank angle loop command (^-,) 
and 

I bank angle loop command indicator (IPHC) 
and . 

Sa^, to rudder loop for turn coordination 


measured roll rate (p) 
measured bank angle (0) 
measured heading (ip) 


SENSORS 


Fig. 4- 13c Aileron Input/Output Interface Variables 




GUIDANCE 


guidance mode requests (MODRG) 1 
sideslip angle guidance cmnd Oq) } 
rudder guidance cmnd (6r^) 


rudder mode change indicator 
(MODRCH) 

BCS mode indicator (MODB) 


rudder channel mode requests 
(MODRR) 

differential pedal (xr) 
pedal trim indicator (IRTRIM) 
yaw SAS gain switch indicator 
(IRGAIN) 


bank angle loop command indicator ^ MOPING 
(IPHC) 


aileron command (Oa^,) 

bank angle command (<£(-.) 

bank angle for turn coordination (<£rp) 


measured yaw rate (r) 

A/ 

measured sideslip angle (/3) 
measured lateral acceleration (ay) 



CONTROL 

OUTPUTS 


■> rudder command (Sr^) 


! rudder loop mode (MODR) 

rudder guidance mode 
(MODRG) 


CROSSFEEDS 


m- none 


CROSSFEEDS 


SENSORS 


Fig. 4- 13d Rudder Input/Output Interface Variables 
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none 


GUIDANCE 


mode change indicators (MODXCH, 
where x = U, E, A, R, B ) 


manual controller commands 
(xs, xf, xe, xa, xr, St) 
manual trim indicators (IETRIM, 
IATRIM, IRTRIM) 


EXECUTIVE 


PILOT 


none 


none (across-the-board 
self-contained BCS code) 


> 


MODING 


BCS 

(BACKUP CONTROL SYSTEM) 
LOOPS 


T 


CONTROL^ 
OUTPUTS * 


STATUS 


CROSSFEEDS 


> 


CROSSFEEDS 


6e c , 6a c> 6r c 

6s C' 6t c. 

fcfj-, to control output 
subroutine 


mode indicators (MODXCH, 
where x = U, E, A, R, B) 


v £>a^ to elevons 


measured pitch rate (q) 
measured roll rate Cp) 
measured yaw rate Cr) 


SENSORS 


Fig. 4-13e Backup Control System Input/Output Interface Variables 




4.6 Functional Block Diagrams for Cruise Phase DFCS 

This section presents a detailed control-engineer's description of the cruise 
phase DFCS, using functional block diagrams. These include all the feedbacks, 
crossfeeds, limiters, switches, and compensation filters currently available with 
preflight parameter specification options. 

A number of simplifying conventions have been employed, particularly in the 
areas of filter and limiter descriptions. These conventions are summarized in 
Fig. 4-14. 

Tables 4-2 through 4-4 indicate the realizable cruise phase DFCS modes of 
operation, and the corresponding settings of all flags. Note for example, that while 
MODS = 3 (speed brake using u-filter for IAS-hold) and MODT = 3 (throttle using 
u-filter for IAS-hold) are both valid modes of operation, they must not occur simul- 
taneously. The mode change Initialization Subroutine (Fig. 4-25 through 4-28) 
handles all the setting of all derived flags, according to the logic indicated in these 
tables . 

Pilot input signal processing for speed brake throttle and flap controllers, for 
pitch and roll sidearm stick, and for pedals, is shown in Fig. 4-15. Pitch, roll and 
yaw manual trim circuits are shown in Fig. 4-16. Outer loop (semi-automatic or 
automatic) error signal generation is shown in Fig. 4-17. Note that the pilot "Select" 
option is implied (except for altitude) simply by "writing" into the appropriate storage 
location, e.g., H REF . 

Speed control, including throttle and speed brake options, is shown in Fig. 
4-18a, which indicates an input from TAEM-guidance interface signal processing as 
shown in Fig. 4-18b. Elevator control loops are shown in Fig. 4-19. Also shown 
in Fig. 4-19 are the generation of <t> r j- F for turn coordination (pitch and yaw rate 
command circuits), the generation of body flap commands (manual override of a 
Mach profile), and an elevator trim-tracker network akin to conventional stabilizer 
trim-tracking. (This trim indication is being studied for trim- initialization on 
emergency downmodes.) 

Figure 4-20 presents the aileron control loops, and Fig. 4-21 presents the 
rudder control loops, including provision for both /3 feedback and "0" feedback (a 
combination of yaw rate command, yaw rate feedback without washout filtering, 
lateral acceleration, and aileron). Note that pedal and trim are available even in 
auto mode (through the gain GRXAUT) for manual decrab maneuvers. 

Figure 4-22 indicates the override capability of the BCS mode, indicated by 
MODB = 1. Details of BCS loops (essentially the corresponding MANUAL or SAS 
loops with independently specified gains, coefficients, and nonlinearities) are given 
in Fig. 4-23 . 
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BLOCK DIAGRAM SYMBOL 


FUNCTIONAL EQUIVALENCE 


GENERALIZED FILTER 


U E 

FILTER 

U EF 

U E 

GUE(1) + GUE(2)Z _1 + GUE(4)Z -2 + . . 

U EF 


GUE 


1 +GUE(3)Z -1 +GUE(5)Z“ 2 +. . . 



TRAPEZOIDAL INTEGRATION 



RECTANGULAR INTEGRATION 



WASHOUT (1 st order) 



Fig. 4-14 Simplifying Conventions for Block Diagrams 
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BLOCK DIAGRAM SYMBOL 


FUNCTIONAL EQUIVALENCE 


LIMITERS 

SYMMETRIC: 



DACMN, MX 



DECMN, MX 



WORKING COEFFICIENTS/FILTERS 



Fig. 4-14 Simplifying Conventions for Block Diagrams (Cont) 
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Table 4-2 Logic for Speed Control Flag Settings 



Input Flags 

Derived Flags 


MODU 

MODUG 

MODS 

MODT 

MODUFT 

MODUFP 

BCS, UBCS 

1 

X 

1 

0 

X 

X 


1 

X 

1 

1 

X 

X 

UMAN 

2 

X 

2 

0 

0 

1 

(IAS - MAN) 

2 

X 

2 

2 

0 

1 

UASB 

3 

X 

3 

0 

1 

2 

(IAS - SB) 

3 

X 

3 

2 

1 

2 


4 


1 

5 

0 

0 

1 




1 

5 

5 

0 

1 




2 

4 

0 

1 

3 




2 

4 

2 

1 

3 




3 

4 

5 

1 

3 




4. 

2 

4 

2 

3 

AUTO, UAUT 



5 

5 

4 

2 

3 




6 

6 

0 

0 

4 




6 

6 

2 

0 

4- 




7 

6 

5 

0 

4 




8 

2 

6 

0 

4 




9 

5 

. . 6 

0 

4 




10 

6 

6 

0 

4 

UATH 

5 

X 

2 

3 

2 

2 

(IAS - TH) 




































































Table 4-3 Logic for Elevator Control Flag Settings 



Input Flags 


Derived Flags 


MODE 


MODEG 




CAS, ECAS 

3 


4 


4 

AUTO, EAUT 

4 


4 

EATT 

5 

EALT 

6 




Table 4-4 Logic for Aileron Control Flag Settings 


Input Flags 


Derived Flags 


MODA MODAG IPC IPHC IPSC 


ABCS 


ASAS 





























































































xf 


xf 


xf. 


PROCESSED PILOT 
CONTROLLER 
DEFLECTIONS 


Fig. 4-15 Block Diagram, Pilot Input Signal Processing 





























PILOT CONTROLLER 
TRIM BIASES 


Fig. 4-16 Block Diagram, Pilot Manual Trim Circuits 
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Fig. 4-17 Block Diagram, Outer Loop Error Signal Generation 


























MODUG 


°G 

(TAEM' 



MODUG 


Fig. 4-18b TAEM Guidance Processing for Speed Control 


4-38 


















































4-41 


IPHC 




Fig. 4-20 Aileron Control Loops 







































1 



MODB 



MODS 




MODB 


Fig. 4-22 


NOTE: 


4-43 




Fig. 4-23 Backup Control System 
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Backup Control System (Cont) 
















4. 7 Functional Flow Charts - Cruise Phase Initialization 

The figures in this subsection represent the detailed logic flow of part of the 
cruise phase DFCS at the engineering/functional level rather than at the coded soft- 
ware level presented in Appendix 4A. Implementation (especially with respect to 
programming language syntax) should not be restricted by overspecification. There- 
fore, in future revisions of this document functional flow charts similar to those in 
this subsection will replace the detailed flow charts in the appendix. 

Figure 4-24 presents the cruise phase segment of the routine which is called 
by the DFCS executive when the logic for a mission phase is to be completely initial- 
ized. Figures 4-25 through 4-28 present the cruise phase segment of the routine 
which is called for any phase or mode initialization; they interface with the structure 
shown in Fig. 4 A- 2. 
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Prepare working registers for integrator 
input gains (premultiplied by TF for 
rectangular or TF2 for trapezoidal), and 
for rate limits (premultiplied by TF); 


fT 4s ma , X T*6t max , T„* xt max (where ; t max = 6t max /GTX), 

SPEED 

T f *GSDDG, 

Tp-GTDDG 

TRIMS | 

Tp'-xe^, Tp 

*xa T , Tp*xr T 

Controller f T max 

Rate Limits { F xe 

Tp*xa maX , Tp* xr max 

BCS | Tp*6t mxB , 

T B *xe B , T B *xa B , T B *xr B 

Misc T^e^, 

(Tp/2) * GAIP 


(load '1 1 for automatic integrator 
IEINOK = IAINOK = 1 activation at gear-retract; 

load 'O' for manual activation) 



Fig. 4-24 Functional Flow Diagram, Load Manipulation Routine 
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Fig. 4-25 Functional Flow Diagram, Speed Initialization Routine 
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■ Type of Elevator 
Control Mode 
Change 


MODECH = 4 [ MODECH = 3 


MODECH = 2 


Load CAS damper gain 
GEQM = GEQMC 
( SAS and a-GUID will 
over -write) 


Initialize elevator 
trim for new mode 


2 (SAS) 


Which mode is requested • 

. 3 f (CAS) 


GEQM = GEQMS (IEGAIN ) 
JEGAIN = IEGAIN 


3 | (a-GUID) j" (c-GUID) 

GEQM = GEQM.A 


I b 4| (AUTO) 

CMODEGlb- 


-Which guidance mode I — 
is requested \& 


5 1 (Q-ATT) 6|(ALT) 

REF = 0 H REF = H 


2 (0-GUID) 


4 (H-GUID) 


Load 

coefficients 
and gains 
for CAS 


COMN. MX=CCCMN. MX CCMN, MX=CCAMN, MX 
GCE(i) = GCEC(i) GCE(i) = GCEA(i) 


GEC = GECC GEC = GECA > co <- I1Iclenl 

T2GEIC = TF2+GEC T2GEIC = TF2*GEC I ga ‘ nE \ 

* GEINTC * GEINTA ) g 


Load 

coefficients 


9 EN (1)= W 2) ‘°- »EN (1 >-H en <2>‘0. 


an N (l) = an^(2) = 0. 
r F = 4 W1 4 0 





C EN* 1) = c 

EN (2) 4 0. 


q N (U 4 q N (2) = 0. 

xe^, = xe - xe_ 


Initialize pitch 
stick filter for 
SAS, CAS; ignore 
deadband and limits 


6t CP = ^C' ^CF = 6f C 


<f> (l)-<b (2)=* =o. 

TN TN TF 


Initialize trim- 
tracker 

e T = e TF = DETRIM 


Initialize flags: 
ICC1, ICC2, 

ICF. ITHC, IHC 
(see Table 4-3) 

—r~ 

Set elevator integrator 
flag IEINT ; clear 
control stick steering 
flag: . IECSS = 0 


MODEG = MODEGR 
MODE = MODER 


Fig. 4-26 Functional Flow Diagram, Elevator Initialization Routine 
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Type of rudder 
control mode 
change 


| (modrch - 4^ ^MODRCH = 3^ 


^MODRCH = 2^ 


— ' 

Initialize rudder | 
| trim for new 
mode 


|r EN <1) = r EN (2) ~ °* 
r Wl = °- 

^N (1> = ^N (2) = °‘ 
^EN (1) = ^EN (2 ^ = °‘ 


Initialize 

filters 


W1 


EF 


6a CN (1) ='GDAC(l)*6a c 
! 6a CN < 2 ) 3 GDAC(4)*5a 



Initialize aileron 
I crossfeed 
washout filter 


Initialize pedal filter; 
ignore deadband and 
limits. 


Is it auto mode? 


GRX = GRXSAS 
IRW = IRWP 
GRRW = GRRWS(IRGAIN) 
JRGAIN = IRGAIN 
MODRR = 2 (SAS) 


Reinitialize yaw 
rate washout 
filter 


Load SAS gains 


MODR = MODRR 



Load AUTO gain 

Is pad load requesting 
’yaw rate feedback 
without washout? 


YES 


Guidance mode 
determined use IRW = MODRRG - 1 
of washout 


IRW = 0 


I No washout 



GRRW = GRRWA 
MODRG = MODRGR 


Load AUTO gain 


Fig. 4-28 Functional Flow Diagram, Rudder initialisation Routine 











APPENDIX 4A 


DETAILED FLOWCHARTS FOR CRUISE 
PHASE DFCS 


Detailed specifications for the cruise phase DFCS, in terms of FORTRAN 
flow charts of an operating version of the DFCS, are presented in Figs. 4A-1 through 
4A-20. Minute functional details are covered, but unfortunately they are often over- 
whelmed by implementational aspects (e. g. , higher-order- language syntax), and are 
thus less-than-ideal for specification purposes. 

Figure 4A-1 depicts the cruise phase preflight (or first-pass) load manipulation 
operations. Figure 4A-2 presents the skeletal structure .into which the remaining 
figures fit. In particular, there are five basic subroutines. These subroutines 
interface directly with the DFCS executive as shown in Section 1. 2. Their functions 
for cruise phase are as follows: 

Initialization Subroutine (Fig. 4A-3 through Fig. 4A-6) 

handles all first-pass (preflight turn-on or inflight phase-change) and 
mode change initialization activity. 

Filter Part 1 Subroutine (Fig. 4A-7 through Fig. 4A-10) 

handles all preparatory calculations for the Control Output Subroutine, 
including pilot input signal processing, filter output node computations, 
control stick steering override detection, etc. 

Control Law Subroutine (Fig. 4A-11 through Fig. 4A-13) 

handles all control output calculations in preparation for the Control Out- 
put Subroutine (Subsection 2. 2). 

Filter Part 2 Subroutine (Fig. 4A-14 through Fig. 4A-18) 

handles all filter wrapup processing and miscellaneous activity, such as 
body flap command calculations, bank-angle filter calculations, and re- 
engaging semiautomatic operation following a control stick steering override. 

BCS Subroutine (Fig. 4A-19 through Fig. 4A-20) 

handles sequentially, the Filter Part 1, Control Law, Control Output, and 
Filter Part 2 calculations for a 4-axis BCS mode of operation. 
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(*) For automatic integrator -activation at gear-retract; 
for manual reactivation load IEINOK = LAINOK = 0. 


Fig. 4A-1 Load Manipulator Routine Cruise Section 
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I 

cn 



Fig. 4 A -2 Overview 
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Fig. 4A-3 Initialization Subroutine, Cruise Phase Section, Speed Subsection 
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2 

SAS 

GEQM = 
GEQMS 
JEGAIN 

l(IEGAIN) 
» IEGAIN 



EAG 


< 

ECG 

" MODEGR ) 

2 



■ 

4[eHG 



5 

EATT 

ICCI = 
THREF 
ITHC = 

2 

’ =* THM 
1 


HH 

EALT 

| ICCI 

= 2 1 

HREF = HM 1 

1 IHC * 

- 


GCE(i) * GCEC(i) 
for i » 1 to 5 




■ 

ICCI =5 3 
IHC * 2 


ICCI = 3 
ITHC = 2 




ICC2 = 1 



THEN(I) * THEN(2) = 0. 


(Cont.) 



QMN(l) * QMN(2) = 0 
XEF = XEM - XEB 

I " 


( ignore XEDB) 


DSCF * DSC 
DTCF = DTC 


GEC = GECC 


GEC = GECA 

^ ^ _ 

T2GEIC = TF2*GEC*GEICC 


T2GEIC = TF2*GEC*GEICA 

| ET = ETF = DETRIM 


ANN(l) = ANN(2) = 0 


RMF = QW1 = 0 


IEINT = 1 

-"—I — 


IF(T2GEIC.EQ.O) IEINT = 0 


IF ((MODER. EQ. 3). AND. (UGEARX.EQ. 1).0R( IEINOK. EQ. 0))) IEINT = 0 


PHTN(l) = PHTN(2) = PHTF = 0 


CEN(l) « CEN{2) = 0 


I 


A } (Cont. ) 


IE CSS = 0 

T 


MODEG-MODEGR 
MODE = MODER 


Fig. 4A-4a Initialization Subroutine, Cruise Phase Section, Elevator Subsection 
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START 


MODECH 
\ =4 ^ 



MODE 


/ 5 


6 

EBCS 2 

ESAS 


DETRIM = GEXB * XETB 
+ DECUB 


4 EAUT 


IF (IEINT. EQ. O) EINT = EINT + GEC-'CC 


DETRIM - GEX^ET IF (IEINT. EQ. 0) EINT=EINT+GEC*GCX *XET DETRIM = DEC 
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XET = DETRIM / GEX 



3 ECAS 

5 EATT 

6 EALT 


XET = 0. 
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4 EAUT 


DEG = 0. 


EINT1= DETRIM 


Fig. 4A-4b Initialization Subroutine, Cruise Phase Section, Elevator Trim Subsection 












Fig. 4A-5a Initialization Subroutine, Cruise Phase Section, Aileron Subsection 
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Fig. 4A-5b Initialization Subroutine, Cruise Phase Section, Aileron Trim Subsection 

























Fig.4A-6a Initialization Subroutine, Cruise Phase Section, Rudder Subsection 
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Fig. 4A-6b Initialization Subroutine, Cruise Phase Section, Rudder Trim Subsection 









Fig. 4A-7a Filter Part I Subroutine, 


Cruise Phase Section Speed Subsection 
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Fig. 4A-7b Filter Part I Subroutine, Cruise Phase Section Speed Subsection, 
(TAEM Guidance Routine) 
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IF (ITHC. EQ. 2) THREF = THG 



THE = THREF - THM 




THEF=THEN (1) + GTHE (1) * THE 
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CC = GCX * (XEP + XETI 
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IECSS = ICC1 = 1 
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HEF = CC/GCH 


THEF = CC/GCTH 

HE = 0. 


THE = 0. 

HEN (2) =' -GHE (5) * HEF 


THEN (2) =-GTHE (5) * THEF 


Fig. 4A-8b Filter Part I Subroutine, Cruise Phase, 
Elevator Subsection, Outer Loops 
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Fig. 4A-9a Filter Part I Subroutine, Cruise Phase Section, Aileron Subsection 
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Fig. 4A-11 


Control Law Subroutine, Cruise Phdse Section, Speed Subsection 
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Fig. 4A-12 Control Law Subroutine, Cruise Phase Section, Elevator Subsection 
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Fig.4A-13 Control Law Subroutine, Cruise Phase Section, 
Aileron and Rudder Subsections 
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Fig. 4A-14 Filter Part II Subroutine, Cruise Phase Section, Speed Subsection 
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Fig. 4A-15 Filter Part II Subroutine, Cruise Phase 
Section, Elevator Subsection 
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(*) additional elements as required (SAS) 

(**) IRCHA1 (indicate request: change aileron integrator status! 


= 1 when AINT-OUT button is depressed; otherwise IRCHAI = 0 


Fig. 4A-16 Filter Part II Subroutine, Cruise Phase Section, Aileron Subsection 
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Fig. 4A-17 Filter Part II Subroutine, Cruise Phase Section, Rudder Subsection 
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Fig. 4A- 19a Cruise Phase BCS Subroutine, Speed, Elevator, and Flap Subsection 


4-78 




















Fig.4A-19b Cruise Phase BCS Subroutine, Aileron and Rudder Subsection 
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XTPPB = XTPB 



Fig. 4A-20 Cruise Phase BCS Subroutine, Wrapup Subsection 
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SECTION 5 


ENTRY PHASE 

b y 

Richard D. Goss and Greg L. Zacharias 


5 . 1 Introduction 

The DFCS design presented in this section is intended to provide stability-axis 
attitude control for the vehicle during entry, where this mission phase begins outside 
the sensible atmosphere (e.g., at approximately 400,000 ft altitude) and ends in the 
high subsonic flight regime. The control logic makes use of both the ACPS and the 
aerodynamic control surface torques, blended together in such a manner as to mini- 
mize ACPS jet activity, while providing capability for simultaneous parallel use 
when necessary. 

The basic design features include gain scheduling and control submodo switch- 
ing (both as functions of selected trajectory parameters) to account for variations in 
the vehicle dynamics and to allow for reasonably constant closed-loop transient 
response characteristics. In addition, except for very low dynamic pressure opera- 
tion, maximum use is made of the aerodynamic control surfaces to minimize ACPS 
fuel. Finally, the controller synthesis, and thus the design presented here makes 
advantageous use of the vehicle's aerodynamic torques resulting from trim attitude 
deviations, allowing the vehicle's rigid body aerodynamics to aid in attitude maneu- 
vers . 

This section is organized into three subsections. Subsection 5.2 provides an 
overview of the entry DFCS organization and operational characteristics, while Sub- 
section 5.3 is a detailed description of the entry DFCS, both functionally and as 
implemented in the present FORTRAN simulator. Finally, Subsection 5.4 is a brief 
summary of simulation experience to date. 

5 . 2 Overview 

Figure 5-1 provides an overview of the entry phase control logic structure. 

The interface, state estimation, and air data routines are described in Section 2 
(Subsection 2. 1), Section 3, and Section 2 (Subsection 2. 2), respectively. The pur- 
pose of the first two is to provide estimates of the stability axis angles and rates, 
along with the guidance commands, while the purpose of the latter is to provide 
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Fig. 5-1 Overview of Entry Control 




estimates of the trajectory parameters (trim angle-of-attack, Mach number, and 
dynamic pressure). 

Longitudinal control of angle-of-attack and pitch rate is provided by elevator 
(tandem elevon) deflections and pitch ACPS jet firings. The elevator logic uses con- 
ventional rate and position feedback (with feedforward trim integration) through 
scheduled gains to maintain relatively constant closed-loop pitch dynamics, while the 
ACPS logic uses phase-plane control to maintain a desired limit-cycle operation. 

Lateral control of sideslip and bank angle is provided by aileron (differential 
elevon) and rudder deflections, and roll and yaw ACPS jet firings. The aileron/ 
rudder logic uses rate and position feedback of the lateral variables (with feed- 
forward trim integration for center-of-gravity offset compensation) to provide atti- 
tude and rate control of both sideslip and bank. The ACPS logic uses phase-plane 
control to maintain a desired limit- cycle operation when there is low dynamic 
pressure, and sideslip rate damping to maintain lateral stability when there is high 
dynamic pressure. 

In both the lateral and longitudinal channels, blending of control by the surfaces 
and the ACPS is provided through a logic which determines ACPS inhibition on the 
basis of the magnitude of the current surface command, and on the basis of the 
vehicle’s Mach number and angle-of-attack. 

5.3 Description of Control Laws 

This section presents the control system block diagrams, equations, system 
parameters, and detailed flow charts defining the entry DFCS design. For conven- 
ient reference, this section is divided into four subsections. Subsection 5.3.1 and 
5.3.2 describe the longitudinal and lateral control laws, while Subsection 5.3.3 
specifies the fixed and scheduled parameters used in the present design. 

5.3.1 Longitudinal Control 

Longitudinal control uses both elevator (tandem elevon) and pitch ACPS jets 
to maintain the commanded angle-of-attack. As shown in Fig. 5-2, which illustrates 
a typical sequence of control activity during entry, there are two control submodes for 
the elevator and one for the pitch jets. 

Early in the entry, when the dynamic pressure is low, the elevator is used 
strictly for trim, so as to minimize pitch disturbance torques and thus fuel usage. 

The trim control law operates in a closed-loop fashion by slowly driving the elevator 
in the direction which nulls the pitch disturbance . The estimate of the disturbance 
torque is obtained by observing the pitch ACPS firing history. During this flight 
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Fig. 5-2 Controller Activity during Entry (Longitudinal) 









regime, the pitch ACPS jets, controlled by an angle-of-attack /pitch-rate phase- 
plane logic, provide attitude and rate control in response to guidance commands. 
There is no jet inhibition in this regime. 

When the dynamic pressure gets high enough, (q^), the elevator is used for both 
trim and transient control. Trim is initially maintained through the use of the trim 
logic of the previous submode. At a slightly higher dynamic pressure, (q^), trim is 
maintained through the use of an integrator acting upon the angle-of-attack error; 
the initialization value of the integrator is obtained from the value generated by the 
previous trim logic. Elevator control of transient errors is maintained through 
feedback of angle-of-attack and pitch rate, with appropriate gains. During this 
flight regime, pitch jet control identical to that of the previous regime may be uti- 
lized together with the elevator, or control with the jets may be inhibited. This is 
done by the blending logic which is essentially a two-sided deadbanded relay with 
hysteresis. The commanded elevator deflection is compared with preset fractions 
of its maximum and minimum values; should the command be "small", it is assumed 
that the elevator has sufficient control authority, and thus the pitch jets are not 
required. Conversely, should the elevator command be "large”, the pitch jet control 
is not inhibited, but is allowed to assist the elevator in controlling the vehicle. This 
blended control continues throughout the rest of the entry, so that the pitch jets are 
always available for control assistance. 

The logic provides for a gradual turn-on of the elevator control so as to avoid 
a switching transient when the second control submode is entered. 

The remainder of this section describes the longitudinal control laws, and their 
digital implementation, in detail. The block diagrams and equations defining the 
control logic for the elevator and pitch jets are given in Subsections 5. 3. 1. 1 and 
5. 3. 1.2. Subsection 5. 3. 1.3 consists of the controller coding summarized in flow 
chart form. 


5.3. 1.1 Elevator Control 

Figure 5-3a is a block diagram of the elevator control loop, showing the two 
control submodes, with switching determined by q^. 

When the dynamic pressure is less than q^, only an elevator trim setting is 
sent to the elevator actuator. This trim setting is determined by integrating the pitch 
ACPS commands (u^.), or » equivalently, by slowly driving the elevator so as to null 
the longitudinal disturbance torques. 

When the dynamic pressure reaches q^, closed-loop control of the angle-of- 


t 

See Subsection 


5.3.3 


for numerical values of the parameters 
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attack is provided by conventional position (a) and rate (a) feedback through a net- 
work which attempts to maintain a constant transient response throughout the flight 
envelope. Specifically, the desired closed-loop transfer function is chosen to be a 
well-damped second-order system, or 



desired 


a 

s 2 + 2 C d “ d s + tJ d 2 


(5-1) 


where (f,^, u^) are specified for acceptable performance (k & is not explicitly speci- 
fied). The gains shown in the diagram may be separated into one of three types: 

(1) prespecified by the desired closed-loop characteristics; (2) dynamic pressure 
dependent (l/q); or (3) scheduled gains which are functions of the aerodynamic flight 
regime (e.g., £^)t Trim control is accomplished in this submode either by the 
trim logic of the previous submode (if q <qj) or by a conventional clamped trim 
integrator acting on the angle -of- attack error, thus providing a steady- state trim 
elevator command in the absence of transient angle-of-attack errors. 

As can be seen from Fig. 5-3a and 5-3b, in order to allow a gradual blending 
in of active elevator control, the elevator command limits are functions of dynamic 
pressure which gradually open up from the trim setting to the physical deflection 
limits . 


5. 3. 1.2 Pitch ACPS Control 

Shown in Fig. 5-4 is the ACPS phase-plane switch logic representation used 
for control of the pitch jets. For longitudinal control, the phase-plane coordinates 
are angle-of-attack error and angle-of-attack rate, or 

(ef, e 2 ) = (a- ff^,, a) (5-2a) 

The phase-plane is separated into three regions (by the switch curves), in which the 
commanded jet firing is either positive (u = +1), negative (u = -1), or zero (u = 0). 
Thus, for the particular application of pitch jet control 

(u, U) = (u , U ) (5- 2b) 



respectively. 

The equations for the switch curves themselves are based on the fuel-time 
optimal solution to the double integrator control problem and are modified by 


* See Subsection 5.3.3 for numerical values of the parameters 


5-7 



TYPICAL TRAJECTORY 






deadband incorporation to ensure practical limit-cycle convergence (see typical 
trajectory of Fig. 5-4). The switch curves are defined by: 

y±: e 1 = =F (& 1 - 2u e 2 ) 

(5-3a) 

r± : e 1 = ± ( fi 2 + 2U e 2 ^ 

where, in this case, the deadbands are pitch attitude deadbands, or: 

**r 6 2> = < s «i- s „ 2 > ,5 ' 3b) 

and is a fuel-time weighting constant greater than unity given for pitch control by: 

<y - a (5- 3c) 

a 

The interface (or blending) logic between commanded pitch jet firings and the 
commanded elevator deflections is shown in Fig. 5-5. The approach taken is to 
inhibit firing when there is "sufficient" pitch control acceleration from the elevator. 
This measure of sufficiency is obtained by comparing commanded elevator (Se^-, 
shown in Fig. 5-3a) with threshold values to determine whether or not "excessive" 
elevator deflection is being called for. The hysteresis path included in the blending 
interface avoids chatter in elevator and ACPS activity due to cross-coupling effects. 
Note that this design includes an additional switch to preclude jet inhibition when the 
dynamic pressure is low and the elevator is used solely for trim. 

5. 3. 1.3 Digital Implementation 

Coding implementation for longitudinal control during the entry phase is 
summarized in flow-chart form in Fig. 5-6. 

5.3.2 Lateral Control 

Lateral control uses both aileron (differential elevons) and rudder, and roll 
and yaw ACPS jets to maintain small sideslip angles and follow commanded bank 
maneuvers. As shown in Fig. 5-7, which illustrates a typical sequence of control 
activity during entry, there are three control submodes for the aerodynamic control 
surfaces, one for the yaw jets, and three for the roll jets. 

Early in the entry, when the dynamic pressure is low, the aileron and rudder 
commands are zero. During this regime, the yaw jets, controlled by a bank angle/ 
yaw-rate phase-plane logic, provide attitude and rate control of the bank angle in 
response to guidance commands. The roll jets are controlled by one of two logics, 
depending on the yaw jet activity. When there are no yaw commands, the roll jets, 
controlled by a sideslip /sideslip-rate phase-plane logic, maintain small sideslip 
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Fig. 5-6a Logic Flow for the Longitudinal Control Channel 
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Fig. 5-6b Logic Flow for Elevator Control 
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Fig. 5-6b Logic Flow for Elevator Control (Cont) 
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Fig. 5-6c Logic Flow for Pitch ACPS Control 
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Fig. 5-7 Controller Activity during Entry (Lateral) 





attitudes and rates. When the yaw jets are commanded to maneuver the vehicle, 
the roll jets are controlled by a sideslip /sideslip- rate phase-plane logic, which is 
similar to the first submode except that compensation for yaw jet torques is in- 
cluded. Again, the purpose of this logic is to maintain small sideslip attitudes and 
rates. There is no jet inhibition in this flight regime. 

Later in the entry, when the dynamic pressure reaches a certain prespeci- 
fied value (q^) , the aileron and rudder commands are still zero. The yaw jets are 
controlled by the same logic as in the previous regime. The roll jets are controlled 
by one of two logics, again depending on yaw jet activity. When there are no yaw jet 
commands, the roll jets are commanded by a deadbanded relay logic which simply 
provides rate damping in sideslip. When the yaw jets are commanded to maneuver 
the vehicle, the roll jets are commanded by the same corresponding roll jet logic of 
the previous regime. There is no jet inhibition in this flight regime. 

Still later in the entry, when the dynamic pressure is higher (>qg), the aileron 
is used for both trim and transient control. Trim is initially maintained at zero. 

At a slightly higher dynamic pressure (q 2 ). trim is maintained through the use of an 
integrator acting upon the bank angle attitude error; the initialization value of the 
integrator is zero. Aileron control of transient errors (in a "reverse" sense, as it 
has come to be labeled ) is maintained through feedback of bank and sideslip attitudes 
and rates with appropriate gains. During this regime, the rudder command is null 
until the Mach number reaches a certain level (7??^), after which it is proportional 
to the commanded aileron deflection from trim. The proportionality variable is 
a constant. Also during this regime, roll and yaw jet control identical to that 
of the previous regime may be utilized together with the aileron and rudder, or 
control with the jets may be inhibited. This is done with the same type of blending 
logic used longitudinally; here, the commanded aileron deflection determines jet 
inhibition. Also as in the longitudinal case, there is a gradual turn-on of the allowed 
aerodynamic control authority. 

Still later in the entry, when the vehicle has crossed a prespecified curve in 
the Mach-alpha plane, only the aileron and rudder are used for transient control, 
the jets being inhibited for the remainder of the entry. Aileron and rudder control 
(in a "proverse" sense to distinguish it from that of the previous submode ) of 
transient errors is maintained through feedback of bank and sideslip attitudes and 
rates (to both surfaces), with appropriate gains. Lateral trim for this regime has 
not yet been implemented. 

The remainder of this section describes the lateral control laws, and their 
digital implementation, in detail. The block diagrams and equations defining the 
control logic for the aileron and rudder are given in Subsection 5. 3. 2. 1, while those 
for the yaw and roll jets are given in Subsection 5. 3. 2.2. Subsection 5. 3. 2. 3 con- 
sists of the controller coding summarized in flow-chart form. 
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5. 3. 2.1 Aileron /Rudder Control 


Figure 5-8 is an overview of the aileron/rudder control logic, illustrating the 
two control submodes, with switching determined by the flagl^^. This flag is 
specified by vehicle controllability considerations (see Subsection 5.4) and is defined 
by a curve in the Mach-alpha plane as shown in Fig. 5-9. 


5. 3. 2. 1.1 Early Aileron /Rudder Control 

Figure 5-10a illustrates the "reverse" control submode. When the dynamic 
pressure is less than q aileron trim is set to zero, as are the aileron and rudder 
commands. When the dynamic pressure reaches qg, aileron trim and rudder com- 
mand are still set to zero, but closed- loop control of bank and sideslip (by the 
aileron alone) is provided by position (|3, <j>) and rate (/3, 0) feedback of the lateral 
state variables, through a network which attempts to maintain a constant response 
throughout the trajectory. Specifically, the desired closed-loop transfer function 
for bank angle response is chosen to have four specified complex poles, or: 



desired 


(s' 


2 C 


l“l 


Ms 2 ' 

S + to^ 2 ) (s 2 + 2 i^gio^s + u 2 2 ) 


(5-4) 


where (^, are chosen for acceptable transient response (k^ is not explicitly 
specified and u , is a function of the airframe dynamics. ). This response is ob- 
tained, in the absence of rudder effectiveness, by making use of sideslip mistrims 
to bank the vehicle .. 


As with the longitudinal controller, the gains of Fig. 5- 10a may be separated 
into one of three types: (1) prespecified by the desired closed-loop characteristics 
(e.g., a^); (2) dynamic pressure dependent; or (3) scheduled gains which are functions 
of the aerodynamic flight regime (e.g., Ig)- For convenient reference, the gains 
of the first category are defined in terms of desired response characteristics with 
the following equation set: 


a 

a 

a 

a 


1 

2 

3 

4 


2 <M 2 + Ml )/( “l W 2> 

2( Ml + M2> 

“l 2+ “2 2+4! ^2 W l“2 


(5-5) 


where, as noted above, the (£., o^) are chosen for acceptable roll response. 
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Fig. 5-8 Overview of Aileron/ Rudder Control Logic 
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Fig. 5-9 Control Submode Switch Parameter Definition 
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NON-LINEARITIES: 


(1) Limiter; unity slope 
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Fig. 5-10a Aileron/Rudder Control Block Diagram (Reverse Mode, I 














Aileron trim is commanded zero until the dynamic pressure reaches qg , and it 
is then determined through the use of a conventional clamped integrator in the bank 
angle error channel. This allows the vehicle to trim to a non-zero sideslip angle 
(by commanding a non-zero trim aileron deflection) in the presence of body-axis 
roll disturbance torques due to lateral displacements in the vehicle's center-of- 
gravity. To avoid interpreting large transients (such as occur during a bank rever- 
sal) as lateral cg-offset induced mistrims, integration of large magnitude errors 
is inhibited by the \(f> switch. The gradual expansion of the aileron command limits 
during this submode as dynamic pressure increases is shown in Figs. 5- 10a and 
5- 10b. 


PARAMETER 



DYNAMIC PRESSURE, q 


Fig. 5-10b Aileron/ Deflection Limit Parameter 


When the dynamic pressure reaches q^ , later in this "reverse" submode, a 
rate command loop becomes effective when large attitude errors are detected ( as 
in a bank reversal). As shown in Fig. 5-10a, when the attitude error is greater 
than r^, the logic is no longer concerned with maintaining a bank attitude, but in- 
stead is configured to drive the vehicle to a desired bank rate. This allows the 
vehicle to complete the maneuver at higher rates than normally obtained with the 
attitude control loop providing the feedback. 

Finally, when the Mach number drops below 77^, the rudder is switched in to 
augment the aileron effectiveness. The aileron to rudder crossfeed gain is a 
constant . 
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5.3.2. 1.2 Late Aileron /Rudder Control 


Figure 5-11 illustrates the "proverse" control submode. Closed-loop control 
of bank and sideslip is provided by position (/3, <j>) and rate (/3, <j>) feedback of the 
lateral state variables, through a network which attempts to maintain a constant re- 
sponse throughout the remainder of the trajectory. Specifically, the desired closed- 
loop transfer matrix describing bank and sideslip attitude response to a guidance 
commanded step is chosen to be: 


(s 2 + a lt s + b lt ) 

( a 12 s + b l2 ) 


/3 


0 

(a 2i s + b 2i> 

(s 2 + a 22 s+ b 22 ) 


*-*G 


0 


where the a^. and b^ are chosen for acceptable transient response. It may be noted 
that in this ^'proverse" submode, the availability of the rudder allows for con- 
siderably more freedom in shaping the transient response, so that both sideslip and 
bank characteristics may be specified, in contradistinction to the previous submode 
(see Eq (5-4)). 


As shown in the block diagram, there are eight gains (r.^ and p^.) which are 
derived functions of both the desired closed-loop parameters (in Eq(5-6)) and of 
four selected scheduled gains. They are given by: 


It 

= a ll T lt + 

a 21 G l T 12 

Pll 

b ll T ll + b 21 G 1 7 12 

12 

= a 12 r il + 

a 22 G 1 T 12 

P 12 

b 12 7 ll + b 22 G 1 T 12 

21 

= a ll T 21 + 

a 21 G 1 T 22 

P 21 

b ll T 21 + b 21 G 1 T 22 

22 

a 12 T 21 + 

a 22 G 1 7 22 

P 22 

b 12 T 21 + b 22 G 1 T 22 


(5- 7a) 


where is a constant and the T-j are scheduled gains. Also shown in the block 
diagram are a constant gain, G^, two additional scheduled gains v ^ and v 2 , and one 


additional derived gain, a function of the 7^ : 

* " T 11 7 22 ' T 12 T 21 


(5- 7b) 


5 . 3 . 2 . 2 Yaw and Roll ACPS Control 


Yaw and roll jet activity is completely inhibited while in the "proverse" aileron/ 
rudder control submode = 1). This section describes the lateral ACPS 

control logic while in the "reverse" surface control submode. 


Yaw jet ACPS control synthesis is based on double integrator modeling of the 
bank angle dynamics, so that the switch logic presented for pitch ACPS control may 
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Fig. 5” 11 Aileron/ Rudder Control Block Diagram (Proverse Mode, I 
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be used. Roll jet logic consists of three models, one based on double integrator 
modeling of the sideslip dynamics, a second based on oscillator modeling, and a third 
which accounts for control axis cross -coupling. 


5. 3. 2. 2.1 Yaw Jet Control 

A double integrator model of the bank angle dynamics, combined with the 
assumption of turn coordination provided by the roll jets, allows for the use of the 
ACPS pitch jet control logic, with the appropriate redefinition of control parameters. 
Specifically, yaw jet control is defined by the phase-plane switch logic of Fig. 5-4 
and Eq (5-3a), where the coordinates are given by: 


(e l* e 2> = [-*lim sin V r ] 

(5-8a) 

and the switch curve parameters are given by: 

(u, U) = (u z , U z ) 

(5- 8b) 

(Sp 6 g ) = (6^, fi^ 2 ) 

(5-8c) 

a = a , 
0 

(5-8d) 


The blending interface logic is discussed below. 


5. 3. 2. 2. 2 Roll Jet Control 

Roll jet control of vehicle sideslip is accomplished by the use of three control 
modes. Two of these use a double integrator model for sideslip; the remaining one 
uses an oscillator model. 


The purpose of the first mode is to provide turn coordination when the yaw jets 
are firing (u z t 0). This is accomplished by the phase-plane logic of Fig. 5-12, 
which is a generalized (0, 0) plane. The figure illustrates a typical state trajectory 
and the resultant limit cycle (with switch curve overshoots). The equations for the 
switch curve u are: 


u: S0 


2~ [(8 m ) 2 - / jif S0 > 0 m 
[( 8 m ) 2 - 0 2 J if S0 < 0 m 


(5- 9a) 
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where 


U 1 = U z cos » T 

U- = U sin a - U cos a 

2x11 (5-9b) 

S = sgn(u z ) 

/3 m = sf 2U 2 6 1 

and where 6^ is defined below in Eq (5-10). 

The purpose of the second mode is to provide direct control over sideslip when 
the yaw jets are inactive (u z = 0) and when the dynamic pressure is low. With double 
integrator modeling, the phase-plane logic of Fig. 5-4 and Eq (5-3a) may be used, 
where the coordinates are given by: 

(ej, e 2 ) = (0, 0) (5 - 10a) 

and the switch curve parameters are given by: 

(u, U) = (u x , U x sin » T ) (5- 10b) 

(«j, 63 ) = (« 0 lt & q 2 ) (5- 10c) 

CT = CT^ (5-10d) 



Fig. 5-12 Roll Jet Control when u * 0 
& z 

It is convenient to label this control mode (used at low q and when u z = 0) as attitude 
hold, since both rate and attitude error are driven toward zero. 
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The third control mode takes advantage of the vehicle's dynamic stability in 
sideslip, which becomes significant once the dynamic pressure becomes sufficiently 
high. What is required is simply rate damping of the sideslip angular velocity as 
shown below, where /3 m is given by Eq (5- 9b). 


J8 • j 


- 

-U x 

/3 m 

u x 


-/§ m 

' U x “ 

- 




Fig. 5-13 Roll Jet Control when u z ~ 0 (q > q^) 

The control modes described above, defining lateral control with yaw and roll 
jets, are summarized in block diagram form in Fig. 5-14. Several points should be 
noted. First, the mode select block which defines the roll jet control mode to be 
used is defined by: 

If u i 0 then Mode = I 

z 

u z = 0, q < q 3 II (5-11) 

u z = 0 q > q 3 III 

where the modes are coordination, attitude hold, and rate damping, respectively. 

The three roll control blocks correspond with the above modes. A second point to 
note is the gating by the blending parameter X , which is defined by logic that is 
identical to that of the pitch channel. Finally, it may be noted that the logic dis- 
allows ACPS inhibition during aileron trim-only (null aileron command) operation. 

5. 3, 2. 3 Digital Implementation 

Coding implementation for lateral control during the entry phase is summarized 
in flow-chart form in Fig. 5-15. 
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AILERON /RUDDER CONTROL 
(Fig. 5- 15b) 


ROLL/YAW AC PS CONTROL 
(Fig. 5- 15c) 



Fig. 5-15'a Logic Flow for the Lateral Channel 
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Fig. 5-15b Logic Flow for Aileron/Rudder Control 
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BLENDING SUBROUTINE 
If 6a c < - or 6a c > -tj, X a * 1 
U-l 2 < 6a c <l 2 , X a - 0 
Otherwise, X unchanged 



Exit to 

yaw/roll AC PS logic 
(Fig. 5-l5c) 

Fig. 5-15b Logic Flow for Aileron/Rudder Control (Cont) 
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to Mode II logic to p^ode III logic 


Fig, 5-15c Logic Flow for Yaw/Roll ACPS Control 
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Fig. 5-15c Logic Flew for Yaw/Roll ACPS Control (Cont) 
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5.3.3 DFCS Parameter Specification 

All parameters used in control computations for the entry phase are either 
fixed or are scheduled on the basis of Mach number , TT[. Table 5-1 tabulates all 
entry- phase fixed control parameters. The scheduled control gains through £ g 
and and v are plotted in Figs. 5-16 through 5-23. The scheduled control param- 
eters through 722 are plotted in Figs. 5-24 through 5-27. As noted in Subsection 

5. 3. 2. 1.1, the lateral gains r. p.., and n are calculated as functions of the 7 . 

i] ij 1.1 

according to Eq (5-7). 

In an initialization pass through the entry control logic, several parameters 
which are functions of pad- loaded values are calculated and stored for future use. 

The Load Manipulation Routine for the entry phase is shown in Fig. 5-28. 

5.4 Performance Summary 

Shown in Fig. 5-29 are three curves in the Mach-alpha plane which define re- 
gions of lateral controllability using differing control modes. Also shown is a 
typical late entry Mach-alpha profile. The curves were derived by conducting simu- 
lated maneuvers with the entry DFCS and the NR 00089B (Feb. version) vehicle at 
various Mach numbers and angles-of-attack. The curves separate desirable from 
undesirable simulated vehicle response. 

Curve A indicates the lower boundary (in the Mach-alpha plane) to which the 
vehicle may be flown using ailerons alone for lateral control. By driving the rudder 
proportionally to the aileron command (with a crossfeed gain of 2 ), the region of 
"reverse" control may be considerably extended farther down the entry profile 
(curve B). Finally, the region of controllability using the "proverse" control mode 
is bounded (above) by curve C. The resultant overlapping region, defined by curves 
B and C, provides the motivation for the switch curve choice shown in Fig. 5-9. 

It should be noted that an on-going effort in determining design sensitivity to 
aerodynamic coefficient uncertainties will reduce the area of the overlap region. 
However, the control design used in deriving the curves of Fig. 5-29 was not "tuned" 
so as to maximize the resultant overlap area. It is anticipated that future efforts 
in this direction will enlarge this region for the nominal aerodynamic coefficient 
values, while counteracting the shrinking effect of off-nominal aerodynamics. 
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Table 5-1 Fixed Parameters 


Symbol 

Value 

Dimension 

U 

X 

1. 65 

deg / sec 3 

U 

y 

1. 32 

deg/sec 2 

U z 

0. 66 

deg/sec 2 

6 al 

0. 75 

deg 

CM 

a 

lO 

1 . 0 

deg 

6 pi 

0. 2 

deg 

6 P2 

0. 4 

deg 

6 tl 

0. 75 

deg 

* <t>2 

1. 5 

deg 

o 

5 

- 

a 

°p 

2 


% 

5 


h 

0. 8 

- 


0. 6 

_ 

2 

f 3 

0. 6 

- 

f . 

0. 8 

_ 

4 

*1 

10 

lbf/ft 2 

^2 

20 

lbf/ft 2 

^3 

5 

lbf/ft 2 


20 

lbf/ft 2 

qg 

30 

lbf/ft 2 


80 

lbf/ft 2 
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Table 5-1 Fixed Parameters (Cont) 


Symbol 

Value 

Dimension 

*1 

6 

-1 

V 

-1. 0 

sec 

-1 

K 0l2 

0. 2 

sec 

-1 

K . 
0 

0. 01 

sec 

-1 

“d 

0. 7 

sec 

-1 

«1 

0. 7 

sec 

-1 

W 2 

0.7 

sec 

^d 

0. 7 

- 

?! 

1. 4 

- 

^2 

0. 7 

- 

r l 

5 

deg 

r 8 

25 

deg 

V 

5 

deg 

r 02 

25 

deg 

r 0'3 

8 

deg / sec 

C RUD 

2 

-1 

a U 

1. 0 

sec 

-1 

a l 2 

0. 0 

sec 

-1 

a 21 

0. 0 

sec 

-1 

a 22 

1. 3 

sec 

-2 

b ll 

1. 0 

sec 

-2 

b 12 

0. 0 

sec 
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Table 5-1 Fixed Parameters (Cont) 


Symbol 

Value 

Dimension 

b 2 1 

1.0 

-2 

sec 

b 22 

0. 8 

-2 

sec 

G 1 

0. 133 

- 

G 2 

29. 0 

lbf-sec 

6e mLn , 6e max 

see 


c max . max 
6a , 6r 

Subsection 2-2 











Fig. 5-20 Gain Schedule for 



Fig. 5-21 Gain Schedule for £g 


MACH NO. 


MACH NO. 
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■ 0.0035 


■ 0.0045 
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Fig. 5-29 Controllability Regimes 







SECTION 6 


ON-ORBIT PHASE 


b y 

Steven R. Croopnick, Donald W. Keene and John J. Turkovich 


6.1 Introduction 

The DFCS (Digital Flight Control System) will control the attitude and the 
translation of the Space Shuttle during the various orbital activities by using either 
the ACPS (Attitude Control Propulsion System) or the OMS (Orbital Maneuvering 
System). For the purpose of control, there are two orbital phases: the ACPS phase 
when the jets are the control effectors, and the TVC (Thrust Vector Control) phase 
when the gimbaled engines are used for control. The ACPS phase is presented first, 
in Subsection 6.2, followed by a preliminary description of the on-orbit TVC phase in 
Subsection 6.3. Digital implementation of the ACPS phase, presented in Subsection 
6.4, and preliminary testing was accomplished using structured programming techniques 
in the HAL language. All symbols relevant to the on-orbit phase may be found in a 
separate glossary at the end of this section. 

6.2 On-orbit ACPS Phase: Engineering Description 

A block diagram of the DFCS during the on-orbit ACPS phase is 
illustrated in Fig. 6-1. The on-orbit attitude maneuver routine or another inter- 
face routine (see Part II) supplies to the DFCS: 

1) a set of reference gimbal angles, [0^] 

2) a set of incremental reference gimbal angles, [A 0^] 

3) the desired angular velocity vector, 

4) a set of bias gimbal angles, ] 

5) a bias acceleration vector, a 

-m 

Every tj agt seconds the DFCS adds [A0^] to [0^] to give a new set of reference gimbal 
angles. The DFCS then compares these to the measured set of gimbal angles, [0 q A L 
and computes the angular errors in control axes. The bias angles are added td these 
errors in control axes and fed to the phase plane switching logic. Also fed to the 
phase plane switching logic are the angular velocity error, the net estimated 
disturbance acceleration, 9' dni . the magnitudes of control accelerations, [&q], and 
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FROM 
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ERRORS 
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GIMBAL 

ANGLES 



ROTATIONS 


Fig. 6-1 On-Orbit ACPS 












the last desired change in angular velocity. The angular velocity error, « , is com- 
puted by subtracting the desired angular velocity, from the estimated angular 
velocity, «, and transforming the result into control coordinates. The net estimated 
disturbance acceleration, , is computed by adding the bias acceleration, a , to 
the estimated disturbance acceleration, and transforming the result into control 
coordinates. From these inputs the phase plane switching logic determines the 
desired angular rate changes, [Au] , in control coordinates. From these desired 
rate changes the jet selection logic determines the appropriate jets to be fired, if 
any, and the corresponding firing times, and then turns the jets on and schedules 
their turn-off. 


6.2.1 ACPS Control Axes 

This subsection presents a technique for decoupling the ACPS control of the Shuttle 
Orbiter into three essentially independent axes. This technique permits one to treat 
the rigid body rotational dynamics of the vehicle, normally represented by a 6^ order 
system, as three independent 2 nu order systems, thus justifying the use of phase 
plane control logic. A similar technique was used in the design of the Lunar Module 
autopilot. The jet selection logic which has been implemented to utilize the de- 
coupled axes is limited to specified jet pairs and provides only attitude control, not transla- 
tion, in its present form, A more general scheme will be used in the future. 


The ACPS thruster group locations for the NR ATP baseline configuration are 
given in Subsection 2.3. Define a +u thruster combination as one or more thrusters 
from group 2 and/or group 12: e.g., jet #14 and jet #31. Also define a -u thruster 
combination as one or more thrusters from groups 1 and 11: e.g., jets #10 and 25. 
The only restriction in selecting the jets is that the torque produced by the +u com- 
bination is opposite (or nearly opposite) in direction to that of the -u jets. Similarly, 
select a +v combination from groups 4 and 10, e.g., #9 and #26; and a -v combina- 
tion from groups 3 and 9, e.g., #13 and #32. Finally, select a +w combination from 
groups 6 and 8, e.g., #5 and #18; and a -w combination from groups 5 and 7, e.g., 
#6 and #17. The suggested thruster combinations give moderate rotational control 
authority while minimizing the net translation. 


The control axes are defined by the accelerations produced by the u, v and w 
jets. More specifically, let 


-BASIS 

-BASIS 

W 

-BASIS 


unit vector in the direction of the +u axis 
unit vector in the direction of the +v axis 
unit vector in the direction of the +w axis 
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* 

J 


orbiter inertia matrix 


— P2U’ — M2U’ — P2V’ control torques produced by the +u, +v, 

— M2V’ — P2W’ — M2W - W truster combinations 

then 


-BASIS 

= unit 

*-i 

(J 1 

( -P2U ’ 

Mm2U )/2) 

-BASIS 

= unit 

*- 1 
(J 1 

(Mp2v " 

-M2V )//2) 

-BASIS 

= unit 

(J 1 

(M P2 w " 

-M2W )/2) 


(The components are given in vehicle axes.) 

The control axes orientations are shown in Fig. 6-2 for the thruster 
combinations given above. Note that the +u accelerations are not exactly opposite. 
Neither are the +v accelerations. 

The transformation from control axes to vehicle (or pilot) axes is defined as: 
* 

-CP = ^-BASIS’ -BASIS’ -BASIS* 

The transformation from vehicle (or pilot) axes to control axes is simply: 

Mpc = * M CP* 

The attitude and rate errors for the autopilot are derived in vehicle axes. In 
order to utilize the decoupling technique these must be transformed from vehicle 
axes to control axes, i.e.. 



6.2.2 State Estimator for the ACPS On-Orbit Phase 

Control of the shuttle vehicle in orbit requires estimates of vehicle angular 
rate, u, and vehicle angular disturbance acceleration, Oj. These quantities are 


6-4 






determined from measurements of the IMU gimbal angles. The procedure used 
is based on Kalman filter theory in which angular measurements are mixed with 
extrapolated state estimates to form a filtered estimate of the state. The process of 
extrapolation and measurement incorporation is described first in terms of a single 
axis. The complete three axis filter as implemented in the DFCS is discussed at 
the end of this subsection. 


6.2.2. 1 Vehicle Model 


In estimation theory it is necessary to define a model of the process being 
observed. In the case of the shuttle vehicle, we will assume that it is a rigid body 
with the following equations of motion for one axis : 




0 

1 

0 




0 


0 


0 

CO 

= 

0 

0 

1 


u> 

+ 


+ 

a 

m 

+ 

n co 

a d 


0 

0 

0 


~l 

1 


0 


0 


% 


( 6 - 1 ) 


where 




co 



= vehicle altitude 

= vehicle angular velocity 
= angular disturbance acceleration 

= control acceleration (produced by ACPS jets) 


a - known bias acceleration 
m 

n co’ n o' = driving noise terms 


6 . 2 . 2 . 2 Extrapolation 


The discrete filtering process can be divided into two parts: extrapolation 
and measurement incorporation. In the extrapolation process the estimate of the 
vehicle state, j£(tg), given at some time, t Q , is integrated forward to the next meas- 
urement incorporation time, t^. This integration includes the effect of any known 
biases or control inputs during the interval tg, t^. The result of the extrapolation 
process is a predicted estimate x'(t^) . In our case x'(t^) is determined by integrating 
Eq. (6-1), i.e.. 


x'(t t ) =4>(t 1 , tg)x(tg) +/ <t> (tj, r) 

*0 


Q c (t) | 


dr 



0 

a (r) 
m 


dr 


0 


(6-2) 
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where 


iV 


0 B^V 

“V 


a d (t 0> 


= state estimate at t 


0 


(6-3) 


and 


®(t 1 , t) = 


(t 1 - t) 
1 
0 


(tj - T)‘ 


(t 1 - T) 


= transition matrix 


(6-4) 


or^(r) is the control acceleration applied by the ACPS jets. A typical plot of a^(r ) 
is shown in Fig. 6-3. It is assumed that the control acceleration can have one 
of three values + I a^l, 0, - lor I . It is also assumed that if a jet goes on, it turns 
at the beginning of the interval, tg. It is further assumed that when the jet goes off 
it remains off until the end of the interval. With this jet model we can now integrate 
the second term in Eq. (6-2)- 


t« 

1 * 

0 


l“ C l SIGNI “c> *JET (1T 

I <><4' T) 
*0 

a c (T) 

0 

dr = 

|a c | SiGN(a c ) t JET 
0 


(6-5) 


where 


A 

A T = t. 


fast 


a (t) is a known bias acceleration which is assumed to be constant over the 
m 

over the integration interval: 


ol(t) = a (t n ) 
m m u 


The third term in Eq. (6-2) then reduces to 


X * (t X - T) 

0 ' 

ot° (t) 
m 

dT - 

“m (t 0»‘ T <TT>' 

“mV 1T 


0 . 


o 


Before grouping terms, define 
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Aw dm (t 0 ) = C& d (t 0 > + “m (t 0 )] * T 
~ SfGNtecUo)) tj ET 
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°C (T) 



Fig. 6-3 ACPS Jet Control Acceleration Model 
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Combining results, we find 


■ ‘W + “ (t 0 )4T + “dm 11 ® 1 -f- + 4 “ cV <‘ T » ,6 - 7a) 

A A 

w / (t 1 ) = Oi (t Q ) + Aai dm (t Q ) + Aa^,(t 0 ) (6 - 7b) 

a d (t^ = oi d (t Q ) (6-7c) 


6. 2. 2. 3 Measurement Incorporation 

For the measurement incorporation process, the extrapolated state is combined 
with the measurement at time t in the following manner: 

W ‘ " 9 b V + K e <«<*„> - ‘ 8 bV> <6 - 8a) 


»«-> - »'<*■> + - 6 bV /«„-*/> 


n 


* p ' 

n 


“d (t n> = “d V + K « (? V - e B (t n )) 'W 


(6 -8b) 
(6 - 8c) 


where 

K q = position gain (constant) 

= velocity gain (constant) 

= acceleration gain (constant) 

9 (t ) = angle measurement at time t 
n n 

t^ = time of the last measurement incorporation 


This measurement incorporation process is similar to the linear Kalman filter 

if the measurement interval, t - t., is fixed. However, for the Shuttle filter, meas- 

n 

urements will be incorporated only if 

l^n* ‘ ®B (t n^ I > ®BT * 6 ' 9 * 

where 

®BT = " tra P" angle 

If this condition is not true, then no measurement is processed and 


W 

= 0 B (t n ) 

(6- 10a) 

“V 

= 

(6 -10b) 

A 

aAt ) 
a n 

= “d (t n> 

(6-10c) 


6-9 



This additional logic helps to reduce the effect of quantization of the gimbal angle 
measurements by "trapping" small- amplitude noise before it reaches the estimate. 

6. 2. 2. 4 The 3 -Axis Filter as Implemented in the DFCS 


The 3-axis state estimator is basically a simple extension of the single axis 
filter described above. Due to the kinematics of Euler rotations, however, the gim- 
bal angle measurements are not simply related to the vehicle body rates. We can 
circumvent this problem by defining a set of "body angle" measurements. The body 
angles are defined as the integral of the body rates: 


e Bi = I dt 
®B 2 = I w 2 dt 
e B 3 = I “ 3 dt 


(6-lla) 

(6-llb) 

(6-llc) 


To a good approximation we can calculate the body angles in recursive fashion as: 


(n) 

1 


9 (n-1) 

B 1 


1 

SIN(9 GA 3 > 

0 


®GA 1 (n) © ®GA 1 (n - 1) 

(n) 

2 

= 

0 R (n-l) 
B 2 

© 

0 

COS<9 GA )COS<8 GA3 > 

SIN( 8gAi ) 


®GA 2 (n) © ®GA 2 (n “ 1) 

(n) 

3 


0T. (n-l) 
B 3 


0 

- SINIS GA 1 )COS<e GA 3 ) 

CQSIKgAj) 


®GA 3 (n) © ®GA 3 (n " 1) 


( 6 - 12 ) 

where 


and 


®GA 1 

®GA„ 

“-3 


outer gimbal angle 
inner gimbal angle 
middle gimbal angle 


© indicates modular addition 
© indicates modular subtraction 


(Modular addition and subtraction of two angles produces a result which lies in the 
range: -180° < 9 £ 180°. ) Measurement incorporation is performed as before (Eq. (6-8)) 
except that the measurements are now the computed body angles (Eq.(6-12)). For the 
x axis. 


if 


«B V' 9 B VI > 9 

X X 


BT 


(6-13) 
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then 

«B (t n> ■ S B V 

X X 

© K «pB «„» © »B <‘n>] 

L X X J 

(6- 13a) 


ai x (t n ) = < (t n } + 

K J?B 

L X 

«„» © 6 b «„>]/«„- > 

X X 

(6- 13b) 


*dx (t n ) = 0 i dx (t n ) + 


A 

<*„> © *B <*„>]/<'» - l t > 2 

X J X 

(6-13c) 


*«r l " 



(6- 13d) 

otherwise. 


V - 1 

* 

= e B V 

X 

(6-14a) 



w x (t n } 

= < (t n> 

(6- 14b) 



“dxV- 

-■«dx (t n> 

(6- 14c) 


In the above equations 

A 

0g = estimated value of the x body angle 


0 B 


S B 




CO 


Of. 


'dx 


“dx 


extrapolated value of the x body angle 
"measured" value of the x body angle 

estimated x-axls body rate 

extrapolated x-axls body rate 

estimated x-axls component of disturbing acceleration 
extrapolated x-axls component of disturbing acceleration 
current time 


n 


t . = time of the last measurement Incorporation for the x-axls 

X 

The y and z axis computations are similar. Note that a separate "trap" Is required 
for each axis. 

After the appropriate jet control Is exercised the extrapolation of the vehicle 
state, which depends to a large extent on the ACPS jet selection and timing logic used 
in the Shuttle autopilot, is performed. It is currently assumed that the jet firings 
produce accelerations about the u, v and w axes defined in Subsection 6.2.1. That is, 
the angular acceleration produced about the vehic’e axes by the jets is the sum of the 
three components: 
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(6-15a) 


“cu " ^ *-P2U — M2U^ 2 




J (M p2V - Mm2V^ 2 


^CW ' ^ ( Mp2W " -M2W^ 2 


(6-15b) 
(6-1 5c) 


where J Ls the Inertia matrix and the M's are the torque vectors defined In Subsection 
6.2. 1. It Is also assumed that the restrictions on the timing of jet firings which were 
given for the single-axis filter apply. If we define 

tjET = ti-rne of firing of the u jets during the interval t^ < t < t^ 

tjET = time of firing of the v jets during the interval tp < t <t^ 

2 

tjET = time of firing of the w jets during the interval t Q < t < t^ 

3 


and 


SGNC, 


sign of the desired u control acceleration 


SGNCg = sign of the desired v control acceleration 
SGNC 3 * sign of the desired w control acceleration 


then by analogy with the results of the single axis filter: 


and 


^-u (t 0 } = 

l-cu' 

SGNC 1 

Vet 

Aw v (t 0 ) = 

l“cvi 

SGNC 2 

^JET, 

^wV = 

IScwl 

SGNC 3 

Vet. 


A 




^dmV = ( “d (t 0> + «m (t 0» AT 


(6- 16a) 
(6- 16b) 
(6-16c) 

(6-16d) 


Note that these are vector quantities. The bias acceleration, a m , is computed by 
the attitude maneuver routine as 


“m = -J'^xJ^) 


(6-17) 


where 0)^ is the desired maneuver rate. With these definitions, the complete three 
axis extrapolation equations can be written as 

AT 


— B '*l' ' i B (t 0> © 


A 

w(t 0 ) t 

t 


AT + Aw d m^V 
JET. 


u xj -L 

+ Aci^j (tQ) (AT ^ (AT - — J 


JET r 


+ Aw w (t 0 ) (AT 


Vet, 

-- S- 5 -*. 


(6- 18a) 
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(6-18b) 


o> ' (tj ) = a)(t Q ) + L w dm (t 0 ) + £« u (t Q ) + ^W v (t Q ) + Aw w (t 0 ) 

^ A 

= « d (t 0 ) (6- 18c) 

Note that these quantities are vectors with components in the x, y, z vehicle (body) 
axes. 

A 

The estimate of disturbance acceleration, c^j , may be accurate enough for the 
powered flight phases, but it may be too noisy for use during drifting flight. If an 
estimate is needed during coasting flight to counteract such disturbances as gravity 
gradient and aerodynamic torques, then rate gyro measurements may be required. 

If not, the estimation process may be disabled by using appropriate filter gains. 
Alternatively, the gain, K^, simply can be reduced during coasting flight. In any 
event, the bias acceleration estimate, a , should be used during coasting flight 
maneuvers to reduce fuel consumption. 

6.2.3 Control Laws for the On-Orbit ACPS Phase 

The object of the phase plane switching logic is to compute changes in angular 
velocity so that by turning on or off the appropriate jets to achieve this velocity change 
the vehicle limit cycles about a desired attitude and attitude rate. The two dimensional 
phase plane is applied to the decoupled control axes described in Subsection 6.2.1. Limit 
cycling is desirable since it would be extremely costly to attempt to maintain the de- 
sired attitude and rate exactly. By designing appropriate switching lines into the 
phase plane, the vehicle can be driven into the desired limit cycle. Two separate 
logics are used depending on the size of the position and rate errors. For large er- 
rors, a switching logic has been designed to drive the state into the small-error con- 
trol region. The small-error switching logic then drives the state into a limit cycle 
about the desired attitude and attitude rate. The distinction between large and small 
errors will be quantified later in this section. 

The small-error switching logic has the following features: 

1) The design is based upon the optimal switch curves for minimizing a 
weighted combination of fuel and response time. The optimal curves 
are separated by a deadband in order to allow for limit cycling. 

2) The appropriate deadband can be chosen in order to 

a) fix the limit cycle period 

b) fix the limit cycle maximum rate 

c) fix the maximum attitude excursion while limit cycling 

3) The average attitude error is minimized by biasing the placement of 
the' switch curves. 
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4) Parabolic switch curves approximate the trajectories which the state 
follows. However, hysteresis logic prevents performance degradation 
due to errors in these approximations. 

The design incorporating these features is described in the following subsections. 

In Subsection 6.2.3. 1 two idealized limit cycles are presented. One limit cycle is 
in the presence of positive disturbing accelerations while the other is in the presence 
of zero disturbing accelerations. For each limit cycle the average attitude error will 
be derived. In Subsection 6. 2.3. 2 the design of switching curves to drive the state into a 
limit cycle with zero average attitude error is described. This design is referred to 
as the small-error control logic. The large-error control logic is described in 
Subsection 6 . 2 . 3 . 3 . 

6. 2. 3.1 Idealized Limit Cycles 

An idealized limit cycle for a positive disturbing acceleration is illustrated in 
Fig. 6-4. The coordinates 9 g and are chosen to be the state, where 0 g represents 
the angular position error and m represents the angular rate error. is the 

attitude variation in the course of limit cycling. In the presence of a positive disturb- 
ance acceleration, , the state is driven along the trajectory from A to B. At point 
B a negative control acceleration, -U^,, is applied to return the state to point A. 


0 = 

e 


u)„ 



DB 


Fig. 6-4 Idealized Limit Cycle for a Positive 
Disturbing Acceleration 
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The average attitude error, "X, of the limit cycle in Fig. 6-4 is found by deter- 
mining the average attitude error for each of the trajectories OB and BO' (called X, 
and Xg) and then "weighting" them according to time to calculate the net average, X. 

Determination of 

X. is, by definition, 

X 1 = T, / J u D* 2dt ,6 - 19> 

1 0 


tp the time to travel from point O to point B, can be calculated from the value of the 

9 abscissa of point B, B a , since 
6 8 e 

B e =^U Dtl 2 (6-20) 

e 

B fi Is found by first rearranging 6 = McoJ and 0 = f 9 (u> ) such that, 

e j. e 0 c, e 

“e ■ Si(« e >. » e ■ g 2 <e e ) (6-21) 


Then, at point B, g 1 (6 e ) = g 2 (0 e >, and ' 



and from Eq(6-20) 


Using t^ in Eq(6-19) 


2X db (U c - U D ) 



3 X DB 



( 6 - 22 ) 


(6-23) 


(6-24) 


Determination of X 2 : 


The time average attitude error for the trajectory segment BO' can be found by 
setting t = 0, at point B., Then, 

h 

X 2 * r, ' C x o +i o l+ ^ (u D- u c )t2,d, 

1 0 


or. 


x r 


X + 


X t, 

iO 2 




(6-25) 


X Q is defined in Eq (6-22) and X Q can be determined from. 
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Wi = 


+/ 2X DB U D (U C “ U D^ 


(6-26) 


To find tg one can use the following relation. 


X o + < U D - U C )l 2 * 0 


from which it follows. 


/ 2X DB U D 


C ' C U D' 


(6-27) 


Substituting the results of Eq (6-22), (6-26) and (6-27) into Eq (6-25), yields: 


X 2 X DB\ 1 ' 3U. 


(6-28) 


Determination of X: 

Finally, the total average attitude error, X, can be found from: 


^ _ x i‘i + x 2 t 2 


(6-29) 


Substituting the results of Eqs (6-23), (6-24), (6-27), (6-28) into Eq (6-29), yields: 




(6-30) 


This is the average attitude error for the limit cycle in Fig. 6-4. Thus, the average 
attitude error would be zero if the intercept of the curve on the right, which is the 
trajectory along which the control is applied, were shifted from point O' to have an 


abscissa of 


DR ' X DB ” X 


(6-31) 


A similar result can be obtained for the case of a negative disturbance acceleration. 
Equations (6-30) and (6-31) are used in the DFCS in order to minimize the average 
attitude error. 


The attitude variation, X^g, could be chosen as follows: 


1) To achieve the cycle period, T: 


set X T 


D 1 t 2 


(6-32) 


2) To achieve the limit cycle maximum rate, X 


set X T 


| u r (x r 

1 C 1 max 

2|U d | (jU c | - |U D |) 


(6-33) 


3) By simply setting X—.-. to the allowable attitude variation. 
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Currently, the final option is employed in the DFCS. 

An idealized limit cycle in the presence of no disturbing accelerations is illus- 
trated in Fig. 6-5. From point A the state is allowed to drift 



Fig. 6-5 Limit Cycle in Presence of No 
Disturbing Accelerations 

to point B where a negative control acceleration, (in the form of a minimum 

impulse jet firing) is applied to drive the state to point C. From point C the state is 
then allowed to drift to point D. At point D a positive control acceleration, +U C , is 
applied to return the state to point A. By choosing DL = DR the state limit cycles 
with zero average attitude error. The generalized limit cycle for positive disturbing 
accelerations is made to reduce to the limit cycle for zero disturbing accelerations 
by selecting 

DL = X + i X DB (6-34) 


6. 2. 3. 2 Small-Error Control Law 

Control regions for relatively small errors are illustrated in Fig. 6-6, defini- 
tions of curves are summarized in Table 6-1, and corresponding control actions are 
given in Table 6-2. All logic in this subsection is formulated in terms of positive dis- 
turbing accelerations. In the actual implementation, (see Subsection 6.4), negative 
disturbances are treated by inverting the signs of the state variables 0 g and 
Biasing the right and left deadbands, DR and DL, to null the average attitude error 
in accordance with the equations derived in the preceeding subsection, defines the 
placement of the switching curves BOUNDARY_II_TO— COAST and BOUNDARY_III_TO_ 
COAST. These curves approximate the trajectory of the state during a firing. When 
the state lies in Region II, control is applied to drive the state to the curve TARGET— 
DISTURB. The vehicle is then allowed to coast. 
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Table 6-1 Definitions of Boundaries of Small-Error Control Regions 


Boundary/ Parameter 
BOUNDARY. I_TO_II 

BOUNDARY- II _TO_COAST 

BOUNDARY- 1 -TO—COAST 

BOUNDARY-III-TO-COAST 

BOUNDA R Y_ IV_TO_III 

BOUNDARY-IV—TO -COAST 

TARGET-DISTURB 

WHERE 


Definition 


u 


2<(u„! - |u n |) 


+ DR 1 


'C 1 1 D 1 

2 

UL 


8 = - 
e 


0 = + 
e 


e = + 

e 


9 = + 

e 


2(|U c i - |U D |) 

2CD{|U G j + |U D |) 
2 


+ DR 


+ DR' 


W, 


2<|U C I + l U D l 

2 


2(iU c | + |U D |) 


- DL 


DL' 


u> 


0 = 

e 


2CD(|U C ] - |U D |) 
u > 2 


- DL' 


0 = + — - 

6 2|U T 


*DB °* 2 X DB 



— I 

DR =X DB -X 

DL=X+^X db 

DR ' = DR + i DB 

DL ' = DL + & DB 

DR" = (X + DR') (jU D 

h /<!° cl - ' u dI 

\<cd<|u c ! + ! 


) + CD(|U C | - |U D |) 


D ’ |U D' 


■ j| + DR ' 
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Table 6-2 Control Actions in Small-Error Control Regions 


Region 

I 

II 

III 

IV 
HS_I 


HS-II 


CS_I 

CS-II 


Control Action 

Drive state to BOUNDARY- I_TQ_COAST 

Drive state to TARGET- DISTURB 

Drive state to zero rate 

Drive state to BOUNDARY- IV_TO_COAST 

If coming from LARGE ERROR CONTROL 
LOGIC or REGION III then drive state to 
zero rate, else coast 

If coming from LARGE ERROR CONTROL 
LOGIC or REGION II then drive state to 
TARGET-DISTURB, else coast 

coast 

coast 


The deadband increment A DB is added to prevent fuel- consuming overshoot when 
the disturbance acceleration is zero or very small. The points 

DR'=DR + A DB (6-35) 

DL' = DL + A db (5-36) 

define the placement of the curves BOUNDARY-L.TO— COAST and BOUNDARY— I V_ TO- 
COAST. 

The factor CD determines the weight given to response time as opposed to fuel 
consumption, by adjusting the steepness of the switch curves that are normally crossed 
when entering the coast region. To insure that BOUNDARY_IV_TO_COAST is always 
less than BOUNDARY_II_TO_COAST and BOUNDARY_I_TO_COAST is always greater 
than BOUNDARY_III_TO_COAST, CD must be between 0 and 1. Currently the value 
of CD is 0. 8. 

Hysteresis has been incorporated to make the most efficient use of disturbance 
accelerations, to prevent switching line chatter, and to prevent tight limit cycling about 
the point DR. The dashed curves within the shaded area and parallel to BOUNDARY_ 
II_TOJCOAST provide hysteresis in the regions HS_I and HS_II. If the state lies 
inside either Region HS_I or HS_II and the sign of the newly computed change in rate 
is not the same as the sign of the last computed change in rate, then no control is ap- 
plied and the vehicle is allowed to coast. However if the signs are the same, then in 
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Region HS_I the state is driven to zero and in Region HS-II the state is driven to 
TARGET-DISTURB. Whenever the state lies inside either Region CSU or CS_II the 
vehicle is simply allowed to coast. 

When the state lies to the right of BOUNDARY_I-TO_ II position error is reduced 
more quickly by driving the state to BOUNDARY— I _TO—COAST and then coasting instead 
of driving the state to TARGET-DISTURB and then coasting. This is particularly im- 
portant as the disturbance acceleration approaches zero. The curve BOUNDARY_I— TO_ 
II is defined by the intersection of TARGET-DISTURB and BOUNDARY_I_TO_COAST. 

When in Region III the rate error is driven to zero, and when in Region IV the 
state is driven to BOUNDARY— IV— TO-COAST to more quickly reduce position and rate 
errors. 

6. 2. 3. 3 Large-Error Control Law 

Control regions for relatively large errors are illustrated in Fig. 6-7, defini- 
tions of the boundaries are summarized in Table 6-3, and corresponding control ac- 
tions are given in Table 6-4. As in the preceeding subsection, all logic is formulated 
in terms of positive disturbing accelerations. In the actual implementation (see Sub- 
section 6.4), negative disturbances are treated by inverting the signs of the state var- 
iables, 0 and to . 

e e 

The purpose of the large-error control logic essentially is to drive the state into 
the small-error control region. Bounds between the small-error and large-error 
regions are set by the values of the parameters RL and RL' (currently 4. 0 and 3. 0 
deg/sec respectively). The rate error bounds are equal to ± RL. The left position 
error limit, -PLL, is defined by the intersection of BOUNDARY-IV- TO-COAST and 
COAST-IV-BOUNDARY. Similarly, the right position error limit, +PLR, is defined 
by the intersection of BOUNDARY_I_TO_COAST and COAST-I-JBOUNDARY. The 
boundaries COAST_IV_ BOUNDARY and COAST_I_BOUNDARY are rate errors equal 
to ± RL ' . 

TARGET-ARL is the average of RL and RL'. TARGET_ARL and -(TARGET- 
ARL) are taken to be the target rate when the state is in Regions II and I respectively; 
i. e. , when the firing is against the disturbance acceleration. When the state is in 
Regions III and IV the targets are the nearest boundary; i.e. COAST-IV— BOUNDARY 
and -RL respectively. 

Hysteresis is incorporated into the large-error coast regions to prevent chatter 
about +RL and COASTLI-BOUNDARY. If the state is driven into Region HLd from 
Region II, then it is further driven to TARGET-ARL. If the state otherwise appears 
in Region HL_I the vehicle is allowed to coast. Similarly, if the state is driven into 
Region HL II from Region I then it is further driven to -(TARGET_ARL); otherwise 
the vehicle is allowed to coast. Should the state appear in either Region CL_I or 
Region CL_II the vehicle is simply allowed to coast . 
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REGION II 


REGION HL_I 
REGION CL_I 
COAST IV BOUNDARY 


TARGET ARL"V^ 


+ RL 


REGION I 


CT! 

I 

to 

to 


-PLL 


SMALL ERROR 
CONTROL REGION 


+ PLR 


0e 


REGION III 


- (TARGET_ARL) 


-RL 




COAST_l_BOUNDARY 

REGION H L 1 1 

REGION CL II 


REGION IV 


Note: DRAWN FOR DISTURBANCE 

ACCELERATION , U Q , POSITIVE. 


Fig. 6-7 


Large Error Control Region 





Table 6-3 Definitions of Boundaries for the Large-Error Control Regions 


Boundary/ Parameter 

Definition 

COAST_IV_BOUNDARY 

co e = RL' 

COAST— I— BOUNDARY 

“ e = - RL ' 

TARGET-ARL 

0J e = (RL + RL')/2 

PLR 

0 = + ^ RL ^ + DR ' 

2CD(|U C | + |U D |) 

PLL 

6 = - ^ RL * - DL' 

0 


2CD(|U C | - |U D |) 


Table 6-4 Control Actions 

in Large-Error Control Regions 

REGION 

CONTROL ACTION 

I 

Drive state to -(TARGET-ARL) 

II 

Drive state to TARGET-ARL 

III 

Drive state to COAST-IV- BOUNDARY 

IV 

Drive state to -RL 

HL_I 

If coming from region II then drive state 
to TARGET ARL, else coast 

HL_II 

If coming from region I then drive state 
to -CTARGET_ARL), else coast 

CL_I 

coast 

CL_II 

coast 
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A final feature of the small-error control logic can be described, now that 
TARGET_ARL has been Introduced. If the disturbance acceleration Is very large, 
part of the curve TARGET JDISTURB may lie below - (TARGET_ARL). In this situa- 
tion a check is made to see if a rate change to drive the state to TARGET_DISTURB 
is greater than a rate change to drive the state to -(TARGET_ARL). If it is, the 
state is driven to -(TARGET_ARL), which means less fuel consumption with the re- 
sult of finishing in a limit cycle about the origin. 

6. 3 On-Orbit TVC Phase: Preliminary Functional Summary 

The Orbiter Maneuvering Subsystem (OMS) is used for major on-orbit maneuvers 
and to establish the deorbit entry trajectory. During powered flight the roll, pitch 
and yaw of the space-shuttle orbiter is achieved through deflection of the two gimbaled 
engines. Computation of gimbal servo commands in response to computed errors be- 
tween commanded and measured attitudes is the function of the on-orbit TVC (Thrust 
Vector Control) "phase" of the DFCS. The following functional summary is an adapta- 
tion of the Apollo Command Module TVC control; it is illustrated in Fig. 6-8. 

1) The steering-loop computations generate incremental attitude commands in 
inertial coordinates and transform them into body coordinates. 

2) The DFCS back-differences the gimbal angles to obtain incremental variations 
over each sampling interval. 

3) The gimbal increments over each sampling interval are then transformed into 
body coordinates and subtracted from the commanded increments which are 
generated by the steering program. 

4) The respective attitude errors are fed to the roll, pitch and yaw compensa- 
tion filters, whose outputs contribute to the commands to the engine gimbal 
servos. These signals bear the burden of controlling rapid transients, deter- 
mining the system bandwidth and the stability of the vehicle slosh and bend- 
ing modes. 

5) As shown in Fig. 6-8b, the total command signal to each engine-gimbal 
servo is made up of the component from the associated compensation filter 
plus another component from a Thrust Misalignment Correction (TMC) loop. 
The latter component serves to bias the total command so that a zero output 
from the compensation filter will cause the thrust vector to pass exactly 
through the center of gravity (eg) when there is no eg motion and no motion 
of the thrust vector relative to the commanded angle. 
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Fig. 6 -8b Block Diagram of the Thrust Misalignment Correction (TMC) Loop 





6) The two major elements of the TMC loop are its summing register, which 
supplies the bias, and its low-pass filter, DELFILTER, which tracks the 
total command signal at autopilot sampling frequencies. The difference 
between the bias and the DELFILTER output is used to increment the sum- 
ming register so as to slowly correct for the offset caused by thrust-to-cg 
misalignment. This slow reset action is roughly equivalent to inserting a 
proportional-plus-integral transfer function between the compensation fil- 
ter and the total command signal. The TMC loop is designed so that its 
dynamics have no effect on the vehicle bending and slosh modes and very 
little effect on the rigid-body stability. 


7) The operation of the TVC digital autopilot cannot be considered separately 
from that of the steering loop which interacts with the autopilot. The 
steering loop operations are shown in Fig. 6-8 for the case of "External AV" 
guidance. This form of guidance is based on achieving a desired acceler- 
ometer-measured velocity change, AV^-,, which is inserted into the steer- 
ing program as an initial value of the computed velocity-to-be-gained vec- 


tor, V g . 


The V g vector is periodically reduced during the course of the 


burn by subtracting the accelerometer-measured increments, AV. 


8) The purpose of the steering loop is to align the vehicle thrust with the 

velocity-to-be-gained. This is achieved by commanding a vehicle turning- 
rate which is proportional to the normalized cross-product between these 
two vectors. 


9) The pitch- and yaw-axis components of the vector cross-product are mul- 
tiplied by the proportionality constant K g(;eer to obtain the commanded 
rates about these axes - and multiplied by the autopilot sampling period, 

Tg, to obtain attitude -command increments, which are supplied to the 
autopilot every Tg seconds. 

10) The roll control channel is designed to provide independent attitude and 
rate control about the roll axis by differential control of the OMS engines. 

6. 4 Digital Implementation (ACPS phase) 

The overall composition of the DFCS as it applies to the on-orbit phase is re- 
iterated in Fig. 6-9 in structured flow-chart form. The procedure MA JOR_TRANSIENT 
provides initial mode and flag information to be used by all mission phases. Currently, 
no initial mode or flag information unique to the on-orbit ACPS phase has been im- 
plemented in MAJOR_TRANSIENT. The procedure INPUT-INTERFACE contains, in 
general, information that is explicitly read by the DFCS from external sources. More 
specifically, when the on-orbit phase is selected, the gimbal angles are read and the 
time of their reading, t n> is recorded. On the first pass of the DFCS for a given 
phase FILTER. AND- PARAMJNIT and PARAM are both called. FILTER- AND— PARAM — 
INIT currently does nothing in this phase. PARAM calls MASS_CALC (See Fig. 6-10) 
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CALL CONTROL 


1 



REFERENCED^) 

PHASE. [tt e 3. [e e ], [« c ], 

C U £)3» [iu] 




ASSIGNED ^ 

[A Ui] 

J 


CALL PART 2_FILTER 




S3 

^ R EFER ENCED^> 

2cu* -CV* 2cw* [t JET^‘ £ SGNC 3 

W* S' Srn* 2* 0 B 




SI 

ASSIGNED ^ — 

1'b- &'■ 2 





ISTART * 0 
IS = IS + 1 


IM * IM + 1 

IM = IM - 1TM INTEGER (IM/ITM) 
IS = IS + 1 

IS = IS - ITS INTEGER (IS/ITS) 




mmm 

REFERENCE!^ 

PHASE 

1 




ASSIGNED ^ 

null 

K 





Pig. 6-9 Logic Flow of DFCS 

















RETURN 


Fig. 6-10 MASS CALC 
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which then takes as an input the torques supplied by each ACPS thrust group and com- 
putes three more or less independent control axes and a pilot -to- control transformation 
matrix. The procedure PART1_FILTER shown in Fig. 6-11 is then entered which 
takes as an input the extrapolated estimates of angular position, rate, and disturbance 
acceleration for each axis and incorporates the measurement of gimbal angles taken 
at time t n . From these updated estimates the procedure CONTROL (Fig. 6-12) calls 
upon: 


1) ERRORS (Fig. 6-13) to compute the errors in attitude and angular rate 
from the desired 

2) PHASE-PLANE (Fig. 6-14) to compute the desired change in angular rate 
to correct these errors 

3) JETjCONTROL (Fig. 6-15) to select the appropriate .jets and control 
their firings, according to provisional logic suggested in Subsection 6. 2. 1. 


The digital implementation of the phase plane switching logic is composed of 
three principle routines: 

PHASE-PLANE 

The procedure PHASE— PLANE selects the appropriate sign to be attached to the 
computed change in angular rate, computes position error intercepts, and also com- 
putes the boundaries between the small error control logic and the large error control 
logic, as illustrated in the logic flow diagram in Fig. 6- 14a. PHASE-PLANE then 
selects either the small-error control logic or the large-error control logic to com- 
pute the desired change in rate. 

LARGE-ERROR-CONTROL-LOGIC 

The procedure LARGE-ERROR-CONTROL-LOGIC computes the change in rate 
to drive the state into the small error control region. As illustrated in Fig. 6- 14b, 
the region in which the state lies is singled out and then the appropriate control action 
is computed. 

SMALL-ERROR-CONTROL-LOGIC 

The procedure SMALL-ERROR— CONTROL-LOGIC computes the change in rate 
to drive the state into a limit cycle with specified deadband about the desired position 
and rate. Current testing has used deadbands of 0. 5 and 5.0 degrees. As illustrated 
in Figs. 6 -14c and 6-14d, the region in which the state lies is singled out and then the 
appropriate control action is computed. 

The function BOUNDARY, illustrated in Fig. 6-14e, is implemented to compute the 
various switching curves. 
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Fig. 6-11 



6 




-FILTER Procedure 










Fig. 6-12 On-orbit Segment of CONTROL Procedure 
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Fig. 6-13 ERRORS 
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IF 6 a +PLR AND u a COAST_I_BOUNDARY 
C e ^ 


X D = “ e ' 


(COAST-I-BOUNDARY - RL) 


IF 0 g £ +PLR AND w e 2 +RL 


IF 0 e £ -PLL AND u) e £ COAST-IV- BOUNDARY) — H — 


IF 6 a -PLL AND u £ -RL 
e e 



IF X D old^°-° 
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Fig. 6- 14c SMALL-ERROR CONTROL LOGIC 
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Fig. 6-15a JET_CONTROL Procedure 
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Inputs: otQ, Uj-, 



Fig. 6- 15b UVWTIME Procedure 
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Fig. 6- 15c OFF1JETS Procedure 



Fig. 6-15d OFF2UJETS Procedure 
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Fig. 6-15g OFF1WJETS Procedure 



Fig. 6-15h OFF2WJETS Procedure 
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After CONTROL has called JETCONTROL to Implement the firing of various 
jets, the estimate of the state is extrapolated for the next pass of the DFCS by the 
procedure PART2_FILTER (see Fig. 6-16). Finally, PARAM is called on passes 
other than initialization passes, and the DFCS pass is terminated. 
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Table 6-5 Glossary for On-Orbit Phase 
MNEMONIC ENGINEERING DEFINITION 



TERM 


ACC2U 

^CU 

Angular acceleration about u-axis produced by Mp 9 ^ or 

ACC2V 

-cv 

Angular acceleration about v-axis produced by Mp 2 ^ 'or' M M2V 

ACC2W 

-cw 

Angular acceleration about w-axis produced by M p2w or M M2W 

ALPHADBAR 

A 

s'd 

Extrapolated disturbance acceleration 

ALPHADHAT 

A 

“d • 

Estimated disturbance acceleration 

ALPHAM 

a 

Bias acceleration 


— m 


basTsu 

-BASIS \ 


basTsv 

-BASIS | 

Unit vectors along u, v, w control axes 

basTsw 

—BASIS ) 


[BIAS-BETA] 

^BIAS-' 

Bias body angles 

CD 

CD 

Scale factor used in phase plane switching logic 

[ CDU] 


Gimbal angles 

[CDUERROR] 

[0 GA e ^ 

Error in gimbal angles from the reference desired 

DELTA 

a db 

Deadband factor 

DELTA-POSITION 

e ' 

e 

Angular position error (within PHASE_PLANE) 

DELTA-RATE 


Angular rate error (within PHASE_PLANE) 

[DELTA_THETA_ 


- 


D] Ue D ] Incremental desired reference glmbal angles 

DELWDM Aw d m Change in angular' velocity vector due to disturbance and bias 

accelerations 


DELWU Aw^ Change in angular rate due to u- control jets 

DELWV- Aoj v Change in angular rate due to v-control jets 

DELWW Aw. w Change in angular rate due to w-control jets 

DL, DL-PRIME1, DL,DL' \ 

DR, DR-PRIMEl, DR, DR' > Switching line angular position error intercepts 
DR-PRIME2 DR " ) 

[DRATE] [Aw] Components of changes of angular rates 
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MNEMONIC 

Table 6-5 Glossary for On-OrbLt Phase (Cont) 
ENGINEERING DEFINITION 


TERM 


DTC 

DTC 


Total firing time of an ACPS control jet for a maneuver 

ISTRTR 

ISTRTR 

Restart required flag (external) 

PHASER 

PHASER 

Mission phase flag (external) 

[FAST] 

[MEDIUM] 

[FAST] 

[MEDIUM] 

Fast, medium, and slow sampling periods for each 
mission phase 

[SLOW] 

[SLOW] 

Slow sampling period for each mission phase 

GAIR 

G AIR 

Filter gain for air data 

inIrtia 

J 


Inertia matrix 

IM 

IM 


medium rate sampling index 

IS 

IS 


Slow rate sampling index 

ISTART 

ISTART 

Initial pass flag for mission phase 

ITM 

ITM 


Ratio of medium to fast sampling periods 

ITS 

ITS 


Ratio of slow to fast sampling periods 

[ JTORQUE ] 

[JTORQUE] 

Array of ACPS jet torque vectors 

K 

K 


Sign of rate change 

KALPHA 

K a 

i 


KOMEGA 

K ^ 
« ( 

r 

> 

Filter gains 

KTHETA 

1 

K e i 

i 


KILL1 

KILLl 

ACPS jet continuity indicator 

[MAGACC] 

[a c J 


Array of magnitudes of control accelerations 

MAG_U_D 

|U D | 


Magnitude of disturbance acceleration 

MAG_U_C 

|u c l 


Magnitude of control acceleration 

MP2U 

-P2U 

Torque in the +u direction by 2 u-jets 

MM2U 

^M2U 

Torque in the -u direction by 2 u-jets 

MP2V 

Mp2V 

Torque in the +v direction by 2 v-jets 

MM2V 

M 

-M2V 

Torque in the -v direction by 2 v-jets 
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Table 6-5 Glossary for On-Orbit Phase (Cont) 


MNEMONIC ENGINEERING 

TERM 

DEFINITION 

MP2W 

M 

— P2W 

Torque Ln the +w direction by 2 w-jets 

MM2W 

M 

— M2W 

Torque in the -w direction by 2 w-jets 

m£p 

]& C p 

Control to pilot transformation matrix 

mI>c 

1&PC 

Pilot to control transformation matrix 

MIN IMP 

fc IMP 

Time of a minimum impulse ACPS jet firing 

NETACC 

[a] 

Net acceleration 

OMEGABAR 

A 

w' 

Extrapolated angular velocity 

OMEGAHAT 

A 

a) 

Estimated angular velocity 

-PLL 

-PLL 

Left angular position error limit 

+PLR 

+PLR 

Right angular position error limit 

P—CEPT 

P 

cept 

Intermediate position error axis intercept 

PHASE 

PHASE 

Phase of the mission 

±RL 

±RL 

Upper and lower rate error limits 

POSITION_ERROR 


Angular position error 

RATE_ER ROR 

tv 

Angular rate error 


REALSEC 

REALSEC 

Scale factor for scaling computer time (= ] 

\ TICK_TO_SEC / 

RUNTIME 

RUNTIME 

Runtime of the computer in computer ticks 

[SGNC] 

[SGNC] 

Signs of the angular rate components 

TARGET. ARL 

TARGET.ARL 

Lower rate limit average 

TF 

tfast 

Fast sampling period 

TM 

'"medium 

Medium sampling period 

TS 

'"slow 

Slow sampling period 

TEMPI 

TEMPI 

Temporary vector of body angle errors in control coordinates 

TEMP2 

TEMP2 

Temporary vector of body rate errors in control coordinates 

TEMP3 

TEMP3 

Temporary vector of acceleration in control coordinates 

[TJET] 

t JET 

Time of an ACPS jet firing in seconds 
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Table 6-5 Glossary for On-Orbit Phase (Cont) 


MNEMONIC ENGINEERING 

TERM 


DEFINITION 


thetabar 

Ib 

Extrapolated body angles 

THETAHAT 

A 

-B 

Estimated body angles 

TNOW 

*n 

Time of present glmbal angle measurement 

[TJETQ] 

t JETQ 

Time of an ACPS jet firing scaled to the computer 

RATE-IN T 

-B 

"Measured" body angles 

[TLAST] 

t, 

1 

[0 D ] 

Time of the last glmbal angle measurement 

[THETA-D] 

Desired set of glmbal angles 

TICK_TO_SEC 

TICK_TO_SEC 

Scale factor for scaling computer time 

[U_D] 

[u D ] 

Array of disturbance accelerations 

W_D 


Desired angular velocity 

X-BAR 

X 

Average angular position error 

X-DB 

X DB 

Deadband of limit cycle 

X-D 

X D 

One component of Au) the desired change in angular velocity 

NO-CHANGE 

X D old 

One component of the last desired change in angular velocity 

XYZERROR 

-Be 

Errors in x, y, z axis rotation angles 


NOTATION CONVENTIONS 

© 

Modular addition 

o 

/ 

Modular subtraction 

/*>-* 

Measured quantity 

A 

Estimated quantity 

/ 

Extrapolated quantity 

— 

Vector 

J e 
*)' 

Matrix 

[ ] 

Array 


6-47 



REFERENCES FOR SECTION 6 


1. Keene, D. , "Control Axes for the Shuttle On-Orbit Autopilot", SSV MEMO 
No. 73-23C-6, April 1973. 

2. Keene, D. , "State Estimator for the Shuttle On-Orbit Autopilot", SSV MEMO 
No. 73-23C-7, April 1973. 

3. Keene, D. and Kubiac, E., Section 6, "On-Orbit Phase", Integrated Digital 
Flight Control System for the Space Shuttle Orbiter , E-2736, March 1973. 


6-48 



SECTION 7 


BOOST PHASE 
by 

Bernard Marcantel (JSC/EJ3) 


7. 1 Introduction 

The baseline launch configuration control system software described in this 
section satisfies the requirement to control the booster from lift-off through main 
engine cutoff at external tank separation. The system is adaptable to various control 
laws, including attitude control, drift minimum control (with and without accelerome- 
ter feedback) and/or load relief control through proper choice of gain schedules. It 
also contains the option of using rudder and thrust vector control combined in the roll 
channel . 

The following paragraphs contain: A description of the control system with 
functional block diagrams of the pitch, yaw, and roll channels; an annotated flow dia- 
gram; and a discussion of the control gains. 

7.2 Functional Description 

The baseline launch control system is a digital autopilot receiving body attitude 
from an inertial reference system, body rates from three sets of body mounted rate 
gyros, and translational accelerations from two sets of body mounted accelerometers. 
The current baseline includes rate gyros at the aft end of each solid rocket booster 
and the orbiter. Normal and lateral accelerometers are located in the forward avionics 
bay and the aft end, of the orbiter. The autopilot also utilizes prestored engine deflec- 
tion commands and acceleration commands. The former are the nominal trim values 

for the expected center-of-gravity and aerodynamic histories, and the latter are the 

7- 1 

normal linear accelerations in the nominal trajectory. The final inputs to the 
control loop are the attitude commands, which may either be prestored or come from 
a guidance system. 

From the inputs, the control system determines and sends out deflection com- 
mands for the engines and the rudder. A general optimal mixing logic has been de- 
signed and may be installed at a later date. It has been designed to distribute the de- 
sired torques efficiently among the control actuators while accounting for failures and 
7-2 

saturation. 

Table 7-1 summarizes the control system inputs and outputs. 
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Table 7-1 Input/Output Summary 


INPUTS 


• Inertial attitude angles 

• Body rotational rates 

• Y & Z translational accelerations 

• Prestored engine deflection commands 

• Prestored accelerations commands 

• Attitude commands from guidance system 


OUTPUTS 


• Deflection commands to orbiter engines 

• Deflection commands to SRM engines 

• Deflection commands to aerosurfaces (rudder) 

Variable control gains have been calculated to satisfy the rigid body response 
requirements. They are stored as a function of time from lift-off; the numerical 
values are presented in Subsection 7. 4. In the future, the design will be extended to 
include consideration of slosh and bending, and these prestored gains will be replaced 
by filters. All filters will be implemented as second order over second order. Fil- 
ter inclusion in two places in each sensor loop recognizes common filtering require- 
ments and reduces the order of the filter. 

Some aspects of the design presented in this section are specific to the first 
stage of boost (prior to SRB separation). However, control of the orbiter during the 
second stage will utilize basically the same control law. The major differences are 
that the rudder will not be used and there will be no acceleration feedback in the sec- 
ond phase. Furthermore, there will only be one set of body mounted rate gyros 
available. These gyros will be located in the aft section of the orbiter. All filter 
coefficients and gains will be changed at separation and vary as the mass properties, 
flexible body and fuel slosh characteristics vary. 

7.2.1 Roll Channel Description 

Referring to Fig. 7-1, the block in the lower right hand corner labeled "Vehicle 
and Subsystem Dynamics" includes actuator and sensor dynamics, possible analog 
prefiltering of sensor outputs, input/output voting, redundancy management, A/D and 
D/A converters. Inputs to the control system from "Vehicle and Subsystem Dynamics" 
are the inertial platform attitude angle, and rate measurements from three sets of 
body mounted rate gyros. Outputs from the control system to the vehicle are deflection 
commands for the engines and the rudder. 


7-2 



7-3 



Fig. 7-1 Baseline Launch Configuration Control System Software Roll Channel 














Following the attitude loop in Fig. 7-1, the platform gimbal angle, 0, is differ- 
enced with the inertial attitude command, 0^, , to produce the inertial attitude error, 

0-gj. This inertial attitude error is transformed to a body roll error, 0gg, by adding 
the inner gimbal error times the sine of the middle gimbal angle (^ E g = S* EI + 0^ sin <fi). 

It should be noted that a limiter has been included in the attitude error loop. 

This limiter reserves some TVC authority for rate damping (even though the attitude 
error signal may be saturated) thus preventing oscillation and complete loss of con- 
trol because of large signal instability. The limiter is set to a value less than eight 
degrees. 

Returning to the "Vehicle and Subsystem Dynamics" block in Fig. 7-1, three 
sets of roll rate gyro signals are filtered then blended and again scaled and filtered 
prior to being summed with the limited attitude error signals. The resulting signal 
is sent to .the rudder and engine mixing logic for distribution to the individual engine 
and rudder actuators. Then each actuator command is limited to eight degrees to 
prevent the roll channel from usurping all pitch and yaw controllability. The rudder 
command is limited to 15 degrees. 

7,2.2 Yaw Channel Description 

The yaw channel block diagram in Fig. 7-2 is laid out in the same manner as 
the roll channel with the "Vehicle and Subsystem Dynamics" block in the lower right 
hand corner, the "inertial Attitude Command" in the upper left hand corner and the 
"Actuator Mixing Logic" on the right. The rudder is not used for yaw control, be- 
cause there is adequate engine control authority to handle yaw-plane torque. 

Three features which were not included in the roll channel appear in the yaw 
channel. Two accelerometer inputs have been added to the attitude loop, a forward 
loop integrator follows the limiter, and a rate command signal is differenced with 
the input from the rate gyros. 

The platform gimbal angle, , is differenced with the inertial attitude com- 
mand, , transformed to body attitude errorC^g = ^ EI cos 0 - 0 EI sin 0 cos (i) 
and scaled. The attitude error signal is then modified by summing with an accelera- 
tion error signal and limiting. The acceleration error is obtained by blending appro- 
priately filtered signals from two body mounted accelerometers and then applying gain 
and filtering compensation according to the control law or bending requirements. 

The accelerometer gain is varied with time to meet the drift minimum or load relief 
requirements. The forward loop integrator trims out steady state errors resulting 
from biases and thrust misalignments. 

The rate command signal is introduced because the normal attitude command 
changes at a predictable rate. The rate command signal minimizes the tendency for 
vehicle attitude to lag behind the attitude command due to rate feedback. If the 
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Fig. 7-2 Baseline Launch Configuration Control System Yaw Channel 























vehicle rolls appreciably the rate command must be applied to the yaw channel. In 
operation, the attitude command Is differentiated, multiplied by the sine of 0, and 
differenced with the rate feedback signal. The rate error signal is then added to the 
attitude error signal and sent to the actuator mixing logic routine. Actuator command 
limiting is performed for the same reason stated for roll control. 

7.2.3 Pitch Channel Description 

Comparison of Fig. 7-3, which is the pitch channel block diagram, with Fig. 7-2 
(yaw channel) reveals that the only difference is the addition of prestored engine de- 
flection and acceleration commands. Both of these command schedules are the nom- 
inal values from a reference trajectory. 

The addition of the prestored acceleration command permits use of accelerom- 
eter feedback with any feasible trajectory to modify the attitude error signal to con- - 
form to the desired control law. 

The prestored engine deflection command schedule provides coarse control of 
the engine gimbal trim in the pitch plane. Center of gravity travel alone could be 
handled by the integral compensation, but aerodynamic pitching moment increases 
sharply in the vicinity of Mach one and decreases again thereafter. A high integrator 
gain would be required to follow these rapidly changing deflection requirements. 

This is undesirable because it degrades the rigid body stability margin; therefore, 
the prestored commands are used to control the nominal part of this variation. 

The rest of the pitch channel is the same as the yaw channel. All engine com- 
mands are sent to the actuator mixing logic routine. 

7. 3 Logic Flow Description 

A flow chart showing the logic flow of the control system is presented in Fig. 7-4. 
7. 4 Values of Control Gains 

The basic control system is independent of configuration, assuming that thrust 
vector control is primary for pitch and yaw, and that the principal guidance commands 
will be in the pitch plane. It includes no provision for use of elevons during boost. 

The control gain and command schedules are, of course, configuration depend- 
ent . They will be affected by modifications to either the vehicle configuration or 
the reference trajectory. This will not impact the design, however, because all con- 
trol gains and commands are continuously variable input functions of time. Some of 
them may subsequently be modified to functions of other trajectory variables, such 
as altitude and velocity, but this will involve only minor software modifications. 

Table 7-2 presents the current numerical values of the control gains. 
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Fig. 7-3 Baseline Launch Configuration Control System Software Pitch Channel 
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Fig. 7-4 Logic Flow (Cont) 
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Fig. 7-4 Logic Flow (Cont) 
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6 4P " APP 4 ® PITCH + ARP 4 6 ROLL + 5 SRM 

6 5P * APP 5 ® PITCH + ARP 5 ' 6 ROLL + 6 SRM 
6 4Y * AYY 4 6 y aw + ARY 4 SpoLL 

«5Y * AYY 5 6 YAW + ARY 5 6 RQll 
64R ■ 0. 707 (64P - 

. 6 5R " °* 707 < fi 5P ' 6 5Y> 

64L * 0. 707 (64Y + 64 P ) 

«5L • 0.101 («- 5Y + e 5 p) 

Limit to ± 8*: 6 4R . 6 5R . 6 4L '. 6 5L 

( END j 



Fig. 7-4 Logic Flow (Cont) 







Table 7-2 Gains for Boost Phase 


Control Gains vs Time ( sec ) 


0. 

25 

30 

■ 75 

95 

120 

150 

Kg 1. 1 

1. 1 

. 88 

. 88 

. 66 

. 66 

. 66 

1.4 

1. 4 

. 56 

. 56 

1. 4 


1. 4 

K ^ .707 



. 707 

. 424 


. 424 

k az 0 

0 

-. 109 

-. 109 

0 


0 

k ay 0 

0 

. 137 

. 137 

0 


0 

Kp .707 



. 707 

. 424 


. 424 

Kq 1.1 



1. 1 

. 66 


. 66 

K r 1.4 



1.4 

. 84 


. 84 

Kp.. Kq., Kr. j 
K Yp K Z . | 

k pa k ya 

s 

1 . 

Currently unity for 

all values of time. 




Mixing Logic Gains 



Orb iter Engines 


APP 

AYY 


ARP 

ARY 

1 


1.4 

0 


0 

-3. 5 

2 


1. 4 

1. 4 


-3. 5 

1. 8 

3 


1. 4 

1. 4 


3.5 

1.8 

SRM Engines 







4 


1. 4 

1.4 


-1.4 

0 

5 


1. 4 

1.4 


1. 4 

0 


ARR = -7.0 


Miscellaneous Fixed Gains 

PLIM 7. 00 degrees 

YLIM 7. 00 degrees 

RLIM 5. 00 degrees 

t • 0.1 
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SECTION 8 


DFCS GLOSSARY 


Phase 

Mnemonic 

Symbol in Text 

Definition 

4 

A 1,2, 3,4 

a l* a 2' a 3’ a 4 

Gains for aileron control 

4 

A 11, 12,21 ,22 

a ll' a 12' a 21’ a 22 

Lateral fixed gains 

3 

ABC ( 1 ) 

d ROLL 

Total roll channel command 

8 

ABC (2) 

6 PITCH 

Total pitch channel command 

8 

ABC (3) 

6 YAW 

Total yaw channel command 

8 

ABE (1 ) 

0 EB 

Body roll attitude error 

8 

ABE (2) 

®EB 

Body pitch attitude error 

8 

ABE (3) 

^EB 

Body yaw attitude error 

8 

AB EF ( 1 ) 

**F 

Scaled roll error 

8 

ABEF (2) 

0 F 

Scaled pitch error 

8 

ABEF ( 3) 

*F 

Scaled yaw error 

2,4 

ACCBOD 

~b 

a 

Linear accel measurement (body coord) 

4 

ACCG 

g 

Acceleration of gravity 

4 

ACCINT (1) 

Ap-* 

Roll integrated angular acceleration 

4 

ACCINT (2) 


Pitch integrated angular acceleration 

4 

ACCINT (3) 

Ar J 

Taw integrated angular acceleration 

4 

ACCL 

a L 

Measured lift acceleration 

4 

ACCLMN 

a LM 

Mean lift acceleration 

4 

ADOT 

• 

a 

a 

Derivative of angle-of-attack wrt time 

4 

ADO (1 ) 

La 

Angle-of-attack rate filter output 

4 

ADO (2) 

La' 

Angle-of-attack rate filter extrapolated value 

2,4 

AG 

a G 

Angle-of-attack guidance command 

4 

AGOLD 

a old 
G 

Previous-pass value of AG 

4 

AI (1) 

Aa 

Angle-of-attack filter input 

2 

AINT 

a INT 

Aileron integrator output 
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DFCS GLOSSARY (Cont) 


Phase 

Mnemonic 

Symbol in Text 

Definition 

2 

AINTI 


Aileron integrator input 

2 

AINT1 


Aileron integrator node 

2 

AINTMX 

max 

a INT 

Aileron integrator output limits(+--max) 

8 

AI V (1 ) 


Platform outer girabal angle 

8 

AT V (2) 

e 

Platform inner gimbal angle 

8 

AI V (3) 

0 

Platform middle gimbal angle 

8 

AIVC ( 1) 


Platform outer gimbal angle command 

8 

AIVC (2) 

0 C 

Platform inner gimbal angle command 

8 

AI VC ( 3) 

0 C 

Platform middle gimbal angle command 

8 

AI VE ( 1 ) 

0 EI 

Platform outer gimbal angle error 

8 

AI VE ( 2) 

0 EI 

Platform inner gimbal angle error 

8 

AIVE(3) 

^EI 

Platform middle gimbal angle error 

4 

ALIM 

°W 

Limited value of angle-of-attack filter output 

2 

AH 

a 

Measured angle-of-attack 

4 

AMEAN 

“t 

Mean (trim) angle of attack 

2 

AN 

an 

Normal acceleration (-zbody, with gravity) 

2 

ANF 

an F 

Normal acceleration filter output 

2 

ANGRV 

an GRV 

Gravity corrections for AN 

2 

ANN (I) 


Normal acceleration filter nodes 

4 

AO(1) 

Aa 

Angle-of-attack filter output 

4 

AO (2) 

Aa' 

Angle-of-attack filter extrapolated value 

8 

APP(I) 

APP. 

i 

Gains for scaling pitch input to the pitch 




actuators 

8 

ARP (I) 

ARP. 

1 

Gains for scaling roll input to the pitch 




actuators 

8 

ARR 

ARR 

Gain for scaling the coll input to the rudder 

8 

ARY (I) 

ARY. 

L 

Gains for scaling roll input to the yaw 


actuators 
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DFCS GLOSSARY (Cont) 


Phase 

Mnemonic 

Symbol in Text 

Definition 

4 

AX, AY , AZ 

U ,U ,U 
x’ y' z 

Jet acceleration in roll, pitch, yaw 

2 

AYM 

ay 

Measured lateral acceleration 

2 

AY MF 

ay F 

Lateral acceleration filter output 

2 

AYMN (I) 


Lateral acceleration filter nodes 

8 

AYY(I) 

AYY. 

L 

Gains for scaling yaw input to the yaw 
actuators 

2 

azh 

az 

Measured normal acceleration (♦zbody) 

4 

B 

b 

Fixed estimator parameter 

4 

B 11,12,21 

,22 b^ 1* ^ 12 * ^21* ^^2 

Lateral fixed gains 

4 

BANK 


Boll (bank) about air-relative velocity vector 

4 

BANKDO (1) 

4 V 

Bank rate filter output 

4 

BANKDO (2) 


Bank rate filter extrapolated value 

4 

BANKG 

Kg 

Bank guidance command 

4 

BANKI (1) 

. 

L( t>v 

Bank filter input 

4 

BANKH 

Km 

Mean bank angle 

4 

BANKO (1) 

4 V 

Bank filter output 

4 

BANKO (2) 

*K 

Bank filter extrapolated value 

4 

BOOT 

0 

Derivative of sideslip wrt time 

4 

BDO (1 ) 

4 

Sideslip rate filter output 

4 

BDO (2) 

4 

Sideslip rate filter extrapolated value 

4 

BDTM 

% 

Hate limit in roll -jet logic 

2 

BE 


Sideslip error filter input 

2 

BEF 

?EF 

Sideslip error filter output 

2 

BEN (I) 


Sideslip error filter nodes 

2 

BG 


Sideslip guidance command 

4 

BI(1) 

tH 

Sideslip angle filter input 
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DFCS GLOSSARY (Cont) 


Phase 

Mnemonic 

Symbol in Text 
. old 

4 

BKGOLD 

4 VG 

2 

BM 


4 

BO (1 ) 

A 

W 

4 

BO (2) 

l}' 

• old 

4 

BRGOLD 

<PVG 

4 

CAT 

COSOrp 

' 2 

CC 

c C 

2 

CCAMN 


2 

CCAMX 


2 

CCCMN 


2 

CCCMX 

min 

2 

CCMN 

C c 

2 

CCMX 

max 

C c 

2 

CE 

C E 

2 

CEF 

C EF 

2 

CEN(I) 


2 

CF 

C F 

2 

CG 

C G 

4 

CLBP 

C ' 
% 

4 

CLDAP 

Sa 

4 

CNBP 

G n 

n fl 

4 

CNDAP 

c' 

2,4 

CRTB 

c b 

r 

4 

CEOD 

C RUD 

2 

CSTAPG 

* 

°G 

4 

DA1,DA2 

6 al* 5 o 2 

- 

DAC 

6a c 


Definition 

Previous-pass value of BANKG 

Measured sideslip 

Sideslip angle filter output 

Sideslip angle filter extrapolated value 

Previous-pass value of bank rate command 

Cos (AMEAN) 

CAS command (QC, AC, or CSTAEC) 
Angle-of-attack values for CCBN 

Angle-of-attack values for CCMX 
CAS values (C* or q) for CCMN 
CAS values (C* or q) for CCMX 
CAS command limit, min (working) 

CAS command limit, max (working) 

CAS error filter input. 

CAS error filter output 

CAS error filter nodes 

CAS feedback (QF, AF, or CSTAEF) 

CAS guidance command (QG or CSTARG) 

Scheduled estimator parameter 

Scheduled estimator parameter 

Scheduled estimator parameter 

Scheduled estimator parameter 

Hef erence-to-body direction cosine matrix 

Rudder crossfeed gain 

CSTAR guidance command 

Deadbands for pitch jet logic 

Aileron command 
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DFCS GLOSSARY (Cont) 


Phase 

Mnemonic 

Symbol in Text 

Definition 

2 

DACF 

8a CF ! 

Aileron command filter output 

2,4 

DACMX 

. max 
8a c 

Aileron command limit ( ♦-max) 

2 

DACN(I) 


Aileron command filter nodes 

2 

DAG 

8a G 

Aileron guidance command 

8 

DAIVC 

®C 

Platform inner gimbal rate command 

- 

DA MAX 

. max 
oa 

Physical aileron limit (♦-max) 

4 

DAT 

‘a, 

Aileron trim estimate 

2 

DATHIM 


Effective aileron trim, for mode change 

4 

DB 1 , D82 

V •« 

Deadbands for roll -jet logic 

2 

DDG 

'»o 

TAEM guidance interface variable 

- 

DEC 

8e c 

Elevator command 

2,4 

DECMN 

min 

6e c 

Elevator command limit, min 

2,4 

DECMX 

max 
8e c 

Elevator command limit, max 

2 

DECO 


Elevator command from speed loop crossfeeds 

2 

DECOB 


BCS version of DECU 

2 

DEG 

6e G 

Elevator guidance command 

- 

DELC 

6EL C 

Left elevon command 

8 

DELL 4,5 

6 4L’ 6 5L 

Commands to SRM left gimbals 

8 

DELOBB 

6 orb 

Prestored orbiter engine deflection command 

8 

DELP(I) 

6 iP 

Commands to orbiter pitch gimbals (i=1,2,l) 

8 

DELP(I) 

6 iP 

Commands to SRM pseudo pitch gimbals (i=4,5) 

8 

DELR 4,5 

6 4R' 6 5R 

Commands to SRM right gimbals 

8 

DELSRM 

6 srm 

Prestored SRM engine deflection command 

8 

DELY(I) 

6 iy 

Commands to orbiter yaw gimbals (i=1,2,3) 

8 

DELY(I) 

6 iY 

Commands to SRM pseudo yaw gimbals (i=4,5) 

- 

DEMIN, DEMAX 

_ min c max 
6 _ ,6 
e 1 e 

Physical elevator lower (upper) limit 

- 

DERC 

6ER C 

Right elevon command 
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Phase 

Mnemonic 

Symbol in Text 

Definition 

4 

DET 

6e t 

Elevator trim estimate 

2 

DETRIH 


Effective elevator trim, for mode change 

2 

DEC 

6f C 

Body-flap command 

2 

DFCF 

6f CF 

Flap command filter output 

2 

DFCPI 

6f CM 

Flap command from Mach profile 

2 

DFCMCH (I, J) 


2xN matrix of (DFCM,MACH) pairs 

2 

DFCMN 

8f“ in 

Flap command limit, min 

2 

DFCMNB 


BCS version of DFCMN 

2 

DFCHX 

, f max 

8f c 

Flap command limit, max 

- 

DFMIN,DFMAX 

cr min rr max 

OI , 01 

Physical flap lover (upper) limit 

4 

DPI STF 

qf 

Blending parameter for elevator use 

4 

DP2STR 

$2 

Blending parameter for aileron use 

4 

DP3STR 

q| 

Blending parameter for rudder use 

4 

DR1 ,DR2 

S 52 Jl* 6 02 

Deadbands for yaw jet logic 

- 

DRC 

6r C " 

Rudder command 

2 

drcmx 

.max 

Sr c 

Rudder command limit (+-max) 

2 

DRF 

6r F 

Rudder loop feedback 

2 

DRG 

6r G 

Rudder guidance command 

- 

DRMAX 

6r max 

Rudder physical limit (♦-max) 

8 

DROLL 

*C 

Platform outer gimbal rate command 

2 

DRTRIM 


Rudder effective trim, for mode change 

- 

DSC 

«s c 

Speed brake command 

2 

DSCF 

«s CF 

Speed brake command filter output 

2 

DSCMN 

, min 
6s c 

Speed brake command limit, min 

2 

DSCMNB 


BCS version of DSCMN 

2 

DSCHX 

. max 
6s c . 

Speed brake command limit, max 

2 

DSCP 


Past value of DSC 
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Phase 

Mnemonic 

Symbol in Text 

Definition 

2 

DSDMX 

• max 
6s 

Speed brake (command) rate limit 

2 

DSG 

8S G 

Speed brake guidance command 

- 

DS MIN , DSM A X 

6s min , 6s max 

Speedbrake physical upper (lower) limit. 

2 

DTC 

6t c 

Throttle command 

2 

DTCF 

6t CF 

Throttle command filter output 

2 

DTCMN 

..min 

c 

Throttle command limit, min 

2 

DTCMNB 


BCS version of DTCMN 

2 

DTCMX 

..max 

6t c 

Throttle command limit, max 

2 

DTCP 


Past value of DTC 

2 

DTDMX 

• 6t max- 

Throttle (command) rate limit 

2 

DTDMXB 


BCS version of DTDHX 

2 

DTG 

8t G 

Throttle guidance command 

- 

DTMIN , DTH A X 


Physical throttle lower (upper) limit 

2,4 

DYNP 

q 

Dynamic pressure estimate 

4 

DTNP 1,2,3 

q r q 2 « q 3 

Dynamic pressure switch levels for elevator. 




aileron, roll -jet mode III 

2 

BINT 

e INT 

Elevator integrator output 

2 

EINTI 


Elevator integrator input 

2 

EINTWN 

min 

e INT 

Elevator integrator output limit, min 

2 

El NT MX 

max 

e INT 

Elevator integrator output limit, max 

2 

EINTI 


Elevator integrator node 

4 

EL 1,2, 3, 4 

^1» i 2 * -^3* ^4 

Temporary blending parameters 

2 

ET 

e T 

Elevator trim tracker output 

4 

ETA 

V 

Lateral scheduled gain (late entry) 

2 

ETDB 

e TDB 

Elevator trim tracker input deadband 

2 

ETDMX 

•max 

e T 

Elevator trim tracker rate limit 

2 

ETF 

e TF 

Elevator trim tracker filter output 
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Phase 

Mnemonic 

Symbol in Text 

Definition 

2 

ETMN 

min 

e T 

Elevator trim tracker input limit, min 

2 

ETM X 

e^ ax 

Elevator trim tracker input limit, max 

4 

FI ,F2 

f l* f 2 

Longitudinal blending parameters 

4 

F3,F4 

f 3* f 4 

Lateral hlendinq parameters 

4 

FILTER 

FILTER 

State estimation flag: 0=no filterfuse measured 

values of state) 

1=f ilter 

4 

FKLAT 

K 

lat 

Scheduled estimator gains: lateral channel 

4 

FKLON 

K lon 

Scheduled estimator gains: longitudinal channel 

4 

G 1 ,G2 

Gj» ^2 

Lateral fixed gains 

2 

GAINT 


CAS aileron integral gain 

2 

GAIP 


CAS aileron integrator gain (G AP*G AI NT) 

4 

GAIR 

g a 

Filter gain for air data (=TF/TS) 

4 

GAMl 1 , 12, 

21,2 y ll* r 12* y 21* y 22 

Lateral scheduled gains (late entry) 

4 

GANNA 

y 

Flight path angle 

2 

GAN (I) 


Normal acceleration filter coefficients 

2 

GAP 


Roll rate to aileron CAS loop gain 

2 

GAPM 

. 

Roll rate to aileron SAS gain (working) 

2 

GAPMB 


BCS version of GAPH 

2 

GAPMS (I) 


Gain-switch-sensitive SAS values for GAP N 

2 

G AX 


SAS roll stick gearing gain 

2 

GAXB 


BCS version of G A X 

2 

GAYB (I) 


Lateral acceleration filter coefficients 

2 

GBE(I) 


Sideslip error filter coefficients 

2 

GCE(I) 


CAS error filter coefficients (working) 

2 

GCEA (I) 


Angle-of-attack values for GCE(I) 

2 

GCEC(I) 


CAS values (C* or q) for GCE(I) 
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Phase Mnemonic Symbol in Text 


Definition 


2 GCFQH 

2 GCH 

2 GCTH 

2 GCX 

2 GD AC (I) 

2 GDFC 

2 GDSC 

2 GDTC 

2 GEC 

2 GECA 

2 GECC 

2 GEF 

2 GEFB 

2 GEIC 

2 GEICA 

2 GEICC 

2 GEINT 

2 GEINTA 

2 GEINTC 

2 GEQW 

2 GEQHA 

2 GEQNB 

2 GEQHC 

2 GEQMS (I) 

2 GES 

2 GESB 

2 GET 


CSTAF g feedback gain 
Altitude loop gain (H to C) 

Theta loop gain (TH to C) 

CAS pitch stick gearing gain 

Aileron cmnd filter coef. (A to R crossfead) 

Flap command filter coefficient 

Speed brake to elevator crossfeed filter coef. 

Throttle to elevator crossfeed coefficient 

CAS loop gain 

Angle-of-attack value for GEC 
CAS value (C* or g) for GEC 
Flap to elevator crossfeed gain 
BCS version of GEF 

CAS elevator integration gain (GEC*GEINT) 

Angle-of— attack value for GEIC 

CAS value (C* or g) for GEIC 

CAS elevator integral gain (working) 

Angle-of-attack value for GEINT 

CAS value (C* or g) for GEINT 

Elevator damper loop gain (working) 

Angle-of-attack value for GEQH 

BCS version of GEQH 

CAS value (C* or g) for GEQH 

Gain-switch-sensitive SAS values for GEQH 
Speed brake to elevator crossfeed gain 
BCS version of GES 

Throttle to elevator crossfeed gain 
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Phase 

Mnemonic 

Symbol in Text 

Definition 

2 

GETB 


BCS version of GET 

2 

GETF 


Elevator trim tracker filter coefficient 

2 

GEX 


SAS pitch stick gearing gain 

2 

GEXB 


BCS version of GEX 

2 

GFX 


Flap controller gearing gain 

2 

GFXB 


BCS version of GFX 

2 

GHE (I) 


Altitude error filter coefficients 

4 

GN(J1,GND2 


lateral scheduled gains (late entry) 

2 

GPE(I) 


Roll rate error filter coefficients 

2 

GP HE (I) 


Bank angle error filter coefficients 

2 

GPHPS 


Psi-loop gain (psi to phi) 

2 

GPHT (I) 


Phi-turn filter coefficients 

2 

GPU (I) 


Roll rate feedback filter coefficients 

2 

GPHB (I) 


BCS version of GPM(I) 

2 

GPPH 


Phi-loop gain (phi to p) 

2 

GPSE(I) 


Psi error filter coefficients 

2 

GP X 


CAS roll stick gearing gain 

2 

GQM (I) 


Pitch rate feedback filter coefficients 

2 

GQHB(I) 


BCS version of GQH (I) 

2 

GQH 


Pitch rate washout filter coefficient 

2 

GRAY 


Lateral acceleration feedback gain 

2 

GRB 


Sideslip feedback gain 

2 

GRCPH 


Phi to RC gain, turn coordination 

2 

GRDAC 


Aileron to rudder crossfeed gain 

2 

GRDACB 


BCS version of GRDAC 

2 

GRE (I) 


Taw rate error filter coefficients 

2 

GREB(I) 


BCS version, of GRE (I) 


o H 
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DFCS GLOSSARY (Cont) 


Phase Mnemonic 
2 GF.M 

2 GRRE 

2 GRRW 

2 GRRWA 

2 GRRWB 

2 GRRWS(I) 

2 GRW 

2 GRWB 

2 GRX 

2 GRXAUT 

2 GRXB 

2 GRXSAS 

2 GSDDG 

2 GSD 

2 GSUE(I) 

2 GSX 

2 GSXB 

2 GTDDG 

2 GTHE(I) 

2 GTU 

2 GTOE(I) 

2 GTX 

2 GTXB 

2 GDE(I) 

2 GXA 

2 GXAX 

2 GXE 


Symbol In Text 


Definition 

Yaw rate filter coef. (elevator loops) 

Yaw rate feedback gain 

Yaw rate washout feedback gain (working) 
AUTO version of GRRW 
BCS version of GRRW 

Gain-switch-sensitive SAS values for GRRW 
Yaw rate washout filter coefficient 
BCS version of GRW 
Pedal gearing gain (working) 

RAUT mode value for GRX 

BCS version of GRX 

RSAS mode value for GRX 

TAEH guidance DDG to speed brake gain 

Speed brake velocity loop gain 

Speed brake values for GUE(I) 

Speed brake manual controller gearing 
BCS version of GSX 

TAEH guidance DDG to throttle loop gain 
Theta error filter coefficients 
Throttle velocity loop gain 
Throttle values for GUE (I) 

Throttle manual controller gearing 
BCS version of GTX 

Velocity error filter coefficents (working) 
Roll stick filter coefficient 
Roll stick parabolic shaping gain 
Pitch stick filter coefficient 
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Phase 

Mnemonic 

Symbol in Text 

Definition 

2 

GXEX 


Pitch stick paraholic shaping gain 

2 

GXB 


Pedal filter coefficient 

2 

GXRX 


Pedal parabolic shaping gain 

2 

GXSX 


Speedbrake controller parabolic shaping gain 

2 

GXSXB 


BCS version of GXSX 

2 

GXTX 


Throttle controller parabolic shaping gain 

2 

GXTXB 


BCS version of GXTX 

2 

HE 

h e 

Altitude error 

2 

KEF 

h ef 

Altitude error filter output 

2 

HEN (I) 


Altitude error filter nodes 

2 

HG 

h g 

Altitude guidance command 

2 

HM 

H 

Measured altitude 

2 

HREF 

h ref 

EALT mode reference altitude command 

2 

IABES 


ABES present (1) ; ABES absent (0) 

2 

IACSS 


Roll cntl stk steer active (1) ; inact (0) 

2 

IA GAIN 


AGAIN switch position (1 thru 4) 

2 

IAINT 


Aileron integrator active (1); inactive (0) 

2 

IANGFV 


Gravity corrections: UNITY (1), COSTH (2) , 




COSTH * COSPH (3) 

2 

IAINOK 


Pilot requested values for IAINT 

2 

IATRIM 


Roll stick trim activity: 0, +1 (rgt) , -1 (left) 

2 

IBFDFK 


Beta feedback (1) ; lat. accel. feedback (0) 

2 

ICC1 


CAS command switch *1 options (0=off) : 


1 = EC AS mode or control stick steering 

2 = EATT or EALT node 

3 = THG or HG mode 

4 = CG mode (QG or CSTARG) 

5 = AG mode 
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Phase Mnemonic Symbol in Text 
2 ICC2 


2 ICF 

2 ICSTAB 

2 IECSS 

2 IEGAIN 

2 IE IN OK 

2 IEINT 

2 IETRIM 

2 IGEARX 

2 IHC 

4 ILAT 


X LAT 


IM IM 

INGUID INGUID 


INSENS INSENS 


2 IPC 


Definition 

CAS command switch «2 options (0 = off) : 

1 = theta loop (THREF or TH3 ) 

2 = altitude loop (HREF or HG) 

CAS feedback: Q ( 1 ) , CSTAR(2), ALPHA ( 3) 

CAS mode options: CSTAR(I), RCAH (0) 

Pitch cntl stk steer active (1); inact (C) 
EGAIN switch position (1 thru 4) 
Pilot-requested values for IEINT 
Elevator integrator active (1) ; inactive (0) 
Pitch stick trim activity: 0 , + 1 (f wd) ,-1 (aft) 

Extend gear (1), retract gear (0) 

Altitude command: HREF (1), HG (2) 

Mode flag for entry phase lateral control 

0 = early entry, "reverse" aileron 

1 = late entry, "proverse" aileron 
Medium-sampling-rate index: 

0 = perform calculations this pass 
Flag specifies guidance command set 

neg=comraands not updated this pass 
+/- 1=BANKG, ALPHG specified 
+/- 2=ROLLG,C!G specified 
Flag specifies source of measurements 

0=provided by simulator (hybrid sim.) 
1=sensors or environment (SSFS) 

2=G5N and sensors or environment 
(SSFS) 

Roll rate command options (0=off) : 

1 = ACAS mode or control stick steering 

2 = AATT or AHDG mode 

3 = PHG or PSG mode 

4 = PG mode 
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Phase Mnemonic Symbol in Text 


2 IP HC 


2 IPHM 

2 IPHHC 


2 IPSC 

2 IQTRN 


2 IRCHAI 

2 IRCHEI 

2 IRGAIN 

2 IRTRIM 

2 IRW 

2 IRWP 

IS IS 

ISTAFT ISTART 

8 - 


Definition 

2 = A ATT or AHDG mode 

3 = PHG or PSG mode 

4 = PG mode 

Bank angle command options (0 = off) : 

1 = A ATT mode 

2 = AHDG mode 

3 = PSG mode 

4 = PHG mode 

PRH valid (1), invalid (0) 

Phi source for yaw rate command options: 

1 = zero (beta/r-wash. fdbk) 

2 = PHI (measured) 

3 = PHC (commanded, if available) 

Psi command: PSREF(I), PSG (2) 

Pitch rate turn coordination option: 

1 = off, use 0 feedback 

2 = use Q washout feedback 

3 = use Q=R*TANPHI feedback 

Request change IAINT state (1), no request (0) 
Request change IEIHT state (1), no request (0) 
RG AIN switch position (1 thru 4) 

Pedal trim activity: 0, + 1 (lft) ,-1 (rgt) 

Taw rate washout active (1) , inactive (0) 

Pad load indicator for IRW 
Slow-sampling-rate index: 

0 = perform calculations this pass 
Initialization pass(1), not initialization (0) 



DFCS GLOSSARY (Cont) 


Phase 

Mnemonic 

Symbol in Text 

Definition 

- 

ISTRTR 

ISTRTR 

Request initial ization ( 1 ) , no request(O) 

2 

ITHC 


Theta command: THREF(I), THG (2) 

- 

ITH, ITS 

ITM, ITS 

Number of fast passes per medium (slow) pass 

4 

IXX 

*xx 

Fixed estimator parameter 

4 

IZZ 

I 

zz 

Fixed estimator parameter 

2 

iTAGAI N 


Local past-value copy of IAGAIN 

2 

JEGAIN 


Local past-value copy of IEGAIN 

2 

JGEARX 


Local copy of IGEARX 

2 

JRCHAI 


Local copy of IRCHAI 

2 

JRCHEI 


Local copy of IRCHEI 

2 

JRGAIN 


Local past-value copy of IRGAIN 

4 

KA, KE 

K , K 
a* e 

Aileron (elevator) deflection limit parameter 

4 

KALF1 

K «1 

Elevator trim gain (early entry) 

4 

KALF2 

. v 

Elevator trim gain (late entry) 

4 

KC1.KC2 

KC1, KC2 

Low-dynamic-pressure constant filter gains 

4 

KROLL 


Aileron trim gain 

8 

KY(I) 


Individual y accelerometer gain 

8 

KZ(I) 

K zi 

Individual z accelerometer gain 

4 

LAHDAA ,LAHOAE 

X a* X e 

Switch set to 1 if pitch (lateral) ACPS jets 
are needed 

4 

HA, HE 

m , m 
a* e 

Slope of KA (KE) curve 

2,4 

HACH 


Mach number 

2 

HODA 


Current aileron loop mode: 

1 = ABCS 2 = ASAS 3 = ACAS 
4 = AAUT 5 = A ATT 6 = AHDG 

2 

HODACH 


Indicate aileron loop mode change 


1 * no change 2 = HODAGR change 
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Phase 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 


Mnemonic Symbol in Text 


Definition 


MODAG 

HODAGR 

MODAR 

NODB 

HODBR 

NODE 


NODECH 


NODEG 

NODEGR 

HODER 

NODR 


NODRCH 


MODRG 


HODRGF 


3 = NODAR change 4 = first, pass 
Current aileron guidance mode (N0DA=4) : 

1 = PG 2 = PHG 3 = PSG 
Requested values for HODAG 
Requested values for NODA 
BCS mode (1) , not BCS mode (0) 

Requested values for NODB 
Current elevator loop mode: 

1 = EBCS 2 = ESAS 3 = EGAS 

4 = EAtJT 5 = EATT 6 = E ALT 
Indicate elevator loop mode change 

1 = no change 2 = NODEGR change 
3 = HODER change 4 = first pass 
Current elevator guidance mode (H0DE=4) 
1 = QG 2 = THG 3 = AG 4 = HG 
Requested values for NODEG 
Requested values for NODE 
Current rudder loop mode: 

1 = RBCS 2 = RSA5 
3 = RCAS (=2=FS AS) 4 = RADT 
Indicate rudder loop mode change: 

1 = no change 2 = NODRGR change 
3 = HODRR change 4 = first pass 
Current rudder guidance mode (N0DR=4) : 

1 = yaw rate feedback 

2 = yaw rate washout feedback 
Requested values for NODRG 
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Phase Mnemonic Symbol In Text 

2 MODRR 

2 MODS 

2 MODT 

2 MODO 


2 MODDCH 

2 MODOFP 


2 MODOFT 

2 MODOG 

2 NODOGR 

2 MODOR 

4 MO 

4 P 11,12,21,22 P n . Pl2* P 21* P 22 

2 PC P C 


Definition 

Requested values for MODR 
Current speed brake mode; 

1 = BCSSB 2 = MANSB 3 = IASSB 
4 = AOTOSB 5 = DSGSB 6 = TAEMS 
Current throttle node: 

1 = BCSTH 2 = MANTH 3 = IASTH 

4 = AOTOTH 5 = DTGTH 6 = TAEMTH 
Current speed control mode: 

1 = OBCS 2 = OMAN (IAS-man) 

3 = OASB (IAS-sb) 4 = UAUT 

5 = OATH (IAS-th) 

Indicate speed control loop mode change: 

1 = no change 2 = HODOGR change 
3 * MODOR change 4 = first pass 
Flag for FIIOP subroutine speed control 

1 = MANOAL or AOTO-dir 2 = IAS-hold 

3 = AOTO (OG-guid) 4 = AOTO(TAEH-guid) 
Flag for velocity filter moding: 

0 * OFILT inactive 1 = SB-ufilt 

2 = TH-ufilt 

Current guidance speed control mode: 

(1-10, see speed control flag table) 
Requested values for MODOG 
Requested values for MODO 
Roll jet switch parameter 
Lateral derived gains (late entry) 

Roll rate command 
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Phase 

Mnemonic 

Symbol in Text 

Definition 

8 

PCEGK 

k pa 

Pitch channel error gain 

2 

PC MX 

~ max 
P C 

Roll rate command limit, ♦-maximum 

2 

PE 

P E 

Roll rate error / error filter input 

2 

PEF 

P EF 

Roll rate error filter output 

8 

PECK 

<D 

Pitch error gain 

2 

PEN (I) 


Roll rate error filter nodes 

8 

PEF 

P 

ER 

Limited pitch channel error 

8 

PERINT 

P ERINT 

Weighted integral of PER 

2 

PG 

P G 

Roll rate guidance command 

- 

PHASE 

PHASE 

Flight phase flag (internal copy of PHASER) 

- 

PHASER 

PHASER 

Flight phase flag (external) : 1=rollout phas 


2 

2 

2 

2 

2 

2 

4 

4 

4 


PHILON 


-CL 

’ion 


2=cruising phase 
4=entry phase 
5 = orbit. TVC phase 
f=orbit RCS phase 
7=insertion TVC phase 
8=booster TVC phase 


PHC 


Bank 

angle 

command 

PHCNX 

. max 

0 c 

Bank 

angle 

command limit, ♦-maximum 

PHE 

0 E 

Bank 

angle 

error / error filter input 

PHEF 

^EF 

Bank 

angle 

error filter output 

PHEN (I) 


Bank 

angle 

error filter nodes 

PHG 


Bank 

angle 

guidance command 

PHILAT 

* CL 

TLat 

Estimator 

state transition matrix: lateral 



channel 

PHILIH 

0 lim 

Limited value of bank filter output 


Estimator state transition matrix: longitudinal 

channel 
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Definition 

2 

PHM 


Measured bank angle 

2 

PHREF 

0 REF 

AATT node reference bank angle command 

2 

PHSEI 

* SEL 

Pilot-selected values for PHREF 

2 

PHTF 

0 TF 

Phi-turn filter output and/or 

bank angle for turn-coordination (a *> 

2 

PHTN(l) 


Phi-turn filter nodes 

4 

PI (D 

Ap 

Roll rate filter input 

8 

PLIM 

PLIM 

Pitch error software limit 

2,4 

PH 

P 

Measured roll rate 

2 

PMDIS 

Pdis 

Roll rate for AATT or A HDG reengage after 
ACSS 

2 

PMF 

Pp 

Roll rate filter output 

2 

PMFB 


BCS version of PMF 

2 

PMNB(I) 


BCS version of PMN(i) 

2 

PMN(I) 


Roll rate filter nodes 

4 

PN 

Pn 

Nominal roll rate 

8 

PRATEC 

PRATE C 

Pitch rate command 

8 

PRATEF 

PRATE 

Scaled pitch rate error 

8 

PRATEM 

PRATE M 

Blended pitch rate 

8 

PRGG(I) 

Kql 

Individual pitch rate gyro gain 

8 

PRGK 


Pitch rate gain 

2 

PSE 

y e 

Heading error 

2 

PSEF j 

\Ir 

EF 

Heading error filter output 

2 

PSEN(I) 


Heading error filter nodes 

2 

PSG ! 

i 

xj, 

G 

Heading guidance command 

2 

PS M 

i * 

Measured heading 

2 

PSREF 

%EF 

AHDG mode reference heading command 
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Definition 

2 

PSSEL 

to 

*SEL 

Pilot-selected values for PS REF 

2 

QCT 

q CT 

Pitch turn coordination q-command 

2 

QCTMX 

max 

q CT 

QCT maximum (0. = min) 

2 

QG 

q G 

Pitch rate guidance command 

4 

QI (1) 

Aq 

Pitch rate filter input 

2,4 

QM 

q 

Measured pitch rate 

2 

QMDIS 

/v 

q DIS 

Pitch rate for EATT or EALT reengage after 




ECSS 

2 

QMF 

q F 

Pitch rate filter output 

2 

QMFB 


BCS version of QMF 

2 

QtlN(T) 


Pitch rate filter nodes 

2 

QMNB (I) 


BCS version of QMN (I) 

4 

QN 

q N 

Nominal pitch rate 

2 

QWF 

q WF 

Pitch rate washout filter output 

2 

QW1 


Pitch rate washout filter node 

4 

R 1 

r l 

Anqle-of-attack error limits 

4 

R8 

r 8 

Roll angle error limits 

4 

R 11,12,21 

* 22 r l 1 * r 12» r 21* 

Lateral derived gains (late entry) 

8 

RAGTVK 

K 0 

Roll error gain 

2 

RC 

r c 

Yaw rate command 

4 



Derivative of bank angle 

4 

RDOTC 

0 C (0 C ) 

Bank rate command 

2 

RE 

r E 

Yaw rate error 

2 

REF 

r EF 

Yaw rate error filter output 

2 

REN (I) 


Yaw rate error filter nodes 

2 

RENB (I) 


BCS version of REN (2) 

8 

RER 

r er 

Limited roll error for TVC 
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DFCS GLOSSARY (Cont) 


Phase 

Mnemonic 

Symbol in Text 

2 

RG 

r G 

2,4 

RfTO 


4 

RI (D 

Ar , 

8 

RLIM 

RLIM 

2,4 

HM 

r 

2 

RMF 

r F 

4 

RN 

r N 

2,4 

ROLL 

0 

2 

ROLLG 

*G ] 

4 

RPHI 1,2,3 

V %2* r 

4 

RR 

— rel 

8 

RRATEF 

RRATE F 

8 

RRATEB 

rrate m 

8 

RRGG(I) 

K Pi 

8 

RRGTVK 

K p 

2 

RWFB 


2 

RW1 B 


2 

RHF 

r WF 

2 

RW1 


4 

S 

S 

4 

SAT 

slnoirp 

4 

SIGMA A,B , R 

V > % 

8 

TAOINT 

T 

2 

TDSDMX 


2 

TDTDMX 


2 

TEMP 


2,4 

TEMPER 

^ t a 

2,4 

TEMPER 

TK a 


Definition 

Yaw rate guidance command 
Reference air density 
Yaw rate filter input 
Roll error software limit 
Measured yaw rate 

Yaw rate filter output (elevator loops) 

Nominal yaw rate 

Computed roll angle 

Guidance commanded roll angle 

Parameters for maneuver rate calculation 

Position vector in reference coordinates 

Scaled roll rate for TVC 

Blended roll rate signals 

Individual roll rate gyro gain 

TVC roll rate gain 

BCS version of RHF 

BCS version of RK1 

Yaw rate washout filter output 

Yaw rate washout filter node 

Pixed estimator parameter 

Sin (AMEAN) 

Gains for jet logic - pitch, roll, yaw 

Integrator gain 

TF * speed brake rate limit 

TP * throttle rate limit 

Temporary local variable (s) 

Air temperature, deg F and dimensionless 
Input air temperature in deg. Kelvin 
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DFCS GLOSSARY (Cont) 


Phase 

Mnemonic 

Symbol in Text 

Definition 

2 

TETDMX 


TF * elevator trim tracker rate limit 

- 

TF 

t f 

Fast sampling interval 

2 

TF2 


TF / 2. 

2 

TFB 


BCS version of TF2 

2 

TGSDDG 


TF * TAEM/S3 loop gain 

2 

TGTDDG 


TF * TAEM/TH loop gain 

2 

THE 

®E 

Theta error 

2 

THEF 

®EF 

Theta error filter output 

2 

THEN (I) 


Theta error filter nodes 

2 

THG 

9 G 

Theta guidance command 

2 

THM 

? 

Measured theta 

2 

THREF 

0 REF 

EATT mode reference pitch attitude command 

2 

THSEL 

0 SEL 

Pilot-selected values for THREF 

- 

TP! 

t m 

Medium sampling interval 

- 

TNEXT 

TNEXT 

Next DFCS entry time 

- 

TS 

T S 

Slow sampling interval 

2 

TXADMX 


TF * roll stick rate limit 

2 

TXATR 


TF * roll stick trim rate 

2 

TXATRB 


BCS version of TXATR 

2 

TXEDMX 


TF * pitch stick rate limit 

2 

TXETR 


TF * pitch stick trim rate 

2 

TXETPB 


BCS version of TXETR 

2 

TXRDMX 


TF * pedal rate limit 

2 

TXPTR 


TF * pedal trim rate 

2 

TXRTRB 


BCS version of TXRTR 

2 

TXTDMX 


TF * throttle controller rate limit 

2 

T2GAIP 


TF2 * roll integrator gain 
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DFCS GLOSSARY (Cont) 


Phase 

Mnemonic 

Symbol in Text 

Definition 

2 

T2GEIC 


TF2 * pitch integrator gain 

4 

U1 ,02 

U l* U 2 

Roll jet logic gains 

2 

0E 

U E 

Velocity error 

2 

OEF 

U EF 

Velocity error filter output 

2 

UEN (I) 


Velocity error filter nodes 

2 

0G 

U G 

Velocity guidance command 

4 

UJ (1 ) , (2 ) , 

(3) V u y* U z 

ACPS torque commands 

4 

UL(I) 

- L 

Vector in longitudinal plane normal to velocity 

2 

on 

u 

Measured velocity 

2 

UREF 

U REF 

IAS-hold mode reference velocity command 

2 

OS EL 

U SEL 

Pilot-selected values for OREF 

2,4 

OX V,OYV,OZ V 

uv x , uv y , uv z 

Unit vectors along the velocity coordinate 
system, in reference coordinates 

4 

VBODY 

V 

-body 

Air-relative velocity in body coordinates 

4 

VMACH1 

M l 
— rel 

Lateral switch parameter 

4 

VR 

Air-relative velocity in reference coordinates 

4 

HI ,H2 

Oil* u 2 

Gains for ailerbn control logic 

8 

WB(1,I) 

RRATE. 

Individual roll rate qyro signal 

8 

HB (2 , I) 

prate l 

Individual pitch rate gyro signal 

8 

HB (3,1) 

yrate l 

Individual yaw rate gyro signal 

4 

*w> 

HD 

, “d 

Gain for elevator control 

2 

X A 

xa 

Local variable for XAM processing 

2 

XAB 

xa B 

Roll stick detent bias 

2 

XADB 

xa DB 

Roll stick deadband 

2 

XADHX 

• max 
xa 

Roll stick rate limit 

2 

XAF 

' Xa F 

XA filter output 

2 

X AH 

/V 

xa 

Measured roll stick deflection 
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DFCS GLOSSARY (Cont) 


Phase 

Mnemonic 

Symbol in Text 

2 

XAMDIS 

xa DIS 

2 

XAMX 

max 

xa 

2 

XAP 

xap 

2 

XAT 

xa T 

2 

XATB 


2 

XATMX 

max 

xap 

2 

XATMXB 


2 

XATR 

xa T 

2 

XATRB 


2 

XE 

xe 

2 

XEB 

xe B 

2 

XEDB 

Xe DB 

2 

XEDMX 

xe max 

2 

XEF 

xe F 

2 

XEM 

xe 

2 

XEMDIS 

xe DIS 

2 

XEMN 

min 

xe 

2 

XEMX 

max 

xe 

2 

XEP 

xep 

2 

XET 

xep 

2 

XETB 


. 2 

XETMN 

min 

xe T 

2 

XETMNB 


2 

XETMX 

max 

xe T 

2 

XETMXB 


2 

XETR 

xe T 

2 

XETRB 



Definition 

Roll stick disengage level for ACSS 
Roll stick limits (♦-max) 

Processed roll stick 

Roll stick trim 

BCS version of XAT 

Roll stick trim limit (♦-max) 

BCS version of XATMX 

Roll stick trim rate 

BCS version of XATR 

Local variable for XEM processing 

Pitch stick detent bias 

Pitch stick deadband 

Pitch stick rate limit 

XE filter output 

Measured pitch stick deflection 

Pitch stick disengage level for ECSS 

Pitch stick minimum (aft) 

Pitch stick maximum (fwd) 

Processed pitch stick 

Pitch stick trim 

BCS version of XET 

Pitch stick trim minimum (aft) 

BCS version of XETHR 

Pitch stick trim maximum (fwd) 

BCS version of XETMX 
Pitch stick trim rate 
BCS version of XETR 
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DFCS GLOSSARY (Cont) 


Phase 

Mnemonic 

Symbol Ln Text 

Definition 

2 

XFB 


Flap controller bias 

2 

XFM 

xt 

Measured flap controller deflection 

2 

XFMNN 

^min 

Flap controller minimum 

4 

XI 1 


Scheduled control gain: longitudinal 

4 

XI 2,3,4, 5,6 

$2* ^3* ^4* ^5* ^6 

Scheduled control gains: lateral , early entry 

2 

XR 

xr 

Local pedal variable 

2 

XRB 

xr b 

Pedal detent bias 

2 

XRDB 

Xr DB 

Pedal deadband 

2 

XRDMX 

xr max 

Pedal rate limit 

2 

XRF 

xr F 

XR filter output 

2 

XRM 

r*j 

xr 

Measured pedal deflection 

2 

XR nx 

max 

xr 

Pedal limits (♦-max) 

2 

XR P 

xr p 

Processed pedal 

2 

XRT 

xr T 

Pedal trim 

2 

XRTB 


BCS version of XRT 

2 

XRTMX 

xr™ ax 

Pedal trim limits (♦-max) 

2 

XRTMXB 


BCS version of XRTMX 

2 

XRTR 

xr T 

Pedal trim rate 

2 

XRTRB 


BCS version of XRTR 

2 

XS 

xs 

Local variable for XSM processing 

2 

XSH 

xs 

Measured speed brake controller deflection 

2 

XSMDIS 

XS DIS 

XSM disengage level for IftS-sb 

2 

XSMMN 

~ min 
xs 

XSM minimum (retract) 

2 

XSMREF 

XS REF 

XSM reference for IAS-sb disengage 

2 

XSP 

XSp 

Processed XSM 

2 

XT 

xt 

Local variable for XTM processing 

2 

XTM 

xt 

Measured throttle controller deflection 
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DFCS GLOSSARY (Cont) 


Phase 

Mnemonic 

Symbol in Text 

2 

XTMDIS 

xt DIS 

2 

XTMHN 

x ~ t min 

2 

XT HREF 

xt 

REF 

2 

XTP 

xt p 

2 

XTPB 


2 

XTPP 


2 

XTPPB 


8 

YACCER 

• • 

y f 

8 

YACIN (I) 

• • 

Y. 

L 

8 

YACSIG 

• • 

y m 

8 

YAGK 

k ay 

8 

YCEGK 

k ya 

8 

YEGK 

K * 

8 

YER 

Y 

ER 

8 

YERINT 

y erint 

8 

YLIM 

YLIM 

8 

YRATEC 

YRATE C 

8 

YRATEF 

YRATE„ 

F 

8 

YRATEH 

yrate m 

8 

YRGG(I) 

k rl 

8 

YRGK 

*r 

8 

ZACCER 

• • 

Z F 

8 

ZACD 

• • 

Z D 

8 

ZAC IN (I) 

• • 

z. 

L 

8 

ZACSIG 

•• 

Z M 

8 

ZAGK 

k az 

4 

ZETAD 


4 

ZETA 1,2 

Ci* €2 


Definition 

XTM disenqage level for IAS-th 
XTH minimum (idle) 

XTfi reference for IAS-th disenqage 

Processed XTM 

BCS version of XTP 

Past value of XTP 

BCS version of XTPP 

Scaled y acceleration error 

Individual y accelerometer signal 

Blended y accelerometer signals 

Y acceleration gain 

Yaw channel error gain 

Yaw error gain 
Limited yaw channel error 
Weighted integral of YER 
Yaw error software limit. 

Yaw rate command 
Scaled yaw rate error 
Blended rate gyro signals 
Individual yaw rate gyro gain 
Yaw rate gain 

Scaled z acceleration error 
Z acceleration command 
Individual z accelerometer signal 
Blended z accelerometer signals 
2 acceleration gain 
Gain for elevator control 
Gains for aileron control logic 
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SECTION 1 


AUTOMATIC TRACKING 


The designs presented in this section were part of the Apollo/Skylab 
flight programs and are reprinted here with only minor modification. They 
are a baseline for the Space Shuttle pointing routines. 

The automatic tracking function applies to any mission situation 
requiring the orbiter to track a designated body (or target vehicle) by 
pointing a specified body fixed vector at the target or to any mission situation 
requiring the orbiter to rotate at a constant rate around a specified axis , 

Typically, automatic tracking can be utilized for rendezvous navigation, local 
vertical tracking, landmark tracking, and passive thermal control. In this mode 
of operation, the autopilot (DAP) is utilized as a rate command system. 

Attitude control, if required, is maintained by the Tracking Attitude Routine (R61) 
with attitude feedback from the navigation functions. R61 operates at a variable 
computation period, At c , which is generally much longer than the autopilot 
sampling period of 0.1 seconds. Initial alignment of the vehicle to the tracking 
attitude is generally performed by R61 using the automatic maneuver routine (R60). 
In the following TV denotes the target vehicle and ORB denotes the orbiter vehicle. 

1 . 1 T racking Options 

There are five automatic tracking options as defined below: 

Option Purpose 

0 Point a specified spacecraft vector at the 
TV. Do not constrain rotational attitude 
about this pointing vector. (Use VECPOINT) 

1 Point a specified spacecraft vector at a 
celestial body. Do not constrain attitude 
about this vector (Use VECPOINT) 

2 Rotate at a specified rate about a specified 
spacecraft vector. 

4 Same as Option 0, but constraining attitude 
about the pointing vector to a specified 
value (OMICRON). 

5 Same as Option 1, but constraining attitude about 
pointing vector to a specified value (OMICRON). 
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1 . 2 Option 0: Orient the Spacecraft Such that a Specified Body- Fixed Point- 
ing Vector Aligns with the Estimated Line of Sight to the TV. 


The pointing vector (SCAXIS) is defined by two angles: 

1) A specified positive angle (y) from the spacecraft 

X axis about the spacecraft Z axis (0 = 7 < 360 deg) 

2) A specified positive angle (p ) about the Y'-axis 
produced in the XY plane as a result of y( -90 = p = 
90 deg). 
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The LOS is determined by advancing the estimated ORB and TV state 
vectors to the current time. Once the coincidence of SCAXIS and LOS has been 
established, R61 establishes a spacecraft angular rate to maintain the coin- 
cidence. For Option 0, R61 provides the RCS DAP with the following quantities: 


15 ^CA 


desired LOS rate plus LOS correction 
rate in control- axis coordinates 


2) IGAj-j, MGAq, OGAq DAP reference angles, 

i,e., center of the dead- 
band. 

3) AGA a vector defining the desired incremental 

changes in the IMU gimbal angles every 
0 . 1 second. 


The desired LOS rate is computed as follows: 

-LOS = -W ~-C 

-correction =([NBSM] SCAXIS * [REFSMMAT] UNIT (U LOS »/At c 

where NBSM uses the DAP reference angles THETADX, Y, and Z for the trans- 
formation from navigation- base to stable-member coordinates. 

If I - correction! exceeds MAXRATE, distribute as follows: 


^correction = [UNIT (u, correction ) ] 


-LOS 


[REFSMMAT] | UNIT ( u LQS ) x < Y w - g c > | + “correction 


-LOS 


-CA = t NBCA l [ SMNB] <J LOS ' 


where is the LOS rate in stable member coordinates. [SMNB] is the trans- 

formation from stable member axes to navigation base coordinates. [NBSM] is 
the transformation from navigation base coordinates to stable member coordinates. 
[REFSMMAT] is the transformation from reference inertial coordinates to 
stable member axes. [NBCA] is the matrix for transforming a vector from 
navigation base to control-axis coordinates, , Yyy > T^-. > and are the position and 
velocity vectors of the TV and ORB in basic-reference coordinates. SCAXIS is 
the specified axis (7, p ) . MAXRATE is the maximum correction rate and At^ is 
a conservative estimate of the R61 computation interval (thus preventing overshoot 
in vehicle attitude by the time of the next pass). 
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The vector AGA is computed as follows: 


AGA = 


cos IGA sec MGA 
-cos IGA tan MGA 
sin IGA 


0 

1 

0 


-sin IGA sec MGA 
sin IGA tan MGA 
cos IGA 


(0.1 a 


LOS 


where 

IGA = inner gimbal angle 
MGA = middle gimbal angle 
OGA = outer gimbal angle 

where (sJ^qq is the LOS rate in stable-member coordinates computed previously. 

1 . 3 Option 1: Orient the Spacecraft Such tha t SCA XIS Aligns with the LOS to a 
Specified Celestial Body. 

For Option 1, the LOS is determined by advancing the estimated ORB state 
vector to the current time . Once the coincidence to SCAXIS and LOS has been 
established, R61 establishes a spacecraft angular rate to maintain the coincidence. 
R61 quantities provided to the RCS DAP under tracking- option 1 are as described 
for Option 0 with the following exceptions: 

-LOS = * ine of sigh* direction, including aberration correction 


and. 


if the target is a star or planet. 


u, 


LOS ^correction' 


if the target is the center of the earth, 

C 


-los =[REFSMMAT] j UNIT <rc )x( -c ) 


-c 


+ w 

-correction 


1.4 Option 2: Rotate the Spacecraft About SCAXIS 

Option 2 is effected by a special routine known as the Automatic Rotation 
Routine (R67). SCAXIS is as defined in Option 0. R67 initializes COF , BRATE , 
[DEL] , [MIS] , and ADB in KALCMANU and then executes the KALCMANU coding 
to issue the following quantities to the RCS DAP: 
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1) UJ 


CA 


desired rotation rate 


2) IGAq, MGAq, OGA^, .... DAP reference angles, 

i.e., center of deadband. 


3) AGA .... a vector defining the desired incremental 

changes in the IMU gimbal every 0.1 seconds. 


These calculations are repeated once per second. 


1 . 5 Option 4: Orient the Spacecraft as in Option 0, but Further Specify an 

Angle (OMICRON) of Rotation About SCAXIS 

SCAXIS is as defined in Option 0; OMICRON is the positive-sense angle 
between planes A and B - where plane A is defined by the negative orbital angular 
momentum vector (-H) and the LOS, and plane B is defined by the spacecraft 
+ Y-axis and the LOS. 

Alignment of the LOS with either the spacecraft Y axis or the momentum 
vector is proscribed since the necessary two-plane condition for specifying 
OMICRON would not exist. 


OMICRON 


PLANE A 



For Option 4, R61 provides the following quantities to the RCS DAP: 

1) . . . . desired LOS rate in control-axis coordinates 

2) IGAq, MGAq, OGAq, .... DAP reference angles, i.e., center of the 

deadband. 

3) AGA .... a vector defining the desired incremental changes in the IMU 

gimbal angles every 0.1 second. 

The desired LOS rate jin control-axis coordinates is obtained as 

follows: 


-LOS " - W '-C 


= [REFSMMAT ] 

"LOS 


UNIT <aLos )x( 2w-2 C > 


u 

-LOS 


(J CA = [NBCA] [SMNB] u LOg , 


Where is the LOS rate in stable member cooridnates; r^, v^., r^, and v^, 
are the position and velocity vectors of the TV and ORB in basic reference coordi- 


nates . 


The vector (AGA) is computed as follows: 

cos IGA sec MGA 0 -sin IGA sec MGA 
AGA = | -cos IGA tan MGA 1 
sin IGA 0 

+ —correction. 


sin IGA tan MGA 
cos IGA 


(°. 1) u LQS 


where 




correction 


9 DX -°° A D 

V - iga d 

0 DZ - MGA D 


0 . 1 

At„ 


^DX * ^DY* ®DZ des i red gimbal angles :s 

At^, = conservative estimate of R61 computation interval 


The desired gimbal angles O^x* ®DZ are computed as follows: 

First, four unit vectors are computed! 

_s A = cos (OMICRON) unit (u LO g X [v c X £ c ] ) 


+ sin (OMICRON) unit 


11 [—LOS X (- 


LOS X -C X 


)] 


s' A = unit ( SC AXIS X ^q) 
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-B~ -LOS 
-B = SCAXIS 

Given these unit vectors, the unit vectors x, y, and z are calculated defining 
the desired directions of the x, y, and z axes of the navigation base with respect 
to the Basic Reference Coordinate System. The unit vectors x, y, and z and the vectors 

-SM’ ^SM’ ^SM from the matrix REFSMMAT are then used to compute the desired 
gimbal angles 0 DX , 0 Dy , & DZ - 

1 . 6 Option 5: Orient the Spacecraft as in Option 1, but Further Specify an Angle 
(OMICRON) of Rotation About SCAXIS . 

SCAXIS is as defined in Option 0; 

OMICRON is as defined in Option 4. 

For Option 5, R6 quantities provided to the RCS DAF are as in Option 4, 
except as follows: 


and. 


-LOS = line of sight direction including aberration correction 


if the target is a star or planet 
“LOS * 0; 

if the target is the center of the earth. 


-LOS ' [REFSMMAT] 


UNIT (r c ) X (v c ) 
I r r I 
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SECTION 2 


ON-ORBIT ATTITUDE MANEUVER ROUTINE 

by 

John J. Turkovich 


2.1 Introduction 


The purpose of the attitude maneuver routine (ATTITUDE_MANEUVER) is to 
automatically reorient the spacecraft during coasting flight. ATTITUDE MANEUVER 
is an adaptation of the R60 routine from the Apollo flight programs to the Space 
Shuttle with the assumption that a 3-gimbal IMU is used. This implementation is 
coded in HAL and described in accordance with preliminary structured program- 
ming techniques. The general scope of variables and structure of ATTITUDE_ 
MANEUVER is illustrated in Fig. 2-1. Each of the modules and variables used in 
the routine are defined in Tables 2-1 and 2-2 respectively (Subsection 2.7). 


The desired attitude to which the spacecraft is to be aligned can be specified 
in one of two ways: 

1) By a set of terminal desired gimbal angles, j^ c J- 

2) By two unit vectors, and y ANB - 

where is the set of desired (or commanded) gimbal angles defined by 




-ANB * s a body-fixed axis vector, u^, with components in navigation-base- 
coordinates, y^giYj is a desired direction vector. Up, with components in stable- 
member coordinates and y^g-^ specifies the direction in inertial space to which 
-ANB * S t0 k 0 A flag switch, FLAG_3AXIS is provided to specify which 

of the two types of inputs is to be used. 

If the vector type of input is chosen, the routine known as VECPOINT is used 
to compute a set of desired terminal gimbal angles, [0 ], based on the values of 
-DSM’ ancl a ANB’ anc * current gimbal angles, [0]. Regardless of the type of 
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input, KALCMANU is called if the automatic option is selected. Otherwise control 
returns to the executive system. KALCMANU computes a set of parameters (u fl A) 
that determine how the reorientation is to be performed. The maneuver is per- 
formed by rotating the spacecraft about the axis of single equivalent rotation, u^, 
by the angle A.^. The vector u r and the angle A^ are completely determined by 
the initial orientation, [{?], and the desired terminal orientation, [0 ]. Once u r and 
A^ are computed, KALCMANU calls the procedure STEERING. 

The procedure, STEERING, cycles every At c seconds and computes, on each 
pass, five sets of quantities to be used by the digital autopilot in following the angu- 
lar path defined by u r , A. The first of these quantities is a set of autopilot refer- 
ence angles, [0^] (or intermediate desired CDU angles), which is updated regularly 
during the maneuver. To achieve a smoother sequence of commands between suc- 
cessive updates, the program generates a set of incremental CDU angles, 

[A0^],to be added to [0^] by the digital autopilot. The steering routine also com- 
putes the component maneuver rates, (in spacecraft coordinates), a set of bias 
accelerations, g , and a set of bias angles, The bias angles, [|3 ], are 

to be used by the autopilot to prevent fuel- consuming overshoot when starting and 
stopping an automatic maneuver. 

The preceding method allows the autopilot to establish and maintain a space- 
craft rotation about the vector u at a fixed rotation rate w . The maneuvers are 

-r —a 

timed in open-loop fashion so that after a predetermined interval [A0^], and 

^BIAS-^ are set to zero ancl tlie autopilot reference angles, [0^], are set equal to the 
terminal angles, [0 c ] . Thus, upon completion of .the maneuver, the spacecraft 
will finish up in a limit cycle about [0 c J. 

2.2 ATTITUDE JVIANEUVER 

A logic flow diagram of ATTITUDE_MANEUVER is illustrated in Fig. 2-2. 

If the vector option is selected by the calling program, the routine VECPOINT is 
used to compute the terminal desired gimbal angles [0 c ]from u^j^g and Uggjyj* 

If the automatic maneuver is to proceed the automatic mode switch, AUTO_MODE, 
must be set "ON". If it is "ON", KALCMANU is called to calculate the rotation 
vector, u r , and the angle of rotation, A^, and to initiate the cycling of STEERING. 

2.3 VECPOINT 

The desired gimbal angles, [0 ], computed by VECPOINT are determined on 
the basis of minimizing the rotation angle required to align u^^g with Ugg^yj* In 
general, this procedure will reduce the fuel required for the maneuver by not 
constraining the orientation of the spacecraft about the pointing axis, defined by 
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Fig. 2-2 ATTITUDE_MANEUVER : Program 
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-DSM* Additionally, the attitude about l?Qg M will be left unchanged. Consequently, 
the crew may specify the attitude about the pointing axis to satisfy additional com- 
munication, visibility, or thermal constraints, if desired. This is done by manually 
reorienting the spacecraft about the pointing axis either prior to or after the maneu- 
ver. If this procedure is done after the automatic maneuver, a final trim will usually 
be required to assure that is aligned with 

To compute the desired gimbal angles, [0j, from y D g M » ^ANB’ and tde cur " 
rent gimbal angles, [0] , several service routines are required. These routines are 
described below. 

2.3.1 ANGLES_TO_DCM 

This procedure converts three gimbal angles, [0],to a direction-cosine matrix, 

$ 

C, that relates the corresponding spacecraft navigation-base orientation to the 
stable-member frame. The formulas for this conversion are: 

c c c - 

1. 1 '“1,2 1,3 

c c c 

2.1 ^2,2 '“2, 3 

C 3, 1 C 3, 2 C 3, 3 _ 

cos 9 . cos 9 
l m 

-cos 9 . sin 9 cos 9 + sin 9 . sin 9 
l m o l o 

cos 9 . sin 0 sin 9 + sin 9 . cos 9 
i m o l o 

sin 9 

m 

cos 9 cos 9 

m o 

-cos 0 sin 9 
m o 

-sin 9 . cos 9 
i m 

sin 9 . sin 9 cos 9 + cos 9 . sin 9 
i m o l o 

-sin 9 . sin 9 sin 9 + cos 9 . cos 9 

i m o i o 


* 

C 


'“l.l 

C l, 2 

C l, 3 

C 2, 1 
C2, 2 

C 

u 2, 3 
C 3, 1 
C 3, 2 


'3, 3 
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where 


6 q = outer gimbal angle 
0^ = inner gimbal angle 

= middle gimbal angle 

Also see Fig. 2-3 for the logic flow diagram of ANGLES_TO_DCM . 

The interpretation of this matrix is as follows: If a , a , a represent the 

x y z 

components of a vector in navigation -base axes, then the components of the same 

vector in stable -member axes (b , b , b ) are 

x y z 

* 

= C(0) 




2.3.2 DC M_T Q_CDU ANGLES 

This procedure extracts the CDU angles from a direction-cosine matrix, <5, 
that relates the navigation-base coordinates to the stable-member coordinates. 
The formulas for this conversion are 


0 - arcsin C„ , 

m 2,1 

9 = arcsin (-C„ ,/cos 9) 

i 3, x m 

[if i is negative, is replaced by ir SGN (0.) - 0. J 

9 = arcsin (-C„ Q / cos 9 ) 

o 2,3 m 

[if C 2 2 is negative, 9 q is replaced by ir SGN (9 q ) - 6 q j 


The logic flow diagram is illustrated in Fig. 2-4. 


2.3.3 ROT ATION_M AT RIX 

* 

This procedure computes the direction- cosine matrix, R, that relates one co- 
ordinate frame to another that is rotated with respect to the first by an angle A 

about a unit vector u . The formula for this matrix is 
-r 

R (u^. A) = I cos A + u r u T r (1-cos A) + UX r sin A 
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Fig. 2-3 ANGLES_TO_DCM : Procedure ([0]) Assign (C) 
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where 


i 


1 0 0 
0 10 
0 0 1 



2 

u 

rx 

u u 
rx ry 

u u 
rx rz 

T _ 

u u 
-r - r 

u ry u rx 

2 

u 

ry 

u u 
ry rz 


u rz u rx 

u u 
rz ry 

2 

u rz 


— 


— 


— 

_ 



0 

-u u 

rx ry 


* 




UX = 
r 

u rz 

0 -u 

rx 



-u 

ry 

u 0 

rx 



u - unit rotation vector resolved into spacecraft 
axes 

A = rotation angle 


The logic flow diagram is illustrated in Fig. 2-5. 

The interpretation of R is as follows: If a x , a^, a z represent the components 
of a vector in the rotated frame, then the components of the same vector in the 
original frame (b^, b^, b z ) are 


b 


a 

X 

❖ 

X 

b 

= R (u r . A) 

a 

y 


y 

b 


a 

z 


z 


The process used by VECPOINT in determining the terminal desired angles, 

[<g. is illustrated schematically in Fig. 2-6 and in logic flow form in Fig. 2-7. 

VECPOINT uses two frames of reference: the L frame and the M frame - as well 

as the stable-member (SM) frame. The L frame is simply a base frame used as 

the starting point for subsequent rotations . In VECPOINT computations, the L 

frame is the initial navigation -base frame computed from the initial gimbal angles, 

[#]. The vector u~ is resolved from the stable-member frame into this frame as 
~ u * t 

-DL usin § the transformation The L frame is then rotated so that u^^ 

and are coincident. The rotation vector, u r , is defined by the cross product 
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Fig. 2-5 ROTATION_MATRIX: 

Procedure (u r ,A ) Assign (R) 



Fig. 2-6 The Coordinate Frames and 

Transformations used in VECPOINT 
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between u^g anci -DL' ^ otat ^ on a bout the cross-product vector will assure that 
spacecraft maneuvering will be minimized and that the rotation vector will have no 
component in the u^g or -DL direction. 

If the two vectors, u^g and u^g, are either parallel or antiparallel, a special 
computation of u^ is performed. When u^g and u^g are a l m ost parallel then [$ c J 
is set equal to [flg]. If the two vectors are almost antiparallel, then a 180° rotation 
is required for alignment. In this case, the rotation vector is chosen so that u f 
is perpendicular to u^g and lies in the plane that contains yysM(L) anc * -XL" 
however, u A g happens to be perpendicular to the plane, then any vector in the plane 
may be used for u , In this event the x-axis is chosen. From u and A the rotation 


ML 


by 


matrix, R.,. , between the M- and L- frames is computed. Multiplying R 

JYL Is 

Crsm determines the direction - cosine matrix Cj^gj^ from which the terminal desired 
gimbal angles, [e ], are extracted. 


2.4 KALCMANU Maneuver- Parameters Calculation 


The procedure KALCMANU is designed to generate commands for the auto- 
pilot to reorient the spacecraft from an initial attitude to some desired attitude 
specified by [$ c ]. Although a number of authors have discussed optimal solutions 
to this problem, the task is a formidable one because of the nonlinear nature of 
rigid-body dynamics and the kinematics of rotational motion. Optimal solutions 
for a particular class of spacecraft configurations are treated in Reference 1. In 
this reference, the criterion for optimality is based upon the product of maneuver 

duration and fuel expenditure. Using the same criterion, previous design studies 

(21 

based upon a number of suboptimal approaches' ' have shown that a reasonable 
design compromise is to rotate the spacecraft about a single inertially fixed axis. 
This technique is employed in KALCMANU. 


The process used by KALCMANU in computing the rotation axis, y r , and the 

rotation angle, A^, is illustrated in Fig. 2-8. From the current gimbal angles 

KALCMANU computes the transformation matrix between the initial spacecraft 

* 

axes and the stable-member axes, Cgg^. This procedure also computes the matrix 

relating the desired terminal spacecraft attitude and the stable-member axes, 

* # 

Ccsm ' The transformation between the desired and the initial attitude is R^g, where 

where 


* =jc fj. sjt 

R CB = C BSM C CSM 

R^g can be partitioned into its symmetric and antisymmetric components as follows. 
For the symmetric part. 


* 

R r 


1 ,* 

2 (R CB + 


CB' 
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By comparing this with the equations described in ROTATION_MATRIX it can be 

seen that * * _ 

R g = 1 cos A m + u r u r (1-cos A m ) 

For the antisymmetric part, 

* 1 * * rp 

R A = * 2 ^ R CB “ R CB^ 


i 


or 



— » 


— 


0 

-u 

u 

* 


rz 

ry 

R. = 

u 

0 

-u 

A 

rz 


rx 


-u 

u 

0 


ry 

rx 



sin A 


M 


The maneuver angle can easily be obtained from R^g as 


a M = cos 


-1 


/Trace (R CB > - 


-) 


\ 


The rotation vector, u , is normally extracted from the antisymmetric part 


of that is. 


u = unit (-R. , R. , -R. ) 

3,3 i, 3 a 1,2 


As sin A-jyj approaches zero, however, the unit operation becomes indeterminate 
and an alternative calculation must be performed. If the maneuver angle is less than 
0.25 deg, then no complicated automatic maneuvering is necessary and therefore 
one can simply set the autopilot reference angle, [0^], equal to the terminal desired 
angles [# c ] and exit. This is done by calling HOLD_ATTITUDE . If A^ >170 deg, 
a method employing the symmetric part of R^g is used to determine u r , as illus- 
trated in Fig. 2-8. The magnitudes of the components of u r are easily determined 
as 



To obtain the signs of each component, it is first necessary to determine the 

* 

sign of the largest component, by referring to the values of R^. This is done to 

avoid the problem caused by truncation as sin A^ approaches zero. If A M is equal 

to 180 deg, then the sign of the largest component may be arbitrarily chosen. The 

signs of the remaining components may be determined by examining the off-diagonal 
$ 

terms in Rg. 
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2 . 5 Steering Through a Maneuver 


The KALCMANU parameter calculations described in the preceding subsection 
solve the problem of what angular path to take in order to reorient the spacecraft 
without excessive fuel expenditure. However, the problem of how to steer the space- 
craft along this path must now be considered. One solution is to iteratively solve for 
u r and A^ explicitly as a function of the current spacecraft attitude and to control 
the spacecraft angular velocity, cj, by firing the RCS jets so that co is parallel to u r . 
This approach, however, leads to several problems. First of all, because of the 
deadband requirements of the autopilot, and control inaccuracies, explicit guidance 
can lead to excessive changes in control commands as the spacecraft approaches the 
terminal attitude. This problem of terminal guidance is characteristic of explicit 
guidance schemes unless special consideration is given to the terminal- control 
problem. Secondly, there is a timing problem since the computation of u r and A.^ 
is somewhat lengthy and would require an excessive amount of computer time. 

Instead, it is better to formulate the problem in terms of perturbations from a smooth 
reference trajectory that satisfies the end-point constraints. The method employed 
is to solve the global problem once via KALCMANU in order to establish a reference 
trajectory (defined by u r , A^) and then have the autopilot solve the micro-problem 
of controlling the spacecraft deviations from this reference trajectory, using small- 
angle approximations. This approach not only avoids the terminal-guidance problem 
and the timing problem but also provides "tighter" path control for satisfying other 
boundary conditions that may be placed upon the selection of the reference trajectory. 

In order to determine the required steering interface between KALCMANU 
and the autopilot, it is first necessary to define the nature of the autopilot control 
variables. To simplify the discussion, the problem of rotation about a single axis 
is considered first. There are two state variables of interest; namely, the attitude 
error, 0 , defined by 

9 = 9 - 0, 

e d 

where 

9 = actual vehicle attitude 

9 , = desired vehicle attitude (constant) 
d 

and the angular velocity of the vehicle, w. The "state" equations are 

0 g (t) = « (t) 

and 

“(t) = f (0 w) 
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where f (0 , to) is the control acceleration (applied by the jets) as determined by non- 
linear phase -plane switching logic. The control law is defined so that 0 g and w will 
be nulled to within the control-law deadbands so as to maintain attitude hold about 
the reference attitude, 0^. 

New state variables can be defined as follows: 


and 


where 


*e<« ' «*> - *d<‘> + S BIAS 


<Og(t) = co(t) - w d 


0 

’d (t o 


*d (t) 

= W +( “d )(t ' t o ) 

“d 

= a constant 

BIAS 

= a constant 

) = fl(t Q ) 

= constant 


The state equations for the new variables are 

<j) = CO 


“ = f (<t> ,10 ) 

e e e 


With these control variables, the autopilot will attempt to null the rate error, u> g , 
thus forcing the spacecraft to rotate at a constant rate. Note that at the beginning 
of the maneuver the desired angle 0 d (t Q ) is set equal to the initial angle 0(t Q ). The 
constant, ®qiaS’* S added to the attitude error to prevent overshoot when starting and 
stopping the maneuver, is computed as 


_ “dJV 


BIAS 


2 I a 


where 


a = the control acceleration 
c 

To interpret the effect of refer to Fig. 2- 9a. This figure shows a 

simplified phase-plane logic and a typical phase-point trajectory; The maneuver 
starts at point A with 0 g = 0 and w = 0. With the application of steering commands, 
the phase-point instantaneously jumps to a biased condition at B where <t > g 1 
and<J g = (u^ is negative for the case illustrated). The phase-point drifts to 

point C, where the control jets are turned on to reduce the rate error toward zero 
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at point D. The trajectory then traverses a limit cycle DEFG until the maneuver 
finishes at point H, where the biases anc * are removed. At this point, the 

trajectory jumps to an unbiased condition at point I. Note that the actual rate is 
approximately equal to the desired rate and the phase-point drifts to point J, where 
the jets decelerate the vehicle back into a limit cycle around the desired attitude. 
The situation without angular bias, ^gj^g< is illustrated in Fig. 2-9b. Note the 
fuel- consuming overshoot in this case. The inclusion of the bias term is based on 
the assumption that the initial rates are small and that the desired rates are 
achievable during the maneuver. For the three-dimensional problem, it must also 
be assumed that is small enough so that the small-angle assumptions are not 

violated. 


A bias acceleration, a , is also computed and sent to the autopilot to better 
approximate the switching curves in the phase-plane design. In the formulation of 
the phase -plane design a simplification is made to the equation of motion 

• * _ 1 * _ i >1= 

w(t) = J M - J fro x J«) 

where 

M = the torque supplied by the jets 

J *M = or , the control acceleration 
* 

-J fro x J<o) = a the acceleration due to the 
-m 

non -diagonal moment of 

❖ 

inertia J. 

The bias acceleration, a , is assumed small and neglected to arrive at the state 
equations 


m - e d (t) 


Q n 

e 


0 1 


e 

- e 


0 


- 


= 




+ 


W(t) - w d (t) 


CO 

e 


0 0 

— - 


u 

-e 


a 

- c 


and allow the design of three independent phase-planes. To refine the determination 
of the switching curves, is sent to the autopilot and is used in estimating the net 
acceleration. 

The above steering logic has been implemented in the procedure STEERING 
which is described below. 


2.6 STEERING 

The procedure STEERING computes five sets of steering parameters to be 
sent to the autopilot. These parameters are 


2-17 



Fig. 2 



9 Maneuver Phase-Plane Trajectories 
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1) [0^]. a se * i nterme diate desired gimbal angles 

2) [A0 d ], a set of incremental desired gimbal angles 

3) a Set ° f biaS S iml:)al an gl es 

4) oo 

- ( j, the desired angular velocity vector 

5) a , a set of bias accelerations 

STEERING operates on a cyclic basis with period At c - The parameters[/3g^g] , 
and a m are computed once on the first pass whereas [0 d ]and [A0 d ]are computed 
every At £ seconds until the maneuver is complete. 

A logic flow diagram is illustrated v in Fig. 2-10. The initial gimbal angles 
are read and the FIRST_PASS switch is set "ON". The incremental angle of rotation 
is computed from the requested magnitude of angular velocity and the time increment 
At c> From this incremental rotation angle and the rotation vector u r , an incre- 
mental rotation matrix AR is determined by calling ROTATlON_MATRIX. Twotiming 
parameters are computed. These are the duration of the maneuver, t^. 


t M ’ A M /ai r 


and the time at which to stop the maneuver, t , 


t = t,„ + t , 
s M cycle 

where t cyc j e has been set equal to the current time. 

Cyclic computation of the steering parameters takes place in open -loop fashion 
until the computed time to stop the maneuver, t , is reached (Fig. 2-10). If the 
AUTO_MODE switch is "OFF" then STOFRATE is called which sets the desired 
reference gimbal angles [0^] to the present gimbal angles and nulls the remaining 
four sets of parameters [A# d ], [# BIA g], and o^, as illustrated in Fig. 2-11. 
Control then returns to the executive system. However, if the AUTO_MODE switch 
is "ON" the direction-cosine matrix between the navigation- base and stable-member 
coordinate system^ is determined by 

# * # 
c bsm n+1 = C BSM n ar 

From this matrix, the gimbal angles projected to the time of the next pass, [0^], 
are calculated. 

On the first pass STEERING sets the FIRST_PASS flag to "OFF" and com- 
putes w^, and 2^ as illustrated in Fig. 2-10. The desired angular velocity 

vector, is computed from the requested magnitude of angular velocity, w and the 
rotation vector, u , supplied by KALCMANU. 
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Fig. 2-11 STOPRATE : Procedure 
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r- 


0 BIAS (U) 


w dc(U) 1 “dc(U) 1 / ^"cfU) 1 ’ 

^BIAS. 

= 

0 BIAS(V) 

= 

“dc(V) 1 “dcIV) 1 / ^'"cfV) 1 ) 



0 BIAS (W) 


w dc(W) l“dc(W) 1 / (2,Q, c ( W)l) 


L -J 


L. — 

a c (U)’ a c(VY 

(W) are con f r °I accelerations as determined by 


pilot . 


On every pass STEERING also updates the incremental desired gimbal angles. 


[A0,], by 


[A$r| ] ( [# N ] " ftfast /At r ) 


fast 


where 


[<g- 


desired reference gimbal angles projected 

At seconds in advance 
c 

present desired reference gimbal angles 
period of the STEERING cycle 
sampling period of the autopilot 


These incremental angles are subsequently added to [# d ] by the autopilot at the start 
of each autopilot pass. [# d ]is then updated to [0g]and [#g]is replaced by [# N ]for 
use on the subsequent pass. The time of the next proposed pass, t cycle » is computed 
by adding At c to the present value of t cycle - If ^cycle * s * ess tllan the terrn i nation 
time of the maneuver, t , another pass is performed. Otherwise, HOLD ATTITUDE 

S 

is called to update the desired reference gimbal angles, [# d ], to the commanded 
terminal set of gimbal angles [# c ]and nulls the remaining steering parameters 
OW^BIAsl Uj, and a This updating is illustrated in Fig. 2-12. Control is 
then returned to ATTITUDE_MANEUVER which in turn transfers control to the 
executive system. 
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Fig. 2-12 HOLD ATTITUDE : Procedure 


2-23 










2-24 


Subsection 2.7 


Table 2-1 ATTITUDE MANEUVER MODULE DICTIONARY 


MNEMONIC 

DEFINITION 

CLASS 

DECLARED 

IN 

ENTRY FORMAT 

INPUTS 

REFERENCED 

OUTPUTS 

ASSIGNED 

* 

LEVEL 

ANGLES 

TODCM 

Converts an array of 
angles to a direction- 
cosine matrix. 

PROCEDURE 

ATTITUDE— 

MANEUVER 

ANGLES-TO-DCM : 
PROCEDURE ([0]) 
ASSIGN (C); 

... 

— 

2 

ATTITUDE 

MANEUVER 

Sequences an automatic 
reorientation during the 
ON-ORBIT phase. 

PROGRAM 

COMPOOL 

ATTITUDE-MANEUVER: 

PROGRAM; 

FLAGJ3AXIS 

AUTO_MODE 

... 

1 

DCM TO 
ANGLES 

Converts direction- 
cosine matrix to an 
array of angles 

PROCEDURE 

ATTITUDE^ 

MANEUVER 

DCMJTO_ANGLES: 
PROCEDURE (£) 
ASSIGN ([0]): 

TT 

... 

2 

GIMRAL_ 

ANGLE_READ 

Reads the gimbal angles 

PROCEDURE 

ATTITUDE- 

MANEUVER 

GIM BA L_ANG LE_R E AD : 
PROCEDURE ASSIGN (GA); 

... 

... 

2 

HOLD 

ATTITUDE 

Sends steering parameters 
to the autopilot to hold 
attitude at commanded 
gimbal angles . 

PROCEDURE 

ATTITUDE. 

MANEUVER 

HOLD_ATTITUDE: 

PROCEDURE; 

[ 0 C ] 

e d ■ 

[ e BIA&)’-d’ ? m 

2 

KALCMANU 

Computes angle of rota- 
tion, A_ and axis of 
M* 

single equivalent rotation, 

y r 

PROCEDURE 

ATTITUDE- 

MANEUVER 

KALCMANU 

PROCEDURE; 

[0 c } DEG_TO_RAD 

V * M 
[®b1 c bsm 

2 

ROTATION^ 

MATRIX 

Computes rotation matrix 

ft from u and A,.. 

r M 

PROCEDURE 

ATTITUDE- 

MANEUVER 

ROTATION_MATRIX: 
PROCEDURE (u r . A) 

ASSIGN <R); 

* 

^identity 

... 

2 

STEERING 

Computes and sends to 
autopilot the steering 
parameters every At c 
seconds 

PROCEDURE 

KALCMANU 

STEERING: PROCEDURE ; 

A M' W ™PC' 
*■ 

$ 

t fast' At c C BSM, 

REALSEC, AUTO 
MODE 

[ fl d l [ A0 d J- [ 0 BIAS^ 

-d’ ?m’ [®B> C BSM 

3 

STOPRATE 

Sends steering parameters 
to the autopilot to hold 
attitude at present gimbal 
angles 

PROCEDURE 

ATTITUDE- 

MANEUVER 

STOPRATE: PROCEDURE; 

Present gimbal 
angles obtained from 
GIMBAL ANGLE 
READ 

C*dJ t A V 

I^BXAS^’-d’ ? m 

2 

VECPOINT 

Computes [#Jfrom u ANB 
and -DSM 

PROCEDURE 

ATTITUDE- 

MANEUVER 

VECPOINT : PROCEDU RE ; 

-ANB’ -DSM* -XL’ ” 

[®J 

2 


COMPOOL being on level 0 
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Table 2-2 ON-ORBIT ATTITUDE MANEUVER VARIABLE-CONSTANT DICTIONARY 


MNEMONIC 

ENGINEERING 
TERM, RANGE 

UNITS 

DEFINITION 

DECLARED 

IN 

ASSIGNED 

IN 

REFERENCED 

IN 

TYPE 

ATTRIBUTES 

A 

0 

1 A 

> > 
1 A 

rad 

Formal parameter. 
Angle of rotation about 
Hr- 

VECPOINT 

VECPOINT 

VECPOINT 

SCALAR 

DOUBLE PRECISION, 
STATIC 

A 

A 

0 <A < n 

rad 

Formal parameter. 
Angle of rotation about 

u . 

— r 

Angle of rotation about 
u 

— r 

ROT AT ION .MATRIX 

ROTATION_MATRIX 

rotation__matrix 

SCALAR 

DOUBLE PRECISION, 
STATIC 

A_M 

fe 

VI 

2 

<< 

VI 

o 

rad 

KALCMANU 

KALCMANU 

KALCMANU 

STEERING 

SCALAR 

DOUBLE PRECISION, 
STATIC 

ALPHAM 

a 

— m 

rad/ 2 
sec 

Bias acceleration used 
by the autopilot 

COMPOOL 

STEERING 

STOPRATE 

HOLD.ATTITUDE 

... 

VECTOR (3) 

DOUBLE PRECISION, 
STATIC, LOCKED 

AUTO_MODE 

AUTOCODE 

... 

Automatic mode switch 

COMPOOL 

COMPOOL 

ATTITUDE. 

MANEUVER 

STEERING 

BIT (1) 

... 

[bias_beta] 

[ 0 BIAS] 

rad 

Bias gimbal angles 

COMPOOL 

STEERING 

STOPRATE 

HOLD.ATTITUDE 

... 

ARRAY (3) 
SCALAR 

DOUBLE PRECISION, 
STATIC, LOCKED 

e 

t 

— 

Formal parameter. 
Direction cosine matrix. 

ANGLES_TO_DCM 

ANGLES JTO^DCM 

— 

MATRIX (3, 3) 

DOUBLE PRECISION, 
STATIC 

<* 

<* 

... 

Formal parameter. 
Direction cosine matrix. 

DCMJTO_ANGLES 

— 

DCM.TO.ANGLES 

MATRIX (3, 3) 

DOUBLE PRECISION, 
STATIC 

C.BSM 

* 

c 

BSM 

... 

Direction cosine matrix 
between present vehicle 
axes and stable member 
axes. 

KALCMANU 

KALCMANU 

STEERING 

KALCMANU 

STEERING 

MATRIX (3, 3) 

DOUBLE PRECISION, 
STATIC 

C_CSM 

★ 

C CSM 

— 

Direction cosine matrix 
between commanded 
vehicle axes and stable 
member axes . 

KALCMANU 

KALCMANU 

KA LCMANU 

MATRIX (3, 3) 

DOUBLE PRECISION, 
STATIC 

C_LSM 

* 

C LSM 

... 

Direct cosine matrix 
between L- coordinate 
frame and stable 
member axes. 

VECPOINT 

VECPOINT 

VECPOINT 

MATRIX (3, 3) 

DOUBLE PRECISION. 
STATIC 
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ON-ORBIT ATTITUDE MANEUVER VARIABLE-CONSTANT DICTIONARY (Cont) 


MNEMONIC 

ENGINEERING 
TERM, RANGE 

UNIT 

DEFINITION 

DECLARED 

IN 

ASSIGNED 

IN 

REFERENCED 

IN 

TYPE 

ATTRIBUTES 

C_MSM 

* 

C MSM 

— 

Direction cosine matrix 
between M- coordinate 
frame and stable mem- 
ber axes , 

VECPOINT 

VECPOINT 

VECPOINT 

MATRIX (3, 3) 

DOUBLE PRECISION, 
STATIC 

DEG_TO__RAD 

jt/180 

rad/ . „ 
deg 

Conversion factor 

COMPOOL 

COMPOOL 

KALCMANU 

SCALAR 

DOUBLE PRECISION, 
STATIC CONSTANT 

{delta jtheta_d] 

M 

rad 

Incremental gimbal 
angles . 

COMPOOL 

STEERING 

STOPRATE 

HOLD.ATTITUDE 

— 

ARRAY (3) 
SCALAR 

DOUBLE PRECISION, 
STATIC, LOCKED 

DELTAAJV1 

aa m 

rad 

Incremental rotation 
angle. 

STEERING 

STEERING 

STEERING 

SCALAR 

DOUBLE PRECISION, 
STATIC, LOCKED 

DELTAR 

AR 

— 

Incremental rotation 
matrix. 

STEERING 

STEERING 

STEERING 

MATRIX (3, 3) 

DOUBLE PRECISION, 
STATIC, LOCKED 

DELTAT.C 

At c 

sec 

Period of the cycle 
of STEERING 

COMPOOL 

COMPOOL 

STEERING 

SCALAR 

DOUBLE PRECISION, 
STATIC, CONSTANT 

FIRST_PASS 

FIRST.PASS 

— 

First pass switch. 

STEERING 

STEERING 

STEERING 

BIT (1) 

... 

' FLAG_3AXIS 

FLAG.3AXIS 

... 

Type of input flag. 

COMPOOL 

calling program 

ATTITUDE. 

MANEUVER 

BIT (1) 

... 

GA 

9 

rad 

Gimbal angles read 
from sensors 

GIMBAL_ANGLE_ 

READ 

GIM BA L_ANGLE_ 
READ 

GIMBAL_ANGLE_ 

READ 

ARRAY (3) 
SCALAR 

DOUBLE PRECISION, 
STATIC 

id.matrix 

* 

* identity 

... 

Identity matrix 

COMPOOL 

COMPOOL 

ROT ATION. MAT RIX 

MATRIX (3, 3) 

DOUBLE PRECISION, 
STATIC. CONSTANT 

* 

INERTIA 

★ 

J 

slug-ft^ 

Inertia matrix of the 
vehicle 

COMPOOL 

COMPOOL 

STEERING 

MATRIX (3,3) 

DOUBLE PRECISION, 
STATIC 

[MAGACC 1 

["cl 

rad/ 2 

sec 

Magnitudes of the U. V,W 
control accelerations 

COMPOOL 

calling program 

STEERING 

ARRAY (3) 
SCALAR 

DOUBLE PRECISION, 
STATIC 

MAGNITUDE_W_D 

u> 

r 

rad/ 0rt ^ 
sec - 

Requested magnitude of 
angular velocity 

COMPOOL 

calling program 

STEERING 

SCALAR 

DOUBLE PRECISION, 
STATIC 

MPC 

* 

M 

PC 

... 

Transformation matrix 
from pilot to control 
axes 

COMPOOL 

calling program 

STEERING 

MATRIX (3.3) 

DOUBLE PRECISION, 
STATIC 
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ON-ORBIT ATTITUDE MANEUVER VARIABLE- CONSTANT DICTIONARY (Cont) 


MNEMONIC 

ENGINEERING 
TERM, RANGE 

UNITS 

DEFINITION 

DECLARED 

IN 

ASSIGNED 

IN 

REFERENCED 

IN 

TYPE 

ATTRIBUTES 

PI 

ir 

rad 


COMPOOL 

COMPOOL 

VECPOINT 

DCM_TO_ANGLES 

SCALAR 

DOUBLE PRECISION, 
STATIC 

* 

R 

* 

R 

— 

Formal parameter. 
Rotation matrix. 

ROTATION. MATRIX 

ROTATION.MATRIX 

ROT ATION.M AT RIX 

MATRIX (3, 3) 

DOUBLE PRECISION. 
STATIC 

* 

R_A 

* 

r a 

... 

Antisymmetric com- 
* 

ponent of R^ B - 

KALCMANU 

KALCMANU 

KALCMANU 

MATRIX (3, 3) 

DOUBLE PRECISION. 
STATIC 

R_CB 

r cb 

... 

Rotation matrix be- 
tween commanded and 
initial attitudes . 

KALCMANU 

KALCMANU 

KALCMANU 

MATRIX (3, 3) 

DOUBLE PRECISION. 
STATIC 

R_ML 

"ml 

... 

Rotation matrix between 
the M and L coordinate 
frames. 

KALCMANU 

KALCMANU 

KALCMANU 

MATRIX (3, 3) 

DOUBLE PRECISION, 
STATIC 

* 

R_S 

R *s 

... 

Symmetric component 
of R cb - 

KALCMANU 

KALCMANU 

KALCMANU 

MATRIX (3, 3) 

DOUBLE PRECISION, 
STATIC 

REALSEC 

REALSEC 

computer 

ticks /sec 

Scale factor to syn- 
chronize with computer 
clock. 

COMPOOL 

COMPOOL 

STEERING 

SCALAR 

DOUBLE PRECISION, 
STATIC, CONSTANT 

T_CYCLE 

cycle 

sec 

Time at which each new 
pass of STEERING 
begins . 

STEERING 

STEERING 

STEERING 

SCALAR 

DOUBLE PRECISION, 
STATIC 

T_M 


sec 

Time of maneuver 
duration. 

STEERING 

STEERING 

STEERING 

SCALAR 

DOUBLE PRECISION, 
STATIC 

T_S 

‘s 

sec 

Time of maneuver 
completion . 

STEERING 

STEERING 

STEERING 

SCALAR 

DOUBLE PRECISION, 
STATIC 

TF 

l (ast 

sec 

Period of the autopilot 
cycle . 

COMPOOL 

autopilot 

STEERING 

SCALAR 

DOUBLE PRECISION, 
STATIC 

[THETA] 

[01 

- ir < 9. < ir 

— i~ 

rad 

Formal parameter. 
Gimbal angles. 

ANGLES_TO_DCM 

... 

ANGLES_TO_DCM 

ARRAY (3) 
SCALAR 

DOUBLE PRECISION, 
STATIC 
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MNEMONIC 

ENGINEERING 
TERM, RANGE 

UNITS 

DEFINITION 

DECLARED 

IN 

ASSIGNED 

IN 

REFERENCED 

IN 

TYPE 

ATTRIBUTES 

[theta] 

[9] 

-it <Q. < n 
— i — 

rad 

Formal parameter. 
Gimbal angles. 

DCM_TO_ANGLES 

DCM_TO_ANGLES 

DCM_TO_ANGLES 

ARRAY (3) 
SCALAR 

DOUBLE PRECISION, 
STATIC 

[THETA_B] 

[ fi B ] 

- n < 6— < jt 

- B i- 

rad 

Gimbal angles at the 
start of KALCMANU 

KALCMANU 

KALCMANU 

STEERING 

KALCMANU 

STEERING 

ARRAY (3) 
SCALAR 

DOUBLE PRECISION, 
STATIC 

[THETA_C] 

J*> o 

lA 

=» 

rad 

Commanded terminal 
gimbal angles . 

COMPOOL 

calling program 
or 

VECPOINT 

HOLD-ATTITUDE 

ARRAY (3) 
SCALAR 

DOUBLE PRECISION, 
STATIC 

lTHETA_D] 

[ fl d] 

-7 T < 0 , < 7T 

- d r 

rad 

Intermediate desired 
gimbal angles. 

COMPOOL 

STEERING 

STOPRATE 

HOLD_ATTITUDE 

autopilot 

ARRAY (3) 
SCALAR 

DOUBLE PRECISION, 
STATIC, LOCKED 

[theta_l] 

[ fl L] 

- v < 0 t < ir 

rad 

Gimbal angles at the 
start of VECPOINT 

VECPOINT 

VECPOINT 

VECPOINT 

ARRAY (3) 
SCALAR 

DOUBLE PRECISION. 
STATIC 

[theta_n] 

[°n] 

-» < e Ni < ir 

rad 

Gimbal angles predicted 

At in advance 
c 

STEERING 

STEERING 

STEERING 

ARRAY (3) 
SCALAR 

DOUBLE PRECISION, 
STATIC 

U_AL 

-AL 

ft 

Prescribed unit vector 

u A in the L- coordinate 
—A 

frame 

VECPOINT 

VECPOINT 

VECPOINT 

VECTOR (3) 

DOUBLE PRECISION, 
STATIC 

U_ANB 

-ANB 

ft 

Prescribed unit vector 
u^ in navigation- base 

coordinates . 

COMPOOL 

calling program 

VECPOINT 

VECTOR(3) 

DOUBLE PRECISION, 
STATIC 

U_DL 

-DL 

ft 

Desired unit vector in 
the L-coordinate frame 

VECPOINT 

VECPOINT 

VECPOINT 

VECTOR(3) 

DOUBLE PRECISION, 
STATIC 

U_DSM 

. — DSM 

ft 

Desired unit vector in 
the stable-member co- 
ordinate frame. 

COMPOOL 

calling program 

VECPOINT 

VECTOR(3) 

DOUBLE PRECISION. 
STATIC 
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MNEMONIC 

ENGINEERING 
TERM, RANGE 

UNITS 

DEFINITION 

DECLARED 

IN 

ASSIGNED 

IN 

REFERENCED 

IN 

TYPE 

ATTRIBUTES 

U_R 

H r 

ft 

Single equivalent rota- 
tion vector. Rotation 
through the angle A . 

VECPOINT 

VECPOINT 

VECPOINT 

VECTOR(3) 

DOUBLE PRECISION, 
STATIC 

lLr 

Hr 

ft 

Formal parameter. 
Single equivalent rota- 
tion vector. 

ROTATION_MATRIX 

ROTATION^ MATRIX 

R OT AT ION_M A T RIX 

VECTOR(3) 

DOUBLE PRECISION, 
STATIC 

U_R 

Hr 

ft 

Single equivalent rota- 
tion vector. Rotation 
through the angle A^. 

KA LCMANU 

KALCMANU 

KA LCMANU 

steering 

VECTOR(3) 

DOUBLE PRECISION, 
STATIC 

TLrprime 

h; 

ft 

Intermediate rotation 
vector used in deter- 
mining proper signs 
of u r . 

KALCMANU 

KALCMANU 

KALCMANU 

VECTOR(3) 

DOUBLE PRECISION. 
STATIC 

U_XL 

“XL 

ft 

Unit vector along the 
navigation-based x- 
axis (1 , 0, 0) . 

COMPOOL 

COMPOOL 

VECPOINT 

VECTOR(3) 

DOUBLE PRECISION, 
STATIC 

U_YSML 

Hysml 

ft 

Unit vector along the 
y-axis of the stable- 
member frame pro- 
jected into the L-frame. 

VECPOINT 

VECPOINT 

VECPOINT 

VECTOR(3) 

DOUBLE PRECISION, 
STATIC, CONSTANT 

UX_R 

UX r 

ft 

Antisymmetric matrix 
whose components are 

taken from u . 

— r 

ROTATION_MATR!X 

ROTATION_MATRIX 

ROTATIONS AT RIX 

MATRIX (3, 3) 

DOUBLE PRECISION, 
STATIC 

W_D 

-d 

rad / 

sec 

Desired angular velocity 
vector in navigation- 
based coordinate. 

COM POOL 

STEERING 

STOPRATE 

HOLD_ATTITUDE 

autopilot 

VECTOR(3) 

DOUBLE PRECISION, 
STATIC, LOCKED 

W_DC 

-dc 

rad / 

' sec 

Desired angular velocity 
vector in control co- 
ordinates 

STEERING 

STEERING 

STEERING 

VECTORS) 

DOUBLE PRECISION, 
STATIC, LOCKED 
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